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SUNNMOr 


Znv«ati 9 at«d in this study ars finits-slsmsnt nodsling altsrnstivss as 
%islX as ths utility and limitations of ths 2-*D structural rssponss conputsr 
cods CIVM-vJET 4B for prsdicting ths transisnt largs-dsflsction slastic-plastic 
structural rssponsss of 2-0 bsam and/or ring structurss which ars subjsctsd to 
rigid-fragmsnt imptict* Thsss Mttsrs wsrs assssssd in a ssqusncs of conputsr 
runs and by comparing predictions with sxpsrimsntal transisnt-rssponss and 
permanent deformation data obtained for 6061-T651 aluminum beams with both 
ends clamped and subjected to perpendicular impact at midspan by a solid 
steel sphere. Investigated next was the applicability of the CZVM-JET 4B 
analysis and code for the prediction of steel-containment-ring response to 
impact by complex deformable fragments from a tri-hub burst of a TS8 turbine 
rotori corresponding experimental measurements were provided for comparison 
by the Naval Air Propulsion Teat Center (NAPTC) • 

In another aspect of this study, dimensional analysis considerations 
were used in a parametric examination of (a) data from engine rotor burst 
containment experiments conducted by the NAPTC and (b) data from sphere-beam 
impact experiments carried out by the NIT-ASRL. The use of the CIVM-JET 4B 
computer code for making parametric structural response studies on both 
fragment-containment structure and fragment-deflector structure has been 
illustrated. A judicious combination of a few selected experiments plus the 
use of a computer code such as CIVM-JET 4B for parametric and tradeoff studies 
are believed to represent an effective procedure for the design of 2-D 
fragment-containment and/or fragment-deflector structure. 

Since the CIVM-JET 4B computer code as actually implemented accoonodates 
large-deflection, elastic-plastic, transient structural responses but only 
small strain, some modifications to the analysis/computation procedure have 
been developed to alleviate this restriction. Some preliminary exploratory 
results thereby obtained are presented. 


SBCTKW 1 

introdoction 

Kngiaa rotor burst fragau^ts may in^^ct against th 
and/or against special protecUve struc^s These s^ ** 
intended either to contain or to divert the ira be 

escape along a "harmless- path- th ^ ^ 

either fragment containnent or fragm^Hl^ITulT'^Of ^ 

i= thl, pr»ilctio„ Of U,' 

or a«a.«or .truces, .troetur.,, " 

1 .WO.. Purthor, .tt«.tio. .taU b. Z' T 

«t«Uc p»t.ou». con,ider.tio„ of ..atil 

d, ^ «Uti..t.rlol 

If th. a. 

2‘ ““ =»^lo ““ 

““ tb, doflectlo^ Of « 

•trootur. wu b, «o,e»bi.u, th. « .11 loctlo 

«lr«rtion> in this cm., th. a4f.in,tio„ is tm,»j tw^- 
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thickness, and subjected to initial-velocity distributions, prescribed 
eacte m a l ly-applied loads, or fragmex^ is^ct has been developed. The 
capabilities and features of these cooputer codes [1-51* have been summa- 
rized in Ref. 6. Discussed later in this report is the application of two 
of these codes for the analysis of some exaxz^le fragment-structure isapact 
problems. 

For structural response conditions wherein the use of a 2-0 idealiza- 
tion is an excessive over-s impl i f ica^on and where one seeks to predict 
the response in greater detail, the structure needs to be modeled as a 
shell (with stiffeners, if present) [7-9] to permit accounting properly 
for the 3-0 shell structural deflections which are present. On the other 
hand, an excessively fine modeling such as the use of 3-0 solid elements 
to represent a single-layer shell, stiffeners, etc. leads to an excessive 
con^tational burden for many purposes. Hence, "shell-behavior modeling" 
serves as a logical next ir^rovement over 2-0 modeling of C/D structures. 
Accordingly, theoretical prediction methods to cca^)ute the responses of 
plates and shells to initial velocity distributions and prescribed exter- 
nally-applied transient loads [7] are being adapted to predict 3-0 struc- 
tural response to fragment iji^ct [10). Such 3-0 shell response predic- 
tions are discussed no further in this report. 

Confining attention to 2-0 structural response predictions, the 
availability of reliable transient response prediction methods for 
fragment- impacted 2-0 C/D structures may be of use (as noted earlier) for 
preliminary design estimates and parametric calculations. Also, the use 
of such prediction methods may enable one to reduce the scope of actual 
iii^ct/response e3^>eriments on candidate C/D structures; of course, it is 
essential to conduct certain types of esqperiments to establish concrete 
qualitative and quantitative behavior. This experimental information also ' 
can be used to test the accuracy and adequacy of proposed prediction 
methods for certain ranges of parameters and conditions. Properly 
conducted experiments of this type are very expensive; hence, the judicious 
complementary use of theoretical transient response prediction methods is 
of evident value. 

* 

Numbers in square brackets [ ] denote references given in the reference 
list. 
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In order to evaluate the accuracy and adequacy of the cited transient 
structural res^nse prediction methodsr various experiments have been 
carried out to provide appropriate experimental data, Co^lex esqperiments 
involving is^ct of typical engine rotor fragments against conta in me n t 
rings have been conducted at the Kaval Air Propulsion Test Center. Sixtier 
impact eiqperiments with an impacting fragment of simpler and fixed sh^pe have 
been carried out at the NIT Aeroelastic and Structures Researdi Laboratory. 

The Naval Air Propulsion Test Center hasj provided data on the 
responses of aluminum and steel containment rings to (1) ixqpact by a 
single TS8 turbine rotor blade axkl (2) tri-hub burst fragment attack from 
a T58 turbine rotor ([11-14]? for exan^le) . In these cases the attacking 
fragment is coo^lex and ux^ergoes a considerable amount of deformation and 
gecxnetry change daring its impact interaction with the contaizment ring. 

Fcrr case (1), theoretical transient response predictions which take into 
account the changing geometry, mass moment of inertia, etc. of a single 
attacking-rotor-blade fragment have been carried out successfully [15]; 
even this simple case required a substantial book-keeping effort. If a 
similar procedure were applied to the T58 tri-hto rotor burst attack case 
(%^ere each fragment consists of a segment of "rigid” rotor disk and 17 
deformable blades), it is clear toat an-impractically-huge amount of book- 
keeping and cooptation would be encountered. Accordingly, each actual 
bladed-disk fragment has been replaced in the computational model by an 
idealized model which may be vistialized as a "hockey puck" of selected 
fixed diameter; frictional or frictionless impact between this idealized 
fragment and the containment ring is accommodated. The choice of an 
a^ropriate size for each idealized fragment involves some study and 
judgment, and is discussed later in this report. 

To remove that judgmental question and to provide cleaner, better- 
defined, and more detailed transient structural response measurements for 
ev 2 iluating the proposed impact- induced-response prediction method, sixtier 
impact e 3 q>eriments on a smaller scale were performed at the MIT-ASRL [16] . 
The fragment selected was a solid steel sphere of one*inch diameter; 
modelixig of this fragment is simple and efficient since to determine its 
space occt^ancy at a given instant in time, it suffices to know its CG 

*Also, see Ref. 37. 
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SECTION 2 


STRDCTUHAI. RESPONSES OP AUKtNCM BEAMS- WITH-OAHPED 
ENOS TO STEEL SPHERE IMPACT 

2.1 Problem DaPln< 

in order to obtain appropriate and detailed two-dinensional (2-0) 
transient structural response data under well-defined impact 
so that a definitive evaluation could be made of the adequacy of the 
approximate collision-interaction analysis employed in the CIVM scheme, 
some simple experiments have been conducted at the MTT-ASRL and are re- 
ported in Ref. 16. Beams of 6061-T651 aluminum with nominal 8-in span, 
1.5-in width, and 0.10-in thickness and with both ends ideally clamiied* 
have each been subjected to midspan perpendicular Impact by a solid steel 
sphere of one-inch diameter as indicated schematically in Pig. 2. impact 
velocities ranged from those sufficient to produce small permanent deflec- 
tion to those needed for threshold rupture of the beam. Spanwlse-oriented 
straxn gages were appUed to both the upper and the lower (inpacted) sur- 
face Of the beam at various midwidth and spanwise locations, in each test, 
transient strain measurements were attempted for 8 of the gages, after 
each test, permanent strain readings were obtalned-from- all surviving gages 
Also, permanent deflection measurements of the beams were made. 

To provide ideally-clamped boundary conditions at the ends of the 
beams, these specimens were machined from a solid block of 6061-T651 alumi- 
num (J.53 X 3.60 X 13.0-in) as indicated schematically in Pig. 3. The indi- 
cated -integral support collar- was in turn bolted securely to a heavy flat- 
ground steel channel support structure with close-clearance 1 / 2 “-13 boIo- 
Krome shoulder screws, this design and arrangement is believed to have pro- 
vided a very close simulation of "ideally-clamped- ends. To reduce the 
hazard of undesired or "premature cracking- at the boundary because of 
stress concentrations at a sharp re-entrant comer, all -inside comers- 

were machined to a radius of 1/8-in which is somewhat greater than the 
nominal beam thickness. 
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Target beans of 6061-T651 aluminum were employed because this materi- 
al has well-knoim mechanical behavior, exhibits very little-straln- harden- 
ing, and is relatively' insensitive to strain-rate effects. A l-inch diameter 
steel sphere was chosen as the iji5>acting fragment since its geometry is well 
defined (and preserved) throughout the impact attack, is devoid of sharp 
edges which might cause atypical local shearing penetration of a target, and 
was of convenient size, shape, and mass for use %d.th the available experi- 
mental launching apparatus so as to achieve the desired levels of fragment 
kinetic energy at initial lapact. The target beam dimensions: thickness 
h - 0.100-in, width w - 1.50-in, and span L - 8.00-in were chosen for vari- 
ous reasons, as follows. A small width w m desired so that the ensuing 
beam response would be of essentially 2-D~ch;racter; however, the aiming ac- 
curacy of the available launching apparatus was limited and, hence, a 1.5-in 
width was selected to permit achieving impacts reasonably well centeredat 
the beam's midwidth-midspan location. A length-to-thlckness ratio of about 
80 was chosen so that the beam model would possess a significant degree of 
bending as well as extensional rigidity. Further, the span was kept large 
compared with the diameter of the attacking fragment so that "local irpact 
effects- would be confined in a region remote from the supported ends of 
the beam. Finally, the 8-in span was also selected-such thaf the time to 
achieve peak deformation would be short enough to prevent the required com- 
puting time for the theoretical transient response prediction schene from 
being excessive, this judgment was based upon much previous computational 
experience [17, for example) on similar structures. With these selected 
dimensions for the "missile" and the target beam, it was anticipated that 
the beam would exhibit essentially 2-D deformation everywhere except in 
the immediate vicinity of the "sphere-beam impact point". Thus, this ex- 
perimental information could serve the dual purposes of providing data to 
check (a) beam 2-D response predictions and (b) narrow plate 3-D struc- 
tural response calculations (in the future), to fragment ia?»ct. The 

reader is invited to consult Ref. 16 for further details of these experi- 
ments. 

An inspection of each beam specimen indicates that except near the 
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^ and tests CB-13, CB-16, and CB-18 are employed, and the pertinent 
data are as follows. 


Specimen 

and 

Test 

CB-13 

C3-18 

C3-16 


Beam Dimensions (in) 
_Thic]c.p h Width, b Span. I. 


Steel-Sphere Data 


1.501 

1.498 

1.495 


8.002 

8.002 

8.002 


Weight 

(grans) 

66.738 

66.810 

66.809 


Initial Impact 
-Velocity (In/see) 


« IP deflection while speclnen 

-18 displayed a state of large peznanent defomatlon, model CB-16 experi- 
enced an even greater response: to Just b«ely beyond threshold rupture for 
the entire cross section. 

Befor, ntt-ptin, to cOrt»r. prediction, with erperi.«ntnl nea.ore- 

oT«' “ 1* 

of finit, osoi to i»d.i the been, thi. i, di.coe.ed 
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e».itivitp Of the b... ..terirt the predicted transient response, of 
steel-sphrte i.p»ct«l benrt CB-13 CB-18 ... lUostr.ted. Ih«, 1„ 

Sobs«n;ion 3.4 cort»ris«.s,.re n«l. betrten stroctural response predic- 
tion. nnd «.p,ri..ntrt ncurenent. for specinsn. CB-13, CB-18, and CB-18. 
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pradlctlons mqt— reasonably with each other. 

Froa the CB~1 calculations » one can conclude that the use of a 20— 
element per half-span FB nodel with 4 DOF/node elsMnts will provide an 
adequate structural aodel fcr preiJicting large-deflection, elastic-plastic 
transient structural responses of this iapulsively-loaded bean structure 
insofar as the nuober of finite eleaenu used to aodel the structure is 
concerned. Bence, it is r eason able to es^ect that a similar 
fineness should suffice for analysing the responses of siaUar to 

steel-sphere impact — provided that the iapact-interaction and associ- 
ated modeling requirements ^u:e not more stringent. 

It should be noted at this point (as discussed also in Appendix A) 
that converged finite element predictions by the assumed displacement 
method depend not only iq>on using an adequately large nmOyer of elements 
but also iq>on esq>Ioying a type of element-assumed-displacement field 
which includes continuity of w and its first spanwlse derivative and con- 
tinuity of the V displacement as well as constant strain terms and rigid 
body modes. These conditions are fulfilled by the present beam elements 
as noted in Appendix A. Sought here, however, is engineering convergence 
(accuracy to perhaps 5 to 10 per cent with the use of a modest of 

rather than strict mathematical convergence (accuracy to 
a fraction of a per cent by using a great number of elements); is 

"achieved- in the present instance by using about 20 or more elements 
per beam half span. 

Next, consider the analysis of beams CB-13 and CB-18 which were 
subjected to steel-sphere impact. To analyze these cases, various niaibers 
of both uniforvlength and nonunlform-length finite elmnents were con- 
sidered for modeling the entire nominal 8-in span of the beam, as indi- 
cated in Fig. 5. Shown for only the half span are the finlre-element 
arrangements examined; the length of each element is shown. For the 
several arrangraents of elements of nonuniform lengths, the intent was to 
use fine meshing near the point of initial iapact and near the clamped 
ends to provide detaUed information where significant bending (as well as 
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staretching) was expected; in regions of saaU spatial variati<ms of behavior 
and where very little bending was expected* coarser oeshing (longer eleaients) 
were expected to provide adequate. inf oraation. Botrever, the presence of 
snail-length finite eleaents. cansM the iini»timtm natural frequency of the 
nathenatical nodel of the entire structure to be nudli higher would be 
the case if one were to use the sane miaher of unifom-length finite eleaents. 
Thus* since the central-difference tiaewise finite-difference operator is 
used for the tinewise solution of the equations of motion in CIVM-JET 4B* 
this meus in turn that one must use a much smaller *<■■«* increment At for a 
nonunifom mesh FB model than for a corresponding uniform-mesh PB model. 
Further* if one were to analyze a contaiment-ring structure that could be 
impacted by* perhaps* several fragments at any ciromiferential locat io n of the 
ring* it would be infeasible and illogical to attespt to model the ring with 
nonuniform finite eleaents — an array of unifom-length finite elanents would 
be clearly the sensible choice. Therefore* in view of thess considerations, 
subsequent finite element modelings and calculations of the steel-sphere- 
is^acted beams utilized only uniform-length finite elements. 


In partic ular as Fig 5 shows* the entire beam was modeled with 
either 23 or 43 unifom-length finite elements, with the midspan of the 
center finite element located for convenience at the midspan of the beam 
itself: the location of the point of initial steel-sphere impact. The 
effects of these t»«> FE modelings upon the predicted* responses of speci- 
mens CB-13 and CB-18 will be examined next. From the earlier "convergence 
discussion" with reference to inpulsively-loaded model C3-1, it is ex- 
pected that the present 43 element model (Z 22 elements per half span) 
will exhibit "converged impact-induced response". It is also of in- 
terest to assess the sensitivity of the predicted iapact-induced re- 
sponse to a similar but coarser FB modeling; namely, essentially one-half 
of the number of uniform-length elements used in the former case. 


fi^^st* it should be noted that, as e;q>ected, the structural be- 
havior of these steel-slq>ere impacted beams was two-dimensional at all 
spanwise stations except for those near (within about 0.8-in of) the 
center of initial irpact; nearer the point of initial inpact, very 
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pronounced three-dimensional structural deformations are evident. There- 
fore, since the CIVM-JET 4B analysis and code pertains strictly to 2-0 
structural response > theoretics-experimental canparisons (see Subsection 
2.4) can be made logicSly oSy at stations |x| > i.o in. However, as 
a matter of curiosity, conparisons between experiment and 2-0 predictions 
will be made in both of these re^ons: |x| > 1.0 in and |x| < 1.0 in. 

CIVH-JET 4B calculations were carried out to predict the in?>act- 
induced responses of specimens CB-13 and CB-18. In each case the impact- 
ing fragment was modeled ^ a rigid 2-0 circular fragment of 0.5-in 
radius, m the impact-interaction portion of the computer code, locally 
perfectly-elastic in?)act was assumed; hence, a coefficient of restitution 
e of 1.0 was used. The Simiinum beam material (p = 0.098 ib/in^ weight 
density) was modeled as beiiaving in an EL-SH fashion by "illustrative 
Fit A". The time increment sire At used was 1.0 and 0.5 usee for the 23- 
element and the 43-element model, respectively. 

Shown in Fig. 6 are the predicted time, histories of the lateral 
displacement w for the 23- and the 43-element case at spanwise stations 
X i 1.2, 1.9, and 3.6 in for specimen CB-13. A similar plot is given in 
Fig. 7 for the more severely iji?)acted specimen CB-18. Both of these 
figures indicate very close agreement between the predicted time his- 
tories of displacements from these two modelings _at_each- of these span- 
wise stations. 

A more sensitive and meaningful comparison of these predictions 
involves axial strains on the upper (0) or the lower (L, or in?)act) 
surface of the beam at various spanwise stations. Such coa?)arisons are 
shown in Figs.. 8 and 9 for specimens CB-13 and CB-18, respectively. 
Summarized in the following tabulation are the comparisons of these two 
predictions (23 and 43 element, both EL-SH) with each other (and/or versus 
experimental data for rough con^arison only) shown in these figures. 



Time Histories of Strain on 
Surface 0 (upper) and/or L( lower) 


Specimen 

Fig. 

Station x(in) 

Predicted 

Measured 

CB-13 

8a 

0 

0* and L* 

- 

It 

8b 

0.6 

0 

0 

H 

8c 

1.2 

0^ 

- 

It 

8d 

1.5 . 

0 and L 

0 and L 

it 

8e 

2.1 

0* 

- 

It 

8f 

3.7 

0 and L 

0 

CB-18 

9a 

0 

0^ and 


It 

9b 

0.6 

0 

0 

n 

9c 

1.2 


0 

It 

9d 

2.1 

0* 

- 

It 

9e 

3.7 

0 and L 

0 


Superscript "a" denotes locations which coincide with the midspan 
station of a finite element.. Superscript "h" pertains to stations x 
which coincide with a nodal station of the PE ioodel; the associated 
strain value plotted is the average value given by the two elements at 
that nodal junction station. Predicted strains at locations riot marked 
by "a** or "b" occur at some station intermediate between the end and the 
midspan of a finite element. 

The inclusion of measured str^Lxn data in these plots is intended 
to provide some intxiitive guidance in assessing the 23-element vs. the 
43-element predictions. 

Figiares 8a and 8b are included for academic interest; at these loca- 
tions (x » 0 and x » 0.6 in), the computer values are legitimate for 2-D 
structural response and hence may be compared with each other. However, 
it is not proper to con^are predictions with measurements at x = 0.6 
since 3-D ^structural behavior is evident in the experiments and the mea- 
sured strain history there must be affected accordingly; nevertheless, 
the ** inappropriate" 2— D predr ctrons are included vs. measurements in 
Fig. 8b. Strains predicted by nodal averaging at station x =* 1.2 in cire 
shown in Fig. 8c for both FE modelings. At x = 1.50 in, the predicted 


13 





strains are at element locations intermediate between midelement and a node; 
in Pig. 8d, 2-0 predictions are cooqpared with measured strains also under 
2-D structural response conditions. Midelement strains for both calcula- 
tions are shown in Pig. 8e for x - 2.1 in; again one observes similar 
qualitative and quantitative behavior between these two predictions. How- 
ever, no measurements are available to aid one’s assessment. Pinally, 

Pig. 8f shows predicted and measured strains in a region idiere the in^r- 
tance of the bending behavior is changing rapidly with spanwise location; ; 
predictions at element locations intermediate between the midelement 
point and a node are "conqtared” with measured values on the upper surface 
* ” 3.70 in. Here the 23-element prediction differs markedly from the 
43-element result; incidentally, the latter is in better but not really 
good agreement with experiment. Near x - 3.70 in the strains change rapid- 
ly with X in this region — as shown in Pig. 8g where the predicted upper- 
surface strain vs. x at time after initial impact TAIl ^ 1020 ysee is 
plotted for the 23- and the 43-element calculation. Pigure 8g demonstrates 
the severe strain gradients which occur near the "midspan iiiq>act region" 

And n6flur clAinped end of tdie beesi. 

More extensive conq)arisons between El-SH predictions of transient 
strains are given in Pig. 9 for the more severe impact 
P^®'^^3ing for specimen CB— 18. The midelement strains predicted at 
X ■ 0 again exhibit pronounced differences between the 23-element 
and the 43-element case. At x - 0.6 in, predicted upper-surface strains 
are shown and pertain to element locations between midelanent and a 
node; also shown are measured strains on the upper surface at x « 0.6 in. 
Again, the 43-element predictions agree better with experiment than do 
the 23-element predictions, although these 2-0 predictions are not 
stricUy comparable with the 3-D measured data. However, at x - 1.20 in, 

2-0 predictions could be compared legitimately with measurements since 
the latter pertain to 2-0 behavior also at that station, except for the 
jrA-lustrative stress-str ain Pit A ; here, better and racner reason- 
able (perhaps coincidental) agreement between the 43-element predictions 
of nodal-averaged strain and experiment is seen. Midelement predicted 

^Comparisons of predictions with measured straps may be found in Pig. 17. 
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strain in the 2-D region is plotted in Fig. 9d at x > 2.1 in. Finally, 
near tne claiaped end at x - 3.70 in. predicted upper surface stralnr for 
CB-18 are con^ared with measureoents; the severe strain gradients in 
this region can readily lead to pronounced theoretlcal-ei^terinental 
discrepancies. The co»?>arison here is slMlar to that noted earlier con- 
cerning specloen CB-13. 

With respect to representing rapid variations of strain with span- 
wise distance at stations very close to the clamped ends, these variations 
can be represented adequately in principle by eaploying a finer finite 
element mesh in those regions. In addition, however, it should be noted 
that: near the clas^d ends of severely-impacted or loaded specimens there 

is visual (orange peel) evidence of the occurrence of transverse shear defor- 
maUen ~ whereas the present finite elements do not include this behavior. 

If transverse shear deformation effects become sufficiently Important, finite 
elements which include this effect could be used but this use. in turn, will 
lead to a much larger highest frequency of the mathematically-modelfed 
structure 117] thereby requiring one to use a much smaller time step size 
At if the central-difference operator were used for the timewise solution. 

It was noted earlier that near the midspan Isqpact point, the 43- 
element predictions differ from the 23-element predictions for each of the 
specimens CB-13 and CB-18. The principal reasons for this involve the de- 
tails of the (approximate) impact-interaction claculaUon method employed, as 
described in the following. In essence, when the fragment impacts the beam, 
a calculation is made to determine the amount of momentum transferred from 
the fragment to a portion of the beam centered about that impact point — 
called the "locally-affected region" which is defined as that which will 
fflqperience an imparted velocity increment. The distance from the impact 
poi.nt to either end of this locally-affected region* is denoted by the 

"effective length" one estimate** is that is the product of the 

simple elastic rod-wave velocity (|) ^^^ 

* ■ • 

^See, for example. Appendix A of Ref. 4 or Appendix B of Ref. 5. 

Another escinate of is discussed later in this subsection. 
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used, where E is the elastic ^^us p Is the mas. per unit eolune of the 
bean nuiterlal, for aluninon ^ 

carter code was en^loyed and this progran uses the central-difference tlae- 
wise solution operator, the tine step sire At nust be less than 2/o) to 
avoid calculation divergence, where is the naxlnun frequency oJ^e 
mathenatical nodel of the structure for _snall-displacenent linear elastic 
behavior. For large-deflection elastic-plastic transient response calcula- 
tions. experience has shown that one nust use a snaller value for At such as 
At < 0.8(2/«^^). Shown in the following tabulation are n . 2/« 

°*®^^^“nax’ ' actajal tinestep size At used for the 23^ and dSelanent 

-odellngs of^C8-13 and CB-18, Included also are the finite elenent length 


and I, 


'eff 




Mo* of d} 

(^nifotn 

Elco-;nt* frad/see^ 

Timm Stop Slzm Vmlimm 

^mtf 

Wt) ^ 
0 UMd 

fin) 

SlM. 

- Longth 

“?’«x “ ^ Mtl 

«aji moM 

! JK22S1 (USCCJ (usae^ 

0.134349x10^ 
« 0.333600x10^ 

1.608 1,287 1 

0.688 0.50 

0*1985 

0.0993 

1 12sL_ 

0*3478 

0.1860 


Note that since (1) in^^ct occurs at the center of the finite elenent which 
is centered at nidspan and (2) > one-half of an elenent length.the 

locally-affected inpact region Includes not only the center elenent but 
also a portion of one elenent on each side of the Inpacted elenent. ihese 
three elenents. therefore, contribute to the nass portion of the bean which 
acquires a distributed increnent of velocity associated with the monentun 
transfer. The nass involved and enq>loyed in the conq>utatlonal logic is 
that at the end nodes of the central elenent for the aasatwhi f;lnite 
elenent nodel and hence represents the nass of 2 con?>lete elenents: all of 
the center elenent and one half of the elenent on each side of the center 
element. Thus, it can be seen readily that for a given impact, the larger 
impact-affected region assumed to be associated with the 23-element nodel 
will experience a considerably smaller velocity increment than will the 
impact-affected region in the 43-elenent nodel; thus, near the Impact 
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d.., ^ 

thM for tt. 23-,l«.„t I»).l (.„, tor mmol,, ri,,. 9. .od 8b for 

»n^*? ’* **’ »P«cll»» C8-18). «s l„ St. vuimt's 

«f r..po«., ,t loctloM r«ot. trc th. -Ir^ot- 

^.ctM „ d. socd for U.O fs- .0. ,s^i«tt cfoulo- 

on. «,«opt oo«r a. cl«,»d ad rt.ro prt«a.c«l badio, „o„t, 
produces significant differences. 


• «ti«te for the iapact-affected half length 

*"eff “ * single iapact can be seen to yield unreasonable 

estimates if one were to employ a veiy large or a very small value for 
At. in the former case, the normal-<Urection momentum transfer produced 
by the impact-induced through-the-thickness (and spreading) stress waves 
will have ceased long before the large At (which may be permitted by 
another timewise finite difference operator) has elapsed; hence, the actu- 
al region receiving this normal-direction momentum transfer or velocity 
increment will be much smaller than the above estimate for t „ assumes 
On the other hand, for very small At. the impact-induced strLs wave 
passing through the thiOmess may not yet have reached the opposite free 
surface - thus, toe momentum transfer envisioned in the calculation model 
could not have been completed physically. Thus, for a single inpact, one 
can make a more rational estimate for from stress wave propagation 
considerations even though the basic structural response and impact/inter- 
action model being used does not take explicit account of this type of 
"across-the-thickness" stress wave propagation behavior. 


Figure 10 depicts in a sinplified fashion two idealized types of 
impacting-fragment geometries: ( 1 ) a 2-0 fragment of circular cross section, 
which contacts the target at a "line" and (2) a 2-D fragment of rectangu- 
lar cross section. For impact of toe 2-0 beam target by the solid cylin- 
drical fragment, a highly simplified picture is shown of the sequence of 
stress waves which propagate "across the thickness direction" of the beam 
For present purposes let it be assumed that these stress waves are simple 
elastic waves. At the initial impact point, a "line source of compression" 
IS produced; a compression wave prSpagates into the plate as a cylindrical 
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front from this "line source" and weakens as it e^^ands because of 

"geonetric divergence". The material particles being traversed by the^ 

coiQ>ression wave are given a material velocity in the direction of the 

propagating wave front. When this conpression front reaches the opposite 

free surfeure of the plats tdtere the stress is required to remain at the 

zero-stress condition, a cancelling tension wave arises and can be viewed 

as emanating from source center SC-2 located at a distance h be^nd this 

free surface. This reflected tension wave also weakens as it propagates 

(and diverges) through the plate and in^arts an addition particle velocity 

to the material in a direction opposite to the direction of travel of the 

propagating tension wave front. When this tension wave front reaches a 

free surface, a corresponding reflected compression wave is generated at 

pseudo-source-center SC- 3* When this new con^ression front reaches the 

opposite free surface, a new tension front is created with pseudo source 

center at SC-4. The wave front strength created initially by each new 

pseudo so\irce becomes weaker and weaker for each successive source — 

roughly inversely proportional to the square of the distance from the 

source center. Hence, by the time that wave front 4T from SC-4 reaches 

the originally impacted surface, one can aurgue reasonably that siabsequent 

waves would be able to induce only slight additional increments to the 

particle velocity of the material traversed by the stress waves. Thus, 

a reasonable estimate of the beam length region which experiences a 

substantial velocity increment in the direction normal to the plate is 

selected as the intersection of wave front 4T with the half thickness 

station of the beam? from this geometric picture one can show that 

i 1.94 h (see Fig. 10a) . The actual stress wave propagation is 

much more conq>lex than described here (see Refs. 22 and 23, for example); 

however, it is believed that the present description is adequate for the 

purposes of selecting a reasonable value for the beam inpact-aiffected 

zone half-length this type of in 5 )act, therefore, it would not 

be unreasonable to choose 1.5h < L _ < 2.0h — independent of the calcu- 

err 

lation rime step size At. 

By similar arguments, the beam impact-affected zone half length 
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^eit impacted by the rectangular section fragnent depicted in 

Pig. 10b Bight be estiaated «* i^/2 > h < t ^^ < £ / ♦ 2h. 

With respect to the 43releBent calculation discussed earlier, the 
present estimate for uould be approximately 1.94h i 0.194 in, hence, 
the total impact-affected zone of i 0.386 In eould auunt to about 

teo beam element lengths (i.e., 2 (.186) - .372 in). This turns out to be 
consistent with the previously-discussed 43-element calculation wherein 
the employed resulted in estimating the impact-affected zone to in- 
clude both nodes of the ia^wcted center element of the asseii4>led beaib 
structtjre. Those two nodes are regarded as receiving the impact-iii?«rted 
Impulse and do represent the mass for 2 beam elements. Thus, the 43^ 
element calculation can be regarded as valid from the collision-inter- 
action standpoint, whereas the 23-element calculation assumes an un- 
reasonably large region of the beam to receive the momentum transferred 
by the attacking fragment in a single impact. This means that near the 
point of intact, one should use shorter elements than in the 23 equal- 
length-element case in order to take into account properly the beam 
region which effectively receives the transferred momentum (and associ- 
ated velocity increment) from a single intact. 

Finally, it is of Interest to note that 43-element EL-SH Pit A 
calculations were carried out for the steel-sphere-impacted CB-18 beam 
by using various At values in the CIVM-JET 4B program. For this modeling 
it was found that 2/w^ - 0.860 usee; hence, the rule-of-thumb choice 
for an acceptable At to avoid calculation instability from 
error growth for large-deflection elastic-plastic transient response 
predictions is 0 . 8 ( 2 / 0 )^) - 0.688 usee. Calculations were carried out 
for At values of (a) 0.5 usee ^ich corresponds to .58(2 / m^^ ) and (b) 

0.75 usee which corresponds to -87(2/0)^). For At - 0.5 the transi- 

ent response calculations were found to be well behaved and converged. 
However, at At = 0.75 usee, the predicted impact- induced response became 
appreciably different (and larger) than for the converged calculations - 
This experience indicates that when carrying out impact-induced nonlinear 
response calculations which utilize the timewise central difference 


< V 
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operator (CIVM-JEr 4B), one must be careful 

the critical At = 2A) for this «« * At weU below 

terns It aun^ / operator when applied to gy,. 

mbs. It appears from this e*oerlai«,<» ..u -4S2£ sys- 

[17 fo- e«-nsi^i - rough rule>of-thtmib 

U7. for example) for selecting Atd.e.. At < 0.8(2/o. n . 

able and should not be -stretched- unduly. ' »a*^^ ^ «PPlio- 

Strain-Rate Effaet« 

strain rate on the uniaxial stress-strain- behTJor^i^^cm*^ 

Siven by the following Cowper-Symonds relation [ 24 ,: ^ 

‘^ok 'd' ' (2.1, 

where and are, respectively, the static and the ^ • 

dependent yield stresses of the kth 

cal sublayer-e is the straintte -‘-i- 

constants. For aluminum Ref ' ° ^ ^ strain rate 

and p = 4. these va7 use of o - 6500 sec'^ 

for iUu.rr.ti... 

Althou,. th, ,trM,-.tr.i. l,*.rt.r of 606i-I65i ■ 

lleved Ka w-f»a.u . wox To51 aJ.uniinum is h^— 

xaevea to be rather insensitive to strain is be- 

oio».r raicuiati... .... .... ..rriT“r Lr “■ 

O-U- uaio, t.. atrua-rat, 

a-SS-SR „rfi.ti.„ ‘’°=“ •adp-4. n.... 

a«i With experiment for specir^ns CB-13 JZTJ:Zi:ZZ: 

tively, are 43-element 

Placement w versus Ume as well as the 

deflections of each specimen as a funcJ^T Permanent 

- f^aacuo. a.4 
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for corresponding EL-SH calculatiom although small strala-rate sensitivity 
is used, strain-rate effects in metals effectively "stiffen* the structure 
ther^y leading to a lesser response.* Rote that for both the CB-13 and 
the CB-18 example, the EL-SH prediction of the permanent deflection is in 
better agremaent with experiment in the 2-D region than is the EL-SH-SR re- 
sult? the experimental result is plotted along the midwidth station y « 0 

of each beam and includes pronounced 3-D deformation evidencs over the region 
|x| < 0.8 in. 

For specimen 43-element EL-SB and EL-SH-SR were 

made wherein steel-sphere perpendicular impact was specified to occur 
initially at the midspan location of the beam. Theoretically, the en- 
suing beam response and fragment motion should remain synaetric. However, 
because of calculation roundoff error, one finds that after many time in- 
crements in the calculation have occurred some response asyanetry and con- 
sequent fragment motion asymmetry appear.* Hence, a second EL-SH-SR 43- 
element calculation was done wherein, at the end of each At time increment 
in the calculation, the fragment was required or "constrained" to move and 
i;emain in the * “ 0 plane? for convenience, this calculation is termed 
herein 43C. This constraint resulted in only small changes in the pre- 
dicted beam response — affecting slightly the time at »*ich peak strain 
was predicted to occur at a given station only near midspan, but not af- 
fecting the value of that peak significantly. These 43-element EL-SH, 

EL-SH-SR, and EL-SH-SR/43C strain predictions at x - 0 are shown in 
Fig. 13a for specimen CB-13. ^ 


EL-SH and EL-SH-SR predictions of transient strain on the i:«>per 
and the lower surface of beam CB-13 at a representative 2-D location 
(X « 1.50 in) are compared with measurements in Pig. 13b? similar com- 
parisons for the upper surface of the beam at station x - 3.70 in are 


shown in Pig. 13c. At both of these locations, the EL-SH-SR calculation 
provides the better agreement with experiment. Finally, Pig. I 3 d shows 
EL-SH and EL-SH-SR predictions versus measuranents of the permanent 
strain as a function of spanwise location x. In the 2-D response region, 
the measur ed permanent strain is less than about 0.7 per cent for 
*Pigure 13c indicates an "exception" to this for the EL-SH-SR calculation 
but this occurred because of inadvertent asymmetric impact — corrected 
subsequently by the EL-SH-SR/43C calculation. ^ rected 
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spqci awn CB~X3i is sssn tiist tiis predicted pezBssenfc u^per** eorf ace steeins 
tend to be soneMfaat larger than the Deasnred values, with the EXt-SH-SR 
prediction being in closer agreenent with eiCr«^i>eat. 

Siailar cos^arison of BX(>^ and EL-SR-SR strain predictions for K«jin 
CB>18 are shown in Pigs. 14a, 14b, 14c, and 14d. Rote In Pig. 14d that the 
neasured p erm a nent »®per-surface strain vs. s in the 2-0 region is les s 
about 1.5 per cent. Again predict^ pezoanent strains m scaewhat larger 
than neasured. 

Overall, the EZ«-SH-SR predi c tions are better in sene regions 
the EL-SB predictions are better in other regions, particxilarlp rM-r the 
point of i m p a c t. Because some predicted asymetric behavior was observed 
for s^^nnetric iaqpact conditions i^en using the EL-SH-SR condition and 
also because of- uncertainty of proper strain-rate parameter values, most ^ 
of the following congxarisons and studies are perfozmed with only the 
EIi-SB material modeling. 

Pinally, it should be noted that there is widespread agreement that 
the mechanical stress-strain behavior of tes^ered aluainaa is relatively 
insensitive to strain rate at rocm teng>erature conditions (coo^ared with 
mild steel [25] or titaniua [26], for exaBg>le), but the degree of -the 
sensitivity is- in doubt. R efer ence 27 indicates that at room tenqperature 
various of the common almainum alloys in the soft ("O") condition 

(such as 1060-0, 1100-0, 6^1-0, 7075-0, for exan^le) display distinct 
rate sensitivity; these same alloys when tai^>ered to the T3 or T6 condi- 
tion exh i b it almost no strain-rate sensitivi^. There is, however, 
scatter of strain-rate data amongst various reference sources [28-30] . 

2.4 Comparisons of Predictions with Experiment 

Examined in Subsection 2.2 was the effect on the predicted deflec- 
tions aj3d stT 2 d.ns of using various finite-element modelings of the steel- 
sphere-impacted beams when the beam materiaO. was treated as behaving as 
an elastic strain-harding (EL-SH) material. Next, in Subsection 2.3 
the effects of including material £train-rate (SR) dependent behavior 
was e;q>lored (the material was treated as EL-SH-SR) by ai 5 >loying 
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iUttstratiw strain-rate constant values 0 ■ 6500 sec ^ and P ■ 4 
alvniniai [251. In both of those subsections, the predictions were 
cooked priaarily with each other. Based qpon those studies and the 
considerations discussed in those subsections, the coobinatlon of 
(1) a 43-elesient aodel to represent the entire beam and (2) material 
EL-SH b^vior was chosen as being a reasonably satisfactory basis for 
ng steel— sfAiere— ia^acted beam 2— D responses for the pv^pose of 
coB^aring*predictlons with experimental measurements: for beam specimens 
CB-13, CB-16, and CB-18, except that now the stressr^train curve of Ref. 21 
for the 6061-T651 beam material should be employed. ; These comparisons 
(experiment, predictions with illustrative stress-steain Fit A, and pre- 
dictions with the new stress-strain fit for 6061-T651 aluminum) are dis- 
cussed now in this subsection. 


For convenient reference, the test conditions and observed effects 
on these beam specimens are sunmarized concisely in the following: 



Steel Sphere Data 


Specimen 

and 

Test 

Weight 

(grams) 

Initial Impact 
Velocity (in/sec) 

Initial Kinetic 
.Energy (in-lb) 


Post-Test 
Condition 
of Beam 

CB-13 

66.738 

2490 

1187 

Moderate perma- 
nent deforma- 
tion 

CB-18 

66.810 

2795 

1489 

Large perma- 
nent deforma- 
tion 

CB-16 

66.809 

2870 

1569 

Specimen rup- 
tured; thres- 
hold rupture 
condition 


Note that beam CB-13 exhibited a moderate degree of permanent deformation 
while the more severely Impacted CB-18 beam displayed large permanent 
deformation. Within about 0.25 in of the center of impact and along a 
roughly circular bemd, some tiny cracks are observed on the V 5 >per sur- 
face of specimen CB-13 but these cracks do not extend through the 
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thickness of the specioen. Siailar but 

are seen «« * 1 . aoBewhat nore pronounced cracks 

do not extend through the thickness of the wJj»en The 

^-P.ct a» cx,« u„ ^ 

.rrrj^rr jrr“^rr r* — 

the bee- «. «. ^ *»<* edges of 

tiie oeaflie Tbe CB*X6 inmacti ^ 

-uQiBct condition is essentially that for •threih«i^ 

fnc«r.- .!«. tt. ««l ^ ^ 

did «ot p„. tte>w, tt. CTKIc-ptod»*l op«ln,. ..-- 

.« «s«v« in tt. b««idtb dir««ion - iMictl., tt. , 

«. th. lo«« ^ clill^d M „ 

d^~ of itiaiain, to:, ^ W~i«». CB^U to orl 

nounced for specimens CB<-18 and CB-16. 

an examination of these specimens, one can see readily that 

,d-B, ^ tarJ nJi- 

"«ens in a spansise region totalling about l s ^ k • 

-Of tB, o-idspan, atatio. of iTtiT.;^ TTl^ 7“ ™ 

1 , *«Faux: IX a 0). iVlso, Very closa 

clan^>ed exxisp a sznall amount of 3 -d defonBa^ 4 n 
si_ ^ ox J D aefornation occurred. Thus, 

rr:rrir:r:r “ 

-..^to ^ata ait. ,t. ^‘“oTt^.^rraror.'"^ 

large strams ), predictions can be compared "legitimately- with measur 
^ts only in those beam spamrise station regions in which predominantir 

• _ I I , . 0 m. For comparisons made for curiosity at locations 


formulation and^lem^LtS^s^^^^-^SOT^f^JT^®**®^^^® rigoroi 
matter are given in Section 6? “® “™ents on this 
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beyond the "zone of validity", one n»«4. 

Of sod. , ^ «' “• *»»-» ootor. 

'•'-^wisons. Also, the eneriaenfcai - ^- 

•-,*0. Of od»«of. od. . :rr:“ ~jr::; 

txons of the structure. ^ ° 


a«.«o, Of u.. ^ a, ,. 

Of ^ ^ ZZZZ T“‘ 

:t»o prodlction i«tai tt«.„hou. u„ “ “ ~ 

-tdrial toet»i„, oo=«x.O T “ 

f-0 E-icSioni <43 iU^t..Z‘lZ»ZT 

«otf.c oia stroi. V „ .tation. , . o 0 6 14 '«»'«» oppar- 

spoei«n, CB-13, CB-W. dd CB-16 (U.tM L, ttl ^ ’’’ “ 

sovority). Th« odioctod ,. , of u»rM.i,„ 

-.4.0, to. cB-x. 

•Oout 5.4 p«r „«t« for ....i.*, Zu toTf ' 

«nen CB 16 than for speciaen CB-18. 

CoB^arisons of the y,, bean-surfaf.« g, ._■ 

=S2H«»t^*or,p.c3«.,S.33™ 

figures listed concisely beloe for convenient refer^er^'” ^ 




It should be notrt that all of these predictions utilize illustrative 
stress-stra<" Pit A and hence, should not properly be canpared with ex- 
periment; later, predictions utilizi n g a corresponding stress-strain fit 
for the actual 6061-T651 beam material will be presented. 

Since specimens CB-13 and CB-18 each retained its structurzO. in- 
tegrity throughout its transient response, one can compare strain pre- 
dictions legitimately with experiments^. mesis\irements for these speci men s 
for the entire history, including the permanent defor m at i on co n d i t i on for 
the 2-D region O.k < Ixj < 3.70 in. In this regard for transient strains, 
see Figs. 13b and’ 13c for model CB^13, and Figs. 9c and 14b for model 
CB— 18; for both specimens it is seen that the 43-elanent EL-SH prediction 
is in reasonably good agreement with e;^eriment and that the strain 
Igvels involved do not exceed the bounds of validity of the present theo— - 
retical formulation and computer program. However, the bounds of validity 
appear to be exceeded near the midspan im^ct point and close to the 
clamped boundaries. Note also from Fig. 8b for CB-13 and Fig. 9b for 
CB-18 where "illegitimate" comparisons are made in the 3-D experimental 
response region at x » 0.6 in, that the 2-D predictions are in rough 
agreement with measurements during the early part of the time history 
but then tend to diverge vs. experiment as time progresses; also at sta- 
tion X => 0.6 in, the strain gages become detached from the speci:.-=n 
at an unknown time after initial impact and the strain trace on t. e oscil- 
loscope record vanishes. Finally, the spanwise distribution of predicted vs. 
experiment permanent upper stir face strain is shown in Fig. 13d for 
specimen CB-13 and in Fig. 14d for specimen CB-18; in both cases ihe ex- 
perimental strain level is small in the 2-D region 0.8 £ |x| £ 3.70 in, 
and the EL-SH-SR prediction is in better agreement with experiment than 
13 the EL-SH prediction. Also, there is somewhat better theoretical-experi- 
mental permarent strain agreement for specimen CB-18 than for specimen CB-13. 

For specimens C3-13 and CB-18, Figs. 11a and 12a, respectively, 
shew 2-D predictions of the vertical deflection w as a function of time 
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after Initial in^ct at spaxnrise stations x « 0, 1.2# and 3.6 in; 

corresponding transient deflection measurenents were not obtained. Hoi#<- ^ 
ever, measurements were made of the spanwise distributions of permanent 
w-deflection; compared in Figs* ilb and 12b for specimens CB--13 and CB*18, 
respectively, are 2*D predictions of the permanent w-deflection versus 
measurements. These plots of measured w-def lection represent at each span- 
wise station x for |x| ^ 1.0 in, the average of measurements taken along 
lines y » -0.5, 0, and *K).5 inch sixu:e each specimen exhibited some perma- 
nent twist because iiig>act did not occur exactly !at the y » 0 Width loca- 
tion (see Table 7 of Ref. 16). For |x| < 1.0 ixir the measiured %r-deflec- 
tion values plotted represent those along y » 0 and the average of those 
at y » -0.5 in and y » -H)*5 in at each |x| station; in this region, 
the 3-D nature of the deflection is evident from these plots. Note that 
in the 2-D region (0.8 < |x| < 3.7 in), the permanent deflection predic- 
tions (43 eluent. Fit A) coo^>cLre much better with experiment for specimen 
CB-18 than for specimen CB-13; a similar but less dramatic result %ms noted 
earlier with respect to permanent strains. Further, the El-SH-SR w~def lection 
prediction appears in each case to agree slightly better with experiment than 
does the EL-SH prediction. However, this might be accidental because only very 
approximate values for the strain rate parameters D and p in the Eg. 2.1 
approximation were et^Ioyed — the authors are unaware of better values 
for D and p which experimental data could confirm for 6061-T651 aluminum. 

Also, remember that aill of the predictions discussed thus far have 
utilized "illustrative stress-strain Fit A". Shortly, CB-18 predictions 
which en^loy a corresponding stress-strain fit for the data of Ref. 21 
for 6061-T651 aluminum (the beam material) will be discussed. 

With respect to the predictions and comparisons given in Fig. 17 
for specimen CB-16, note that the predictions assumed that CB-16 re- 
mained structurally intact throughout its transient response, whereas 
beam CB-16 fractured across its entire cross section at essentially the 
midspan impact attack station near its transient response peak at a 
time estimated to be between about 550 and 650 microseconds after initial 
impact. After through fracture occurs at midspan, there is rapid 
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unloading of tha •mambrane stata" of tha baan on both halvas of tha baaai 
this is saan aost raadily by comparing tha aaasurad strain traca on tha 
uppar surfaca at stations x ■ 1>S0 in for CB“18 vs» that for CB~16« 
Similarly, tha maasur^ strain traces shown in Pig. 17 for spaciman 
CB-16 depart progressively from predictions after about 600 microseconds 
from initial impact. 

Recall that "illustrative stress-strain Pit A" was made from the 
uniaxial static tensile stress-strain data vs of Raf. 19 by a 
piecewise linear fitting of the Ref. 19 data. Now for 6061-T651 aluminum, 
a similar data fitting has been made cl the data in Pig. 29a of Risf. 21j 
the resulting approximation is teemed Pit B and is defined by the follow- 
ing stress-strain pairs: pai« 0 in/in) » (44,000 psi, 0.0044); 

(46,500 psi, 0.035); and (54,000 psi, 0.175), and serve to identify the 
quantities defining each of the three elastic, parfectly-plastic sublayers 
of the associated mechanical sublayer model. The cited Ref. 21 data and 
the Pit B (as well as the points of Pit A for the Ref. 19 data) are shown 
in Pig. 18. The piecewise linear Fit B was made to tha Ref. 21 data 
shown on Fig. 18 whereas Fit A was made to similar data of Ref. 19 for 
a slightly different aluminum alloy condition. Note that although these 
fits were made to different sets of data, they are in reasonably close 
agreement with each other up to values of about 5 to 6 per cent; be- 
yond this level, these fits differ somewhat but are still. close with 
respect to the accumulated areas under these curves. Hence, one would 
expect the transient response predictions from using the "more realistic" 
Fit B to differ only a small amount from those obtained from using the 
previous "illustrative Fit A" for the conditions of either specimen. 

CB-13 or specimen CB-18. 

Predictions of steel-sphere-impact produced structural response 
for specimen CB-13 were made by using, as before, 43 equal-length 4 DOF/ 
node elements for the entire beam and stress-strain EL-SH Fit B. The 
resulting predictions of w-deflections and strains are compared with the 
earlier EL-3H Fit A predictions and with experimental measurements for 
specimen CB-13 in Fig. 19. As expected, these two predictions are seen 



to differ only slightly from each other? similar qualitative and quantita- 
tive comparisons were n=,*ad at other spanwlse stations (not shown). Hence 
corresponding Pit B calculations were not made for specimen CB-13. 

To Illustrate the effect of a slag>ler approximation to the stress- 
strain behavior of the aluminum material, an elastic, perfectly-plastlc 
(EL-PP) approximation was employed wherein the -yield stress was taken as 
45,000 psl, this approximation, of course. Involves no strain hardening. 
Shown also In Pig. 19 are transient response predictions for specimen 
CB-18 when the EL-PP approximation was used. It Is seen that these pre- 
dlctxons for w (Pig. I9a) and for upper surface strain at x « 1.2 in 
(Pig. 19c) are in very close agreement with the previous Pit A and Pit B 
EL-SH predictions. Only at the midspan iai>act station x • 0 (Pig. I9b) 
does the EL-PP prediction differ significantly from the earlier pre- 
dictions, but here the predicted strains are larger than those for which 
the present analysis and con^)uter code are valid* 

Prom the CB-13 and CB-18 comparisons of 2-D predictions of impact 
induced transient and permanent strain as well as transient and perma- 
nent w-displacement. it is seen that predictions agree better with mea- 
surements for the more severely impacted CB-18 beam than for the less 

severely impacted CB-13 specimen. The former exhibits more pronounced 

3-D deformation, necking, and shear deformation than does the latter; of 
course, none of these effects are included in the 2-0 predictions. Since 
specimen CB-13 behaves in a manner more closely approximating that of the 
idealized 2-D prediction model than does CB-18, one would expect better 
agreement between CB-13 predictions and measurements than between CB-18 
predictions and measurements in regions of 2-0 response behavior, but 
this is opposite to the present findings. 

If one assumes that the 2-0 predictions will tend to overestimate 
the structural response for "gentle impact conditions" such as those of 
specimen CB-13, a similar overestimate might be expected also for 
severely-impacted specimen CB-18. However, in the latter case, the more 
pronounced roles of 3-0 deformation, necking, and transverse shear 
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defonnatloa would tend to produce larger structural deflections than other- 
wise would occur; this tendency, therefore, would bring predictions and 
experisiental measurements of beam response into closer agreement. This is 
in accord with the con^arisons observed. However, the question of tihy the 
predicted 2-D intact-induced structxiral response is too large for cases of 
dominan^y 2-D structural behavior remains to be answered adequately. 

'to analogous situation was observed in the studies of Ref. 20 \^erein 
2-D predictions were carried out for aluminum beams with clanted ends and 
which were "geometrically identical" to the CB-13 and CB-18 beams. Those 
specim^s, termed CB-1 and CB-4, were intulsively loaded [16] by the sheet 
explosive loading technique (SELT) over their entire width and for a total 
spanwise length of 1.8 in centered at midspan; a uniform initial velocity 
of 6660 and 10,590 in/sec for beams CB-1 and CE-4, respectively, was im- 
parted to each beam in the explosive-covered region. The 2-D predictions 
for the more severely loaded CB-4 specimen are in better agreement with 
oxperiment than are those for the less severely iopulsed CB-1 specimen. 

to assessment of the uncertainties associated with (1) the ex- 
plosively- ia^wrted iii 5 >ulse for the CB-1 and CB-4 specircns [16] and (2) 
the pre-impact steel sphere velocity for the CB-13 and CB-18 tests [16] 
indicates that those uncertainties are far too small to contribute sig- 
nificantly to the removal of the discrepancies between 2-D pre^ctions 
and experiment for the "gently loaded" speciaans CB-1 and CB-13. For 
exan?>le, for specimen CB-13, it was estimated that at worst the sphere's 
pre-inq>act velocity might have been 2400 in/sec rather than the 2490 
in/sec veaue cited in Ref. 16 from calibration firings tests. Shown in 
Figs. 20a and 20b are 43— element EL— SH Pit A predictions of transient 
and permanent deflection, respectively, for these two sphere-ia?>act velocity 
vcQues. It is seen for each calculation that the predicted permanent de- 
flections exceed those measured. 

While greater strain-rate sensitivity than assumed and a greater 
energy absorption in the 3-D region near the impact zoart (leaving less to 
be absorbed by the remaining structxire in 2-D deformation) would tend to 
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acwunt £or sob* o£ this overprediction, one cannot cite these as decisive 
factors. The Inclusion of transverse shear defomation and/or the use of 
a larger nuBber of the present finite eleeents would tend to increase the 
severity of the predicted response. Thus, as of now, a decisive quantita- 
tive «- pi «««<•< nn for the cited discrepancies resalns elusive. 

2 e 5 Comnents 

2.5el Utility and- Limitations of the 2-D Analysis 

The prediction method under discussion applies to the large- 
deflection elastic-plastic transient structural responses of siaple struc- 
tures which undergo strictly two-dlBensional (2-D) displaceaents when sub- 
jected to fragaent inpact or to iopulse loading, but only saall to Boder- 
ate strains (< 10 per cent) are accoBBOdated. The present analysis and 
associated coBputer prograBS do not take large strain conditions into 
account (properly) ; it would be tiaely and useful now to extend the analy- 
sis to include large strain behavior but the tiae and effort required to 
acconplish this task is a Batter for future consideration. 

The present 2-D analysis of structural response to fragaent iapact 
pertains to relatively low speeds of iapact in the direction of the noraal 
^^0 surface of the iBpacted 2-D structure. Accordingly , transient 
structural response rather than through-the-thickness severe stress-wave 
response doBinates; hence, under sufficiently rigorous impact attack, the 
structure can undergo severe deformations which results typically in a 
tensile- type fracture near the fragment iapact point for a ductile metal 
target structure [12, 16, for example). At these iapact velocities one 
does 'not observe (local stress-wave-induced) shear-plug failures. wliich 
are conmonly seen for much-higher-speed iapacts [31] . 

The present analysis assumes both the target structure and the im- 
pacting fragment to act in a strictly 2-D manner. Conceptually, this im- 
plies that the structure is a narrow beam or ring and that the attacking 
fragment has the same geometry at all stations in a plane parallel to the 
"spanwise axis" of this impacted structure; in particular, the attacking 


fragment is idealized as consisting of a non-deformable solid circular 
cylinder which extends across the entire width of the beam or ring, 
clearly, this idealized model does not match the physical conditions 
present in the steel-sphere- ispacted beam «q>erlmenta used "to assess the 
accuracy" of the present prediction method: the fragment is not a 2-D frag- 
ment and the observed structural response behavior is not strictly 2-D in 
character; thus, a rigorous theoretical-ej^erimental con^>arison cannot en- 
sue. However, with recent improvements [321 in the experimental techniques 
and apparatus for conducting small-scale ^act experiments having the 
scope and objectives of those of Ref. 16, ;it is now feasible to conduct 
sisixlar beam intact experiments with a solid cylindrical 2-D fragment rather 
than the one-inch diameter steel sphere en^loyed in the Ref. 16 experiments, 
and thus to obtain more nearly 2-D structural response data for ia5)acted 
simple structures. ' 

Clearly, the present 2-D analysis is inherently incapable of pre- 
dicting the 3-D deformations e;q>ected and observed to have occurred in a 
region centered about the "inq>act point" in steel-sphere-impact speci- 
ments CB-13, CB-18, and CB-16. For the 2-D predictions, the Idealized 
fragment was defined to be a solid 2-D cylindrical fragment of one-inch 
diameter and to have the same mass and isqjact velocity as the steel sphere 
in each case. Note that the present 2-D 43-element Q-SH predictions for 
the peak strain at the midspan intact station (x « 0) for cases CB-18 and 
CB-16 is about 16 per cent and 16.4 p^ cent, respectively; these 2-D 
predicted strains exceed the range of validity of the present analysis 
formulation and computer code implementation and may be ve^ different from 
actual peak extensional stredjts which are present in the clearly 3-D 
deformation region there. Experimentally, both beam CB-18 and CB-16 
are close to the rupture threshold for these 3-D deformation conditions 
of strain and stress; a reliable quantitative characterization of the 
biaxxal or triaxial strain state to define this rupture threshold for this 
6061-T651 aluminum material is lacking. Severe necking of the material 
near the fracture station of specimen CB-16 is evident. It is uncertain 
also at which extensional strain level one would observe threshold rupture 
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1: 5.2 Modeling and A nalysis Consideration.. 
it be assumed that: 

(a), a specific 2-D type structure is to be subjected to impact^ 

by a 2^ fragment as accommodated in the CIVM-JET 4B com- 
puter program. 


<b) the uniaxial static stress-strain curve of the struc- 
tural material has been specified and fitted by piece- 
wise linear segments for use in the mechanical sublayer 
model (also that values are specified for the strain 
rate parameters D and p), and 

(c) the sire, mass, and kinetic energy of each attacking 
fragment are given. 

The next matter to be considered is the selection of the number and 
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size of the finite elements to be used to model the structure. The elements 
must be numerous enough and of smaU enough size: (a) to provide a con- 
verged solution at least in an engineering sense and (b) to acconmodate with 
reasonable accuracy the "physically plausible" ia^wict-affected effective 
length region ■ Impacted structure to which mrn—n^tm ia 

transferred during a single fragment impact. The latter means the 
length of each finite element near an io^ct station should not — 
about 2h (i.e., twice the thickness); otherwise, the j 

velocity increment to iiig>act-affected zone of the structure will become j 
in?»lauslbly small and the resulting predicted response will be less severe 
than realistic in this region. 

Further, if one has a structural configuration or arrangement 
wherein it is apparent that Important levels of strain will be present 
and will be changing rapidly with spanwise location, relatively many short- 
length elements should be used in these regions to model the behavior there 
properly, whereas fewer larger-size elements could be used in regions of 
spatially slowly changing strains. Bowever, as pointed out earlier (Sub- 
2.2; also see Section 3), if one is considering a containment ring 
which may be ispacted by fragments at many locations (which may also change 
-as time progresses) around the periphery, the clearly-logical choice is to 
use only equal-l^gth elements. 

If one uses CIVM-JET 4B which enqiloys the timewise central differ- 
ence operator, the solution time step size At must be chosen to be 

~ 0-8(2/g)„^^) where is the maximum natural frequency of the 

mathematical model of the structure for linear h^avior. One finds that 

for a given total number of 4DOF/node finite elements, the smallest u 

. , . max 

occurs tdien equal-length elements are employed; this, in turn, permits the 

allowable At to be used for avoiding cooq>utational blow— up. The 

same number of nonunifonn elements will require one to use a smaller At 

(and hence to incur more computing expense) in order to avoid conputa- 

tional disaster. in most cases, therefore, the use of equal-lenoth 
V 

Alternatively, one can enploy a laurger At while avoiding blow-«p of conpu- 
tational roundoff error by using a different type of timewise finite-differ- 
ence operator (see Subsection 2.5.3). 
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-elements to model the structure will be advisable. 


As described in Refs. 4 and 17, the approximate impact-inter- 
-action-model eo^loyed makes use of ln^mlse-momentum and energy consider- 
ations to predict the momentum transferred to the impact-affected portion 
of the structure from each impact by a fragment. The structure, in turn, 
has been represented by finite elements whose mass matrix properties are 
evaluated from kinetic energy expressions using an assumed velocity dis- 
tribution in each element conslstrat with this basic assumed displacement 
field of the element; accordingly, the resulting mass matrix is called the 
"variationally consistent mass matrix" - this is a symmetric non-diagonal 
matrix. Reference 17 e;^lored the use of this non-diagonal mass matrix 
for these "momentum transfer" calculations and encountered erratic be- 
havior (see Ref. 17 for further details); hence, a diagonalized mass was 
constructed I17J and led to plausible- behavior. Thus, the CIVM-JET 4B code 
[4J has utilized a diagonalized mass matrix model for the momentum trans- 
fer calculations. Finally, the momentum transfer prediction can be 
carried out by assuming- the local ia 5 >act to be perfectly elastic, per- 
fectly inelastic, or intermediate between these conditions by assigning 
as input a value for the coefficient of restitution e of, respectively, 
e - 1, e - 0, or 0 < e < 1; the Ref. 17 studies show results for various 
e values and conclude that the perfectly-elastlc case e - 1 is the most 
reasonable engineering choice for most of the types of low-sp^d intact 
conditions under consideration h^e. The user of CIVM-JET 4B, however, 
can specify any value of e (0 £ e £ 1) as he wishes. 


— Ose of CIVM-JCT 4B versus CIVMtJET 5B 


Both Of these computer programs permit one to predict the large- 
deflection, elastic-plastic transient responses of 2-0 structures which are 
subjected to low-speed impact attack by one or more idealized fragments 
CIVM-JET 4B is restricted to the analysis of single layer (single material) 
structures and employs the timewise central difference operator for the 
transient solution calculation, whereas CIVM-JET SB can deal with multilayer 
nniltimaterial 2-D structures and uses a different timewise finite-difference 
operator (the HOubolt operator) for the timewise solution. This latter 
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operator is unconditionally stable and thus permits one to use a very large 
time-step sire, if desired, without encountering coa^utational roundoff 
-error blow-up r this ’^advantage” was the principal motivation for the 
adoption of the Houbolt operator in CIVM-JET 5B. However, the occurrence 
of is^ct and the attendant accounting for the in^act-ioparted velocities 
to the ** locally affected" portion of the structures results in effectively 
changing the timewise o^rator at and for a short time following each im- 
pact from the BOubolt operator to an operator idiich is no longer uncondition 
ally stable [10]. Therefore, the allowable At is also restricted althou^ 
not as severely as in CIVM-JET 4B. 

Some elaboration: on this "correction for impact" matter follows. The 
CIVM-JET 4B program calculates the exact time of contact between the frag- 
ment and a local region of the impacted structure. This time will always 
occur at some sub-time increment which is smaller than the time-step size 
en^loyed for the solution procedure. The use of the centr 2 Ll difference 
operator allows an accurate correction to be applied to the structiire during 
the remaining portion of the time-step increment. These corrections are 
done in a manner which allows a correct distribution of energy and momentum 
in the system even though a "nonuniform time step" has been employed. The 
Houbolt operator will not advance the solution properly if a nonuniform time 
step correction is isposed; hence, a different ispact interaction procedure 
had to be introduced. The CIVM-JET SB program performs an "approximate 
time of contact correction" to the structure displacement field; an im- 
pact is detected to have occxirred within, a time step by noting the presence 
of an overlapping of the fragment and the structure. The corrections are 
then applied in CIVH-JET SB as though the initial contact occurred at the 
beginning of the time step. A response prediction iteration is then per- 
formed at this time to gxiarantee a cocplete correction, and then the struc- 
tural response solution continues through time until another fragment- 
structure overlap is encountered. Since the applied iapact-interaction 
model in CIVM-JET SB is an approximation of the model used in CIVM-JET 4B, 
a oonparison between the two solution techniques shotzld prove useful for 
future inpact- interaction models. Tt should be noted that the "exact time 
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of contact" solution gives a aore "realistic" analysis of the la^act-inter- 
actlon, but it becones very tla>e-consuaiag as the nuaber of. elesients used 
to n od el the is^acted structure increases. Incidentally* the application 
of this "exact tine of contact" is^ct>interaction solution technique 
tdiai using plate elements (10, for ex^le] for analyzing plate or panel 
response to fragment impact would be prohibitive. 

factor that influwces..the largest allowable At in CIVM> 

JET SB is the fact that the BOubolt operator is of implicit type; there- 
fore* the internal forces associated with large deflections and plastic 
behavior depend upon the solution being sought. These terms* hence* 
have been approximated by extrapolating linearly from known solutions at 
the immediate two earlier time Instaijth without iterating to convergence 
within each time step. Therefore, these pseudo-forces may become pro- 
gressively less and less accurate particularly if leurge At time steps 
are used. . 

Since CIVM-JET SB requires much more storage and computing for each 
step of calculation than does CIVM-JET 4B* to be cost competitive 
with — 4B, one must use for —SB a At at least 4 t im es larger than the 
allowable At for CIVM-JET 4B. Note that CIVM-JET 4B uses the unconven- 
tional form of the equations of motion while CIVM-JET SB enploys the 
equations of motion in their conventional form which involves more conputing 
per step; both of these formulations [4* S* 17J are given concisely in 

j)ppendix~A for the reader’s convenience. 

Figure 21 compares for the C»-18 experiment predictions of the CIVM- 
JET 4B and CIVM-JET SB computer programs using the EL-SH Fit A material 
model and 43 equal-length finite elements. The predictions for the CIVM- 
JET SB model are given for two time step sizes: At « O.S ysec and At- 2.S 
usee. The first time step size is identical to that used in the CIVM-JET 4B 
calculation which uses the central difference temporal operator. The second 
fima, step size represents the point at which the CIVM-JET SB model becomes 
cost effective compared with the CIVM-JET 4B model. 



Predictions are plotted for four spanwise locations (x « 0.0, 0.6, 1.5, 
and 3.7 in) which include the irpact region, the mid-region ^ich is 
dominated by 2-0 deformation, and the region close to the cleuiped boundary. 
In each region the solutions are qualitatively similar to each other; for 
engineering purposes these are also quantitatively the same. 

A CIVH-JET SB calculation was carried out by using a At size which 
was about 10 times that permitted by the central' difference method used in 
CrWf-OTET 4B; this solution exhibited instability when iB?>act and correc- 
tion was encountered for the CB-18 exaiiq>le. It is; not known for the CB-18 
conditions the largest At permissible for well-behaved calculations when 
ClVMrJET SB is used; 4 times and 10 times larger At than permissible for 
CIVM-JET 4B showed, respectively, well-behaved and unstable behavior for 
CIVH-JET SB calculations. 

The Houbolt operator as en?>loyed with the CIVM-JET SB iapact 
modeling is seen to be an alternative solution technique to CIVM-JET 4B, 
but care must be exercised that a cost effective time step be determined 
without encountering a region of instability. 

The dVM-JET SB program «;as not used in further intact correlations 
and studies performed in this report. The SB program con?>ares favorably 
with the 4B program for predictive accuracy; homver, the CIVM-JET 4B pro- 
gram had been exercised for a longer period of time and has provided - " 
plausible results for several intact-interaction wjdelings. These con- 
siderations led to the decision to study jet engine ij^act problems and 
the int>act-interaction parameter variations (see Section 4) with the 
CIVM-^JET 4B program. 

The preceding discussion has dealt with the principal considera- 
tions for deciding whether to use ClVM-JET 4B or CIVM-JET SB to predict 
int>act-induced transient structural responses of 2-D beam and/or ring 
structures. Either code can be used to analyze single-layer structures 
but only CIVM-JET SB can be es^loyed to analyze multilayer, 
structures. Earlier it was pointed out (see also Subsection A.2 of Ref .4) 
that a diagonalized mciss matrix modeling of the structure is needed for 
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the prediction of ncnentum transfer from the fragment to only the trans- 
lational degrees of freedom of the intact-affected region of the struc- 
ture. Now, some discussion is given of various available ctoices for this 
diagonalized, mass matrix modeling, and their consequences. 

First, recall that the present beam and/or ring finite elements 
have 4 DOF/node; these are the in-plane and normal^direction transla- 
tional displacements v and w, respectively, and the "extension" X 
"rotr^-ion" i|( (see Appendix A and Fig. A.l). In deriving the mass matrix 
for a single finite element from the kinetic energy expression or frpm 
the Principle of Virtual Work, one obtains a non-diagonal mass matri^ 
termed the "variationally consistent mass matrix .one obtai n s mass matrix 
entries on the diagonal associated with each of the 8 DOF's (4 at each 
end), as well as the off-diagonal "mass coupling" terms. The use of ^e 
consistent mass matrix as well as a mass matrix diagon al ized in various 
ways has been studied (17, 33-3C]; both types of matrices are used 
widely with reasonably similar transient response results. Where feasible 
or where considerations such as ia^act-induced m o men tum transfer calcula- 
tions are needed, the conq>utationally more efficient diagonalized mass 
matrix is preferred. Further, it was found that the predicted transient 
responses of 2-D structures to transiently applied loads »*en a diagonal 
mass matrix was used were affected only slightly by different choices 
for -the diagonal entries associated with the non-transl ation al_degrees 
of freedom of the element as long as the diagonal entries associated 
with the translational degrees of freedom of the element were kept the 
same [17, 33-36]. 

For the case of single-layer structures of uniform thickness, one 
can readily define the diagonal mass matrix entries associated with the 
translational degrees of freedom v and w at each nodes this is sin?)ly 
one-half of the mass of the element. The choice of the proper diagonal 
mass matrix entries corresponding to the non-translational and y de- 
grees of freedom at each node is less apparent; these entries, however. 



are generated automatically in the consistent nass matrix calculation. For 
a beaa element which consists of several variable-thickness layers of dif- 
-ferent-material, the "intuitive selection" of even the "translational en- 
tries" in a fully diagonalized mass matrix becomes a much more diffi- 
cult and complex undertaking; however, the cailcxilation of all of the mass 
matriic entries in the consistent mass matrix remains very simple and automatic. 
Therefore, it is convenient and natural to ask how one might select a diagon- 
alized mass matrix by making use of the automatically-generated consistent mass 
matrix. This matter has been e;^>lored by Key and Beissinger [33], Clough i35], 
and Hinton [36]; there is general agreement on the selection of those entries 
associated with the translational degrees of freedom, but various selections 
have been suggested for the diagonal entries associated with the remaining 
degrees of freedom. 

The epproech suggested by Key end Bexssijiger [33] end by Hinton 
[36] mey be summarized as: the procedure of lunping (mass matrix 
diagonalization) recommended in view of the infinite possibilities is to 
con?)ute the diagonal terms of the consistent mass matrix and then to scale 
these terms so as to preserve the total mass of the element; that is# so 
that the sum of the mass entries at both ends of the element for each 
translational degrees of freedom shall equal the total mass [36]. A 

orocedure has been proposed by Key and Beissinger [33]: . the trans* 
lational entries for the diagonalized mass matrix are determined from-the 
consistent mass matrix by adding the corresponding off-diagonal transla- 
tional term to the diagonal term ~ this does result in identifying the 
total beam mass as the sum of entries- at each end of the— element for a 
given translational DOF. However, no clearly superior means for choosing 
the non-translational entries for the diagonalized mass matrix from the 
consistent mass matrix is apparent. One possibility suggested by Kay and ^ 
Beissinger is that these non— translational diagonal entries be scaled 
from the consistent mass matrix so that the highest natural frequency is 
matched with the highest natural frequency of the structure modeled with 
consistent mass (and stiffness) matrices. Since no clearly superior 
transient response results are demonstrated for such’ a choice and since 
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the preservation of the same largest natural frequency uould impose the 
same At limit if one were to conq.ute the transient structural response 
by using- thei-timewise central-difference operator, this is regarded as an 
unattractive choice. Instead, a plausible choice for these non-transla- 
tional entries of the diagonalired mass matrix is desired such that the 
highest natural frequency of th4 modeled structure becomes "as low as 
possible" — this would have th^ effect of increasing the permissible 
time step size At for central-difference operator calculations, with no 
detrimental effects when one employs an implicit timewise operator such 
as the HOubolt operator. 

Three candidate mass-matrix diagonalizing schemes are considered 
here. These are described concisely in the following. For clarity, 
consider a single initially-straight beam element with the following 
4 degrees of freedom (q^^) at each end: 


End 2 (or Node 2) 


< 1 . 


V 

w 


T?nd 1 (or Node 1) 


« V 

, . dw 

= an 

^ = ln ^8 “ ^ 

Mass Matrix Diaqonal ization Method 1 -- - - 

This is an "intuitive engineering method" whi^h can be used readily 

only for single-layer elements of either constant or linearly-varying _ 

thickness. The same translational masses are assigned to DOF's and q^ 

and another value to DOF's q^ and q^; the sum m^^^^ + m^^ = + m^g - 

total mass of the beam. Non-translational DOF's q 3 and q^ are assigned 

n and m as being the same "estimated effective rotary inertia" for 
33 44 


the beam; a similar assignment is made for m^^ and m^g. See Appendix A 
of Ref. 4 For a fuller discussion. Concisely, these diagon al mass matrix 
entries are (non-diagoneJ. entries are all zero) : 
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where the thickness-dependent constants and are given by 
(2h2 + h^) 

S ” 3 (hj^ + h^) 

hj + 4 h h, + h, 

• Q m ■■■ ■ ■ 

"^2 36(h^ ^2^ 


In these expressions p is the mass per unit volume of the beam element 
b is the width of the beam, I is the length of the beam element and h^ 
and h^ are the thickness values for the beam at ends 1 and 2, respec- 
tively. For present purposes, this description is considered to be 
adequate; see Ref. 4 for further details. 


Mass Matrix Diagonatlization Method 2 


This method defines a diagonal mass matrix by using the 
entries already computed from the consistent mass matrix m^^ , as 
follows : 
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Hence# this method is readily applied to multilayer# multimaterial vari-^ 
able thickness structures. 

Mass Matrix Diagonaliration Method 3 

This toethod is similar to Method 2# being different only with 
respect to how the non**translational degree -o£r freedom entries are 
calciilated. All diagonalized mass matrix entries are calculated by the 
rule 

D C ^ C 
• n.i + 

Renee# 

D C ^ C 

D C ^ C 

®22 “ “22 ®26 

0 C ^ C 

“33 “ ®33 “37 

D C ^ C 

44 44 48 

These three mass matrix diagonalizations as well as the consistent 
mass matrix %rere used to calculate the maximum natural frequency of 
clamped-ended beam CB-4 when modeled with twenty 4 DOF/node elements per 
half-span ~ iaqposing syinnetry conditions at midspan. The results are 
as follows: - - 


Mass Matrix 

<u 

max 

(rad/sec) 

(At) - — ^ 
cr-^ 0) 

(Usec)®“ 

Consistent 

.7678x10^ 

0.26 

Diag. Method 1 
(Intuitive) 

.2165x10^ 

0.92 

Diag. Method 2 

.42IOxlo’ 

0.47 

Diag. Method 3 

.9874x10^ 

0.20 


Since only diagonalized mass matrices are employed in the calculations of 
structural response to fragment impact# it is seen that Method 1 permits 
the largest allowable At to be used if one were to eoploy CIVM-JET 4B 
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which uses the central difference timewise operator; however. Method I has 
been developed only for single*layer elements. Hence, for multimaterial, 
multilayer variable-*thickness elements. Methods 2 or 3 could be used; of 
these . Method 2 appears to be the. preferred choice since it will permit the 
use of a larger At than will Method 3. 

Some further cos^^arisons are given in Table 1 for some predictions 
of large elastic-plastic transient response of inpulsively-loaded beam CB-4 
The half span was modeled with 20 equal-length 4D0F/node elements, the 
EL-SH Fit A stress-str 2 dji behavior was used, and a uniform lateral velocity 
of 10,707 in/sec was applied to the center node and to the next four 
nodes. Sho%m in Table 1 are predicted quantities (maximum Gatissian sta- 
tion strain, maximum nodal strain, elastic energy, kinetic energy, plastic 
work, and midspan w-displacement) at 100 ysec after the (sudden) initial 
velocity was applied ~ for the use of: (a) the consistent, (b) the 
Method 2 diagonalized, and (c) the Method 3 diagonalized mass matrix. 

It is seen that the response quantities predicted at t « 100 ysec are 
close in all three cases; the two predictions which utilize the diagon- 
alized mass matrix are very close to each other. Hence, it is recom- 
mended that until better mass diagonalization schemes are developed. 

Methods 1 and 2 should be adopted where applicable. These are the 
omthods employed, respective^, in CIVM-JET 4B and CIVM-JET SB. 



SECTION 3 


CONTAlNHENT-aiNG EESPONSE TO T58 TCRBINB 
ROTOR TRI-HOB BURST ATTACK 

3 . 1 Problam Definition 

In th« spln-chaobnr facility of th« Naval Air Propulsion Test Center 
(MAPTC) numerous experiments have been conducted wherein various types of 
actual and s^lated engine rotors have been rotated at high rpm and caused 
to fail iji various ways. High-speed Dynafax photographs nominally at 
35,000 pictures per second have recorded the motions of the resulting frag- 
ments and: their impact and interaction with a variety of containment-ring 
structurea [12-14 and 37. for exiunplel. Some of these containm^it rings 
were of single-layer construction} others were of multilayer multimaterial 
construction. In some eases, strain gages were applied to the containment 
ring, and the resulting transient strains and permanent strains were recorded 

Selected for discussion and analysis in this section of this study 
is NAPTC Test 201 [13.141 in which a 4130 east steel* containment ring of 
0.625-in thickness. 1.50-in axial length, and 15.00-in inside diameter and 
weighing 12. 83. pounds rested horizontally on smooth support wires and en- 
circled a T58 turbine rotor which was caused to fail in 3 equal 120-degree 
segments at about 19.859 rpm and to impact against this steel containment 
ring. Circumferentially oriented strain gages were attached to the ring 
at various circumferential locations at the ring-s midwidth position. 

Transient strains were tape recorded for 10 such gages; transient strains 
from 4 of these 10 gages were also photographically recorded from swept 
oscill . rxjpe displays. In addition after the test, permanent strains indi- 
cated by the four gages which survived intact were recorded. 

Given in Table 2 are the weight and geometric data defining the con- 
tainment ring, the rotor-burst fragment properties, and the test conditions 
for NAPTC Test 201. Note that the I.D. of t.he ring was 15.00 inches while 

* 

From NAPTC* s ACIPCO 2 billet [141. 
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the 0.0. of the TS8 turbine rotor was 14.00 inches; this 0.5-in clearance 
between the ring and the blade tip is untypically large but was selected to 
permit better. photographic_clarity-and-de£inition. Each fragment consisted 
of a "120-degree" sector of the rim with 17 attached blades; the distance 
from the axis of rotation of the rotor to the CG of the fragment was 
2.797 inches. At the rotor burst tpm of 19,859,. the translational velocity 
at the CG of each fragment was 5816.7 In/sec «fhlle the blade-tip speed was 
14,557.2 in/sec. The resulting total kinetic energy of the 3 released frag- 
ments at burst was 908,820 in-lb, of which 476,766 in-lb was translational 
and 432,054 in-lb was rotational. Bence, each fragment at burst had 
nominally 158,922 in-lb of translational and 144,018 in-lb of rotational 
kinetic energy. 

Shown in Pig. 22a is a schematic of the TS8 power turbine rotor modi- 
fied to undergo a tri-hub burst, and Pig i 22b shows a sketch of a typical 
fragment before in^jact. A photograph of one of the 3 fragments after the 
test is shown in Pig. 22c. The quality of the Dynafax photographs (taken 
at 35,000 pictures per second) of the iapact of the 3 fragments and their 
interaction with the ring happened to turn out, unfortunately, to leave 
much to be desired; however, shown in Pig. 23 are a few sample photographs: 
(a) pre-test, (b) after rotor burst but before in?>act of the fragments with 
the ring, and at about 200, 514, and 1286 microseconds j^fter initial impact 
in parts (c) , (d) , and (e) , respectively. Pigure 24 shows the permanently 
deformed containment ring; note that position marks have been provided at 
5-degree intervals along the ring's midthickness location. Also shown in 
Pig. 24 is the post-test condition of the 3 attacking fragments. Note that 
only one of the fragments still possesses many of its original blades. The 
Dynafax pictures obtained show that this denuding occurred late, after the 
fragments fell by gravity below the plane of the containment ring and sub- 
sequently struck the heavy steel sides of the test chamber, resulting in 
greatly enhanced "fragment damage". Of course, very severe deformation and 
curling of many blades of each fragment can be seen to have occurred trtiile 
the containment ring is still suffering the in 5 )act attack. This attack con- 
tinues until well after the peak response of the containment ring has been 
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reached) only later do the fragments "escape the ring" by falling- below 
"the containment plane". 

Figure 25 indicates the circumferential locations of stredn gages 
which were-used-to-measure~transient and permanent strains on the outer sur- 
face of the ring. These locations are defined by stations 1 through 72 
marked at 5-degree interrals around the ring; station 1 is assigned 6 ■ 0°, 
station 5 represents 0 - 20°, etc. Note ^t permanent strains, as follows, 
were obtained from intact strain gages eifter the test at only four locations 

Station 6 (deg) ^ Permanent Strain (%) 


9 

40 


- 0.95 

13 

60 


- 3.23 

33 

160 

3 - 

- 1.07 

37 

180 


- 3.31 


Static uniaxial tensile stress-strain meuurements for coupons of 
4130 cast steel from NAPTC’s National Forge billet were supplied [141; 
these data were analyzed and used for the theoretical analysis of NAPTC 
Test 201 since according to Ref. 14 the Test 201 ring material is almost 
identic 2 d to that of the National Forge billet. Accordingly, these stress-, 
strain data were approximated by piecewise linear segments defined by the 
following stress-strain (o,e) pairs: a,e ■ (0 psi, 0 in/in) ; (80,950 psi, 
.00279).; (105.,.300- psi, .0225); and (121, POO psi, .200). This c,€ fit is 
used~in the mechanical sublayer material-behavior model ud is referred to - - 
as elastic, strain-hardening (EL-SH) . When strain- rate "effects are taken 
into account, the perfectly-plastic yield stress of each (kth) mechanical 
sublayer for strain rate e is taken to be that given by Eq. 2.1. Since 
measurements have not been made of the strain-rate dependence of the a,t 
behavior of the steel used in the NAPTC Test 201 ring, it has been assumed 
that .its strain-rate dependence is approximately that of mild steel [25] 
for which. D = 40.4 sec ^ and p = 5. ' 
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outer ria of the disk part of the fragment* At this stage ^e "effective 
radius” of the attacking fragment is approaching the cited "minimum plausible 
value” of r^ - 2*555 in* These choices are depicted schematically in Fig* 27* 

The use of an implausiblylhrge idealized rigid circular fragment 
• would clearly constrain the containment ring and unrealistically limit the 
amount of its deformation in the impact/contact region of each fragment* 

On the other hand, the physical dimensions of the rigid disk portion of the 
fragment and the region occupied by severely curled over and mutually con- 
tacting blades make the ' choice of an r^ smaller than choice (a) to be clearly 
illogical* Hence, for :a severe threshold-containment type of fragment attack, 
one expects that the u4e of an idealized fragment of fixed r^. ^ 2*555 in 
would result in a quite' reasonable prediction of the maximum response of the 
ring; however, it is in^ssible physically for this type of id e alized model 
to give a reasonable approximation of the actual transient response history 
of the containment ring — actual is^ct %rill occur much e 2 u:lier and, 
initially, much more gently than for the idealized fragment* 

Figure 28 depicts the geometric test and modeling data for the 4130 
steel containment ring subjected to tri-hiib T58 turbine rotor burst in NAPTC 
Test 201* The ring is represented by 48 elements and EL-SH-SR material be- 
havior is taken into account* Three equal-size idealized fragments are used; 
all three fragments are assumed to inpact initiaJ.ly simultaneously at equi- 
distant circxjmferential locations* 

The effects of idealized- fragment size can be illustrated most conveni- 
ently perhaps by conparing the extreme deformed- ring configuration for each of 
these 3 cases, as shorn in Fig* 29a; the time at which t^t extreme deformed 
state occurs is different for each of these 3 cases. For these" 3 extreme de- 
fored states. Fig. 29b compares the circumferential distribution of the cir- 
•cumferential -extension strain on the outer starface of the ring. Clearly, 

the extr^se response of the ring is greatest for the smallest and is lea^ for 
the largest idealized fragment size. In these three calculations, it was 
assumed that there was zero friction between each fragment and the ring* 

It may be of interest to examine the nataire of the transient strains 
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which are predicted to occur at various circumferential locations of the 
containment ring. As a convenient means of identifying circumferential 
locations. Fig. 30 shows a predicted deformed ring configuration (and the 
3 idealized fragments) at 1000 microseconds after initial impact; shown ailso 
are the element identification manbers. Predicted inner-surface and outer- 
surface circumferential transient strains are shown in Fig. 31 at elements 1, 
4, 6, 9, 11, and 47 for the case r^ * 2.555 in and y - 0. One may interpret 
the mean of the inner siarface and outer-surface strains as representing 
roughly the "membrane strain** portion, and the deviation from this mean as 
the part arising from bending. 

The effect of using a friction coefficient U value of 0.3* is illus- 
trated in Fig. 32a where the deform^ ring configurations for U » 0.3 and 
y * 0 at 1200 microseconds after initial in?>act are shown for r^ • 2.555 in. 
Since these deformed ring profiles are shown with respect to the fixed y, z 
coordinates, it is seen that they are "slightly rotated" with respect to each 
other but the deformation severity is comparable. Shown in Fig. 32b at that 
same time instant is the outer-sxirface distribution of predicted circumferen- 
tial strain; here also, it is seen that the effects of friction between the 
idealized fragment and the containment ring have little effect upon the magni- 
tude of the peak strains predicted at this instant. Hence, most of the calcu- 
lation cases -carried out assimied . frictionless impact. 

Summarized in Tablo-4 are the pertinent data which characterize each 
of the various CIVM-JET 4B calculations made to predict the response of the 
present containment ring. The time increment size used in every case was 
2.5 microseconds. - ^ - 

3.3 Comparisons of Predictions with Experiment . 

For reasons noted in Subsection 3.2, predicted time histories of im- 
pact-induced containment ring strains cannot be compared rationally with 
measured transient strains for this ring- fragment inpact problem. Also, 

*This choice as a plausible and reasonable value is indicated by the studies 

reported in Ref* 15. 
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between the fragment and the impacted structure has very little effect upon 
the predicted transient response of a containment ring. 

Of the various modeling parameters and values which the analyst must 
select, that having the greatest effect on the peak predicted structural re- 
sponse of the containment ring is the idealized fragment size itself. 

P'inally, two reminders concerning the limitations of the present CIVM- 
JET 4B analysis arid code are pertinent. First, the analysis applies strictly 
only to the two-d^nsional type of structural response; hence, cases in- 
volving an import^t degree of 3-D structural response cannot be, represented 
by the present analysis and code. Second, the governing equations for this 
analysis apply to ^large-deflection and large-rotation elastic-plastic transient 
structural response but the strains themselves must be small. An upper limit 
on the strains for which this analyses is valid cannot be specified precisely; 
however, predicted strain values exceeding roughly 6 to 10 per cent should be 
regarded with caution. Clearly, rupture threshold predictions for ductile 
metal C/D structures will involve substantially larger strains; accordingly, 
an extension to accommodate large strain behavior properly will be necessary 
to achieve that end rationally cuid reliably. 



SECTION 4 


COJTAINMENT STR0CTOKB PAWVMETRXC EFFECTS ON 
STRDCTORAL RESP<»SB INDOCEO BY FRAGMENT IMPACT 

Sections 2 and 3 of this report have been devoted to discussing 
numerical methods for predicting;, 2-d transient structural responses of, 
respectively, (1) clamped-ended beams subjected to rigid-fragment (steel 
sphere) impact and (2) containment rings subjected to impact attack by 
ideaUsed engine rotor fragment^. The capabiUties and limitations of the 
associated computer codes CIVM-^ 4B and CIVM-JET SB were reviewed together 
with guidelines for selecting an;, appropriate structural modal so as to obtain 
reliable engineering predictions of fragment-impact-induced structural 
response (within the range of applicabiUty of the prediction method as 
actually implemented in these computer codes). Also, predictions were 
ccapared with representative experimental data. 

Since the scope and complexity of actual fragment attack against 
containment structure usually greatly exceeds that which is feasible to 
include in mathematical/numerical prediction models, it is often useful to 
employ a judicious combination of (1) mathematical/numerical predictions and 
(2) actual experiments — in order to study in a cost-effective manner the 
effects of each of the numerous variables upon, for example, the (least) 
well, SI? and/or cost of a containment structure which is to be selected to 
achieve threshold containment"^ of a given fragment at^ckj similar results 
for a series of different postulated fragment attacks may also be of interest. 
The use of mathematical/numerical predictions alone is hasardous because of- 
the possible inadvertent omission of important ingredients present in the 
actual physical problem. On the other hand, the use of experiment alone is 
inadvisable because of the time and large cost per test as well as the great 
number of tests required to explore the effects of each of the potentially- 
important parameters and variables of the problan. 


That is, the dividing line between fragment containment and non-containment. 
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Zn this section the iise of dimensional analysis in studying the effects 
of various pertinent parameters upon the containment threshold response of 
each of various 2-*d containmwt structures to fragment impact attack is 
explored with reference to; (1) mathematical/numerical predictions and 
(2) actual experiments* The dimensional analysis is discussed first in 
terms of a rather general situation involving fragment attack against a 
broad class of containment structures* Then the analysis is specialized and 
applied to (a) NAFTC experimental studies (12^13,371 of 2-d containment 
ring responses to specific types of engine- rotor- fragment impact attack and 
(b) MIT-ASRL experimental (16] and numerical studies of steel-sphere impact 
attack against a simple ductile metal beam* Finally, the effects of scale 
are discussed; that is, the question of determining the physical parameter 
values for a "similar contetinment ring" required to achieve threshold contain- 
ment when the "similar fragment attack" emanates from a rotor of N times the 
size of a given reference rotor of the same material is examined. 


4. 1 Dimensional Analysis Considerations 


There are two basic mutually- complementary methods that can be used to 
quantify the containment threshold capability of a given containment structure 
which is subjected to a given fragment impact attack condition, and to 
identify the major parameters affecting a -particular-structure’s containment 
potential: 

(1) A mathematical model representing the governing eq uations of the 
large-deflection elastic-plastic responses of the containment 

» structure and the ' impacting fragments, and of the impact- interaction 

behavior can be devised. These equations can then be rendered 
dimens ionless^ .These dimensionless parameters can be varied in 
'( magnitude and th^ subsequent solutions of the mathematical model 

will determine the effects that these parameters have on the sought 
threshold-containment conditions. 

(2) A physical approach can be used in which all of the "physical 
parameters" thought to be involved in the impact- interaction and 
response are identified and formulated into a set of dimensionless 
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preliminary containment studies [39K a set of dimensionless parameters can 
be constructed from a set of pertinent variables by application of the 
Buckingham^Pi theorem. If n is the number of variables present and r is the 
number of least-order dimensions represented by these variables (i.e., mass, 
length, time, etc.), a complete set of m dimensionless parameters can be 
constructed where m « n - r. This set of dimensionless parameters will by 
no means be unique; but, by physically examining and justifying each dimen- 
sionless parameter, a meaningful set of m parameters can be constructed. The 
effect of each of these parameters, can be studied individually by conducting = 
a series of (1) experiments or (2) mathematical solutions, holding all other i 
parameters constant. Of course, one could and should conduct partial 
factorial experiments to examine the possible presence of synergestic effects; 
This procedure is permissible if one knows that there are no synergestic 
effects between the parameters. In the absence of such knowledge, some type 
of factorial experiment would be appropriate. 

4. 1. 1 General Description 

It is necessary first to identify all of the principal variables 
influencing a general impact- interaction and response of a cont 2 uner or 
deflector structure which is subjected to fragment impact attack. Before 
these variables are assembled into a useful set of dimensionless parameters, 
the variable list should be reduced by eliminating those variables whose 
effect on the overall interaction is physically reasoned (or can be shown) 
to be insignificant. Once a ** complete** set-of_priiziary„variable^ is identi- 
fied, they can be combined based on physical relationships into~i convenient 
set of dimensionless parameters. 

The variables describing an ** impactrintej^ction** can be divided into two 
distinct groupsV (11 thos'e^ variables unique to the ’impacted structure, and 
(2) those vari2d3les unique to the impacting fragxnent'CS) . These variables 
can be further subdivided into geometric, )cinematic/environmental, and 
material variables. 

The containment structure is defined geometrically by its length, width, 
^d thickness; if the structure is conceptualized as a complete (containment) 
ring, the radius of the ring becomes a primary variable rather than the - 
length. If the container is a section of a doubly-cxirved 3-D body, two radii 
of curvature will become primary geometric quantities. For purposes of 
simplicity, the containment structures treated here are regaunled as being 
three-dimensional thin-shell bodies of revolution (i.e., rings, cylinders, 
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or spheroidal sections, see Pig. 35), or as plates or beans. A thickness 
measure can be defined for any of these structures (here all are assumed to 

bo of constant thickness); appropriate lengths, widths, or radii also can be 
defined. 

For Illustrative purposes the attacking fragments are considered to be 
portions of a jet engine turbine rotor,. Pig. 36. The fragment could be 
either a single blade, a segment of the hub with attached blades, or a 
bladed-disk sector. All three of these potential fragments can be described 
geometrically by a width (i.e,, axial projection which is assumed to be 
constant) , the distance from the rotor axis of rotation to the fragment CG, 
and the length of the fragment from its OC to the blade tip, see Fig. 36. 

The movement of the fragment relative to the container must be defined 
geometrically, since the container is initially motionless in the reference 
space used for the analysis, the fragment will possess all of the initial 
translational and angular velocity of the system. The angular velocity of 
the fragment is a primary variable while the translational velocity of the 
fragment is defined by the angular velocity and the location of the fragment 
center of gravity, r„. 

CG 

geometric a nd kinematic variables and their dimensions are 
summarized below for both the fragment and the containment structure (or 
container) : 


Geometric and Kinenatic Variables 

Quantify 

Descripti(5n 

Units 

u 

0 

r (or 1' ) • 
e c 

Major Radius (or Spanj - - 

IB 

c 

•H 

r (or w ) 

2c c 

Minor Radius (or Width) 

■ 

o 

o 

h 

c 

: ThicJcness 

■ 



Length from Fragment CG to Blade Tip 

B 

■ U 

^CG 

Distance from Axis of Rotation to 

L 

e 

o 

g 


Fragment CG 

L 





U 

Cu 

''f 

Width (Axial Projection) 

L 



Angular Velocity 

T*^ 


n 

Numter of Fragments ' 


1 TOTAL: 8 Variables 
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The material propertiea of both the fragment and the container are not 
as readily quantified as are the gecnetric and kinenatic quantities. Both 


the fragment and the container have their own niAt eyiai mMtt m p 

o 

(assumed constant for each one) per unit initial undefomed volume. Each 

material can be quantified conveniently by an elastic modulus, a yield stress, 

an ultimate stress, a maximum strain, and strain rate parameters. These 

material values are obtained from standard tests conducted %rith a coupon 

of the structural material. The values a and a are defined in terms of 

o u 

engineering stress (also ci^-led the First Piola-Kirchhoff stress) a , 

E 

(Og - P/A^, 


where P « applied force and ■ initial cross sectional area) 


because the maximum engine^ing stress conveniently identifies the maximum 
load that the structure can bear (and this occurs immediately prior to 
necking). The strain are defined in terms of extensional 

strain, Ej^ (%rtiere ^ ° , £ « deformed length, and - initial length) 

because this is the conjugate to the engineering stress such that the total 
internal strain energy per unit initial volume U (proportional to the energy 
per unit mass 0 ) of the test coupon is: 

“ fi 


=1 ^ ■ k z ' XZ " T “ ■« (w) “ 


(4.1) 


The ultimate stress, is defined here as the maxin»™ stress attained 


during the tensile test (0 


P^/A^) and the corresponding extensional strain 


is called "the maximum strain, e^". Strain rate dependency is approximated 
by a strain-rate-affected yield stress, 0^(5), as defined from empirical data 
(Refs. 24 and 25) by: 

<r 


= <r(l +1 



(4.2) 


where D and p are eo^irical material constants and e is the strain rate. 

To complete the set of material variables (assuming isotropic materials 
for both the fragment and the container) are Poisson's ratio and the coeffi- 
cient of thermal expansion* Also, the impact-interaction is influenced by 
two material qunatities which "couple" both the fragment and the container 
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material together: the coefficient of friction and the coefficient of 

restitution, each of which is dimensionless. 

These **material property variabies" are sunnariz^ below. 


Material Property Variables 


Quantity 

Description 

Units 


Pc 

Material Density 

FT /L 



Elastic Moduli'^ 

CM 

< 


0 

Yield Stress 

T/T? 


c 



o 



u 

a 

Ultimate Stress 

t/-i} 

di 

c 


c 

u 



•M 


•- 




Ultimate St'^^in 

LA 

c 

0 




u 



-1 


0 

Strain Rate Constant 

T 


c 




p 

Strain Rate Constant 



c 




V 

Poisson's Ratio 



c 

a 

c 

Coefficient of Thermal Expansion 

0-^ 


pf 

Material Density 

FT A 



Elastic Modulus 

PA^ 


-’'f 

Yield Stress 

FA^ 


o 





Ultimate- Stress . . 

F/L? 


u 







c 

o 


Ultimate Strain 

LA ' ‘ 


m 



(9 

U 

Df . 

Strain Rate Constant 


- 

Pf 

Stratin Rate Constant 

— 


'"f ■ 

Poisson's Ratio 

— 




-1 



Coefficient of Thermal Expansion 

0 ^ 


u 

Coefficient of Friction 

“ 


e 

Coefficient of Restitution 

— 


TOTAL: 20 Variables 
















To define completely the impact interaction, several variables 

associated with the' surrounding ** environment" are needed* The local 

\ 

gravitational acceleration, the local temperature, and the local time frame 
must be included in the analysis,* and are listed below* 


Environmental Quantities 

Quantity 

Description 

Units 

9 

Acceleration of Gravity 

L/r^ 

T 

Temperatxire (above normal) 

0 

t 

Time 

T 

TOTAL: 3 Variables 


To this point, variables which characterize the containment structure 
and the attacJcing fragments have been discussed but no parameters have been 
chosen to describe the consequences or the results of the postulated fragment* 
impact attack. With respect to the containment structure itself, one must 
define the basic objective; two evident possibilities are: 

(1) to prevent fragment penetration — in this case one is interested 
in defining the dividing line between containment and non-contain- 
ment of the fragments or 

(2) to limit the maximum intrusion of the intact (non- ruptured) contain- 
ment structure into the surrounding space (to’ avoid inflicting 
da m a g e to critical components, controls, etc.) — in this case one 

might select the maximum deflection z of the c^tainer to be a 

c . 

parameter of interest. 

Also, it is conceivable (although xinlikely) that one would be interested in 
describing the consequences to the deformable fragment of its deplorable 
impact attack against the containment structure* In this case one might 
choose one or more parameters to describe the maximian deflection z^ of the 
fragment. Accordingly, one might add the following two quantities to the 
variable list: 
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in the -material- category, through are already diaensionleas, 
wlule ir^ through consist_of^os of container-to-fragment mass density, 
static yield stress, ultimate stress, and elastic n»dulus. The -environmental' 
category, includes gravitational to' centrifugal acceleration 

effects, and thermally-induced strain ratioed to ultimate strain for the 
container and the fragment. Included in the geometric set, through ir 
are the (dimensionless) number (n) of fragments (tacitly it irassumed thft 
the bursting rotor fragments into n equal-size pie-shaped bladed-disk frag- 
ments) as well as pertinent geometric ratios of characterizing dimensions of 
the containment structure and the fragment (s). Next, the hinediatic set, tt 
through may be regarded as indicating the severity of the impact atta^ 
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__ term may be interpreted as the ratio of 
fragment CG in time t to the waU thickness of th 
1T25 the CG velocity of the fragment i. 


the CG velocity Of the fragment is ratioed to"-^Il^ 

^ain rate constants and for the container and the frLenr"'”^ """ 

tively. Next, ir and tt„ may be reaarded , fragment, respec- 

26 27 regarded as representino the 

riTn-iixr: - - -- 

tip tiftto ,„«,P “ “• 

•Jiaensionless "eff^ parameters- ratioine ^h ' ’^29 “® 
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.Further if one assiaes (1) that the containment structure is at essentially 
a uniform tes^erature even though ^t teq>erature be "elevated" so that 
there are no thermally-induced stresses and (2) that whatever 
are carried out include the appropriate material property data (E , a , a 

®ob' °c' **c' ''c’ «l«vated tem^ature condition, one can%llminatr 

both the temperature T and from co n sid e ration as variables. With the 
e lim in ation of the 9 rigid-fragment "variables" and these latest 3 variables 
(g, T, and a^) , there are now only 21 x^iables left and these are expressed 
in terms of the 3 fundamental units* le^th, force, and time. Accordingly, 
one is now left with a "reduced set" ck (33-9-3) -3-18 independent dimension- 
less paramete r s (S) , as renumbered and^: listed below; 


REDOCED SET OF DIMENSIONLESS PARAMETERS 
(Risid Fragment, Near Earth, No Teo^rature Gradients) 


Material 

Geometric 

Kin««atic 

®1 * ^nc 

83-0 

®15 ■ ‘VcG^’/^'c 

8, - V 
2 c 

89 - 

®16 ■ 

S3 - Pc 

\o ■ V'a:**,' 

hi ■ '“c 

84 = u 

^u’V^'f 1 

hs • ‘cA 

= « 

S 12 “ V^c 


Bg p^/p^ 

®13 " **c^c 


“ V'^c 

0 




Einally, it should not^ that this is a possible but a non-unique set 
of characterizing dimensionless pa^et^s. Other canbiaations of these (or 
equivalent related) variables could be constructed to give an equally-valid 
set of 18 dimensionless parameters. 
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Next in Subsections 4.2 and 4.3, related more-restricted categories of 
problems will be discussed involving certain experiments and/or numerical 

predictions. Also, an alternate set of dimensionless parameters will be 
considered. 

— Dimensional Cons iderations Applied to the NAPTC Rotor-Burst^ 

Containment Experiaentis 

As a part of the NASA study on Rotor Burst Protection, the Naval Air 
Propulsion Test Center (NAPTC) has conducted a series! of experiments (13,37, 
for example! in which actual aircraft engine rotors have been rotated at high 
speed, caused to burst into various fragment patterns/ and the subsequent 
fragment impact against containment rings has been studied. These experi- 
ments were corducted in order (1) to investigate the associated impact- 
interaction phenomena for containment structure impacted by typical complex 
deformable fragments, (2) to obtain an empirical data base to aid in 
containment-structure design, and (3) to obtain data which would be useful: 

(a) in guiding the development of and (b) in evaluating numerical methods 
for predicting the responses of bodies involved in fragment/structure impact, 
interaction, and response. 

It is convenient for discussion purposes to consider two separate groups 
of the NAPTC experiments. In one group-simple- containment rings of a given 
material were used in an attempt to determine the containment threshold for 
rotor burst attack from a GE TS8 turbine rotor spinning at 20,000 rpm, for 
cases in which the rotor was caused to rupture into 2, "3, or 6 equal-size 
pie-shaped bladed-disk fragments. _ The second group consisted of similar 
fflcperiments but involved a different and larger turbine rotor (a PSW J65) 

with an 8500 rpm rotor-burst condition. Each group is considered separately 
iii the following. 

4.2.1 T58 Rotor Burst Containment . y 

The T58 . turbine rotor employed in these containment, studies as depicted 
in Fig. 37 has a 14-in diameter, about a 1-in axial lengei. 51 blades — each 
with a 3.5-in length, weighs 10.8 pounds, and was modified to rupture into 
2, 3, or 6 equal-size fragments at a nominal 20,000 rpm. 





In these containment experiments, there was a single source of 
— the T58 rotor ~ hence, the fragment material parameters were not varied; 
also, the basic geometric parameters of the fragment source were not varied, 
only the number n of ^l^aLze fragments used was varied, and the geometric 
parameters associated with each are Implied when one specifies the value of 
n. With respect to the containment structure , it was intended that a single 
containment material (4130 cast steel) be employed. The only variables 
employed were the axial length (or width) w^ and the radial thickness h of 
the containment structure, the latter quantity being varied as a means of 
defining experimentally the containment threshold. 

In view of these facts, it is convenient and useful to sunmarize the 
fixed quantities and the variables needed to characterize the containment 
ring and the fragments. Further, becauje of the restricted set of variables , 
involved, it is useful to select a slightly different descriptive set of 
quantities for the fragment attack from the "general set" considered in 
Subsection 4.1. since the fragment geometry is unchanging except as implied 
by the number n of fragments, it is useful to describe each fragment by 
= mass of the fragment 

I^ ■ mass moment of inertia of the fragment about its CG 

Uf » angular velocity of the fragment (and rotor) at rupture 

Vj = translational velocity at the CG of the fragment at release 

'cG “ distance from the axis of rotation to the CG of the fragment 

Automatically implied is the distance from the fragment CG to the blade tip. 
Accordingly, listed below are the characterizing quantities which remained 
fixed-and those that were varied in this group: of NAPTC experiments; 


CONTADmENT RING 

FRAQffiNT ■ 

Fixed Quantities 

Fixed puantities 

Material 

Material 

Inner Surface Radius , r 

“f' '^f 

Variables 

Variables (Only as Implied bv m 

Radial Thickness# h 

m..# I - ]r ' 

c 

Axial Width# w 

f f' CG' ^f 

. c 
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(4.4) 


(^o) — i (hc/r J J f n) 

result variables 

where W denotes the weight of the containment ring. Alternatively, if the 
c 

actual threshold ruptxire condition is being sought and described (thereby 

reolacino (£ ) )« one can represent the containment threshold by the folloi#- 

c max ^ "■ — — — — — 

ing functional dimensionless description : 



Although this characterization is given entirely in terms of dimensionless 

quantities, it may be more useful and clear to describe the threshold contain* 

ment condition in terms of only attack kinetic energy (WS)^ and ring weight 

instead of (KE) /(W r) since both (KE) and r are held fixed. Other 
o c o 

'* dimensionless-result displays" will be presented for similar reasons. 

Summarized for convenience in Table 5 are the parameter values identify- 
ing the T58 turbine-rotor-burst fragment attack for n»2,3, and 6 *-agments. 
Table 6 contains a summary of the NAPTC tests of the T58 turbine rotor bursts 
against 41-30 cast steel containment rings. Note that while the ring thickness 

h was varied. in tests using n*2,3, or 6 fragments, the effects of ring axial 
c 

width w were studied to a meaningful extent only for. the n»3 case. Finally, 
c 

note that the 4130 cast steel material used for these containment rings came 
from two different suppliers and the associated mechanical properties turned 
out to be somewhat different (Table 6). This fact, should be kept in mind in 
assessing the resulting experimental data; most tests involved the NF material 
but a portion of the tests (that is, some of the tests for n®3 only) involved 
the "different" ACIPCO material. 

In the cited NAPTC tests there was no determination of the circumferential 
strain on the ring (either extreme peak transient or extreme permanent); 
however, the permanent deformed ring configuration was photographed for 

'*'lf desired, one could use (W^r)/(KE)^ rather than (KE)^/(W^r) to obtain a 
small rather than a large dimensionless number. 
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doomentatlon. The principal result from'eadi of these tests was that the 
ring either contained (c) or dfd hot contain (HC) the attacking fragnents. 

Thus, unless conditions for C and NC are very close to each other, it is not 
possible to display results such as those indicated functionally by Eq. 4.4 
or by Eq. 4.5. Instead, one could evaluate and display the ratio of the total 
t ranslational pre-impact kinetic energy (KE)^^ to the weight of the contain- 
ment ring (called herein the STPE - specific translational fragment kinetic 
energy) as a function of the number of n of attacking fragments while holding 
constant the ring axial width w^, in this category ofi teats, the ring thick- 
ness h^ would be varied, and would automatically isq>ly a value for the ring 
weight W^. A plot of this type is shown in Pig. 39 for the fixed value of 
ring width w^ - 1.0 in*; the symbol "0" indicates that the fragments were 
contained while the symbol "X" means that the fragments were not contained. 

Note that Pig. 39 is similar to Pig. 6 of Ref. 37 where the ratio of the 
total fragment kinetic energy (KE)^ to the ring weight - called SCFE 
(specific contained fragment energy) — is shown. In this report the STPE is 
chosen as a basis for discussion since rotor burst containment ring performance 
depends almost entirely upon the value of the translational portion of the 
fragment kinetic energy attack; the rotational portion is of comparatively 
insignificant importance. A "curve" has been faired through these data to 
identify (or estimate) roughly the dividing line between containment and 
non-containment; additional data would, of course, improve the reliability of 
this faired curve. One should expect this faired curve to approach an STPE 


asymptote as the number of equal-size fragments is increased indefinitely; 
thxs «rould correspond to a uniform membrane strain state in the ring emd, 
hence, would ^^able the ring to absorb a maximum amount of strain energy 
before "simultaneous outer-surface n^turing" would occur at all circumferen- 
tial station^ of the ring. - On the other hand, the least noaber (n»2) of 
equal-size fragm^ts produces severe local bending in addition to the membrane 
straining; the result is’ that this locally-seyere bending-membrane outer- 
surface extensional strain reaches rupture levels at the least total 
translational kinetic energy attack value. At all other values of n. 


* 


There are insufficient data to permit making similar plots for 
values of* containment ring width w . 


other fixed 
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this threshold eontaiiment curve lies between the n-2 value and the afore- 
men^oned asyn^tote. Also> only tor n«2 is disk sharp~edge gouging of the 
ring evident. . 

*^*^**^^*•9 to Pig. 38r the anounts of translational rotational kinetic 
energy possessed by fragments of varying sizes originating from the sane rotor 
are shown. The three-fragment in^ct situation has more total triuislational 
kinetic energy than does the two-fragment attack. On an individual basis, the 
impact of one tri-hub fragment is more severe than the in^ct of one , bi-hub 
fragment; however, the fragment-ring interaction of the three fragment attack 
lessens the severity of the impact- induced response. As the total number of 
rotor fragments approaches infinity, the deformation of the containment ring 
^roaches that of pure membrane stretching. As the number of fragments 
present decreases, the severity of local bending deformation increases. This 
is readily visualized in Fig. 40 (from Ref. 37) where a comparison is made of 
the final deformed shapes of identi c al containment rings subjected to either 
two- or three- fragment la^ct. The two- fragment lz 5 )act produces fewer regions 
where bending is present, but the level of bending will tend to be greater in 
these regions than in corresponding regions for rings impacted by more frag- 
ments. The hypothesis offered here is that a 1/3-rotor fragment — 

more translational kinetic energy than a 1/2-rotor fragment, but the impact- 
induced ring response (taken as a collective action) appears to be more severe 
for the- two-fragment attack. 

A limited experimental study was made also of the effect of various ring 
widths w^ (or the ratio on T58 turbine rotor fragment contai^ent — 

only for the 3-fragment case. Fig\ire 41 shows the STFE as a function of con- 
tainment ring width w^ for the~3- fragment case of T58 turbine rotor burst; 

These <tota indicate an .optimum "energy-coping capability" at a ring width w 
of 1.00 in, which is equal to the axial projected width of the attacking 
fragments. While w^ = 0. 5 in and w^ ■ 1.0, in clearly. involve 2-d structural 
responses of these rings, the for^ner exhibits larger bendi^- induced strains 
compared with the latter , and hence reaches a ^rupture strcu.n level at a 
smaller attack specific translational kinetic energy veilue than the latter. 

On the other hand, having w^ significantly larger than w^ will lead to 3-d 
structural response behavior in the immediate vicinity of the points of impact 
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used. For the present rough determination of 5^* the NP material is appro^‘ 
uted as being "rigid, linearly strain hardening" to n^turei hence. 
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Then 0 is found by multiplying this value for U by the initial volume of 
c 

the coa^rising the containment ring. 

The soUd curves shown on Pig. 42 for n-2, 3, and 6 were specified as 
follows. Since the total fragment mass is known and the nominal rupture rpm 
is given as 20,;000, the total translational kinetic energy is known and is 
gon sta n » for each value of m for a given containment ring, the ring mass 


and 0 (for the NF material) are both determined — thus defining the soUd 

curves associated with the NP material for eau: value of n. Similarly, the 

dashed curve was established for n»3 for the AClPC^supplied 4130 east steel 

material; for this material, analogously the following very different value 

for U was iised: 
c 



+ 38,Z70 


in 


-lb 


In 


which means that the ACIPCO material is much tougher than the NP material. 

It should be noted that the curves shown on Pig. 42 indicate the nominal 
or intended rotor-burst condition of 20,000 rpm. In actuality rotor rupture 
occurred at a slightly different rpm in a given test. For a given ordinate 
value, rotor rupture actually occurred at a higher (lower) rpm than nominal 
for those points lying to the right (left) of a given curve. Thus, one can 
see readily how near the nominal condition that the actual rotor bursts 
occurred. 

An examination of Fig. 42 reveals the following two trends: 

(1) There is a ratio of ring mass to fragment mass below which there 
are few or no containments. This level is approximately 0.7. 

(2) There 'is a ratio of total fragment tremslational energy to maximum 
energy absorbed fay the containment ring above which there are few 
or no containments. 
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The first trend cited above occurs at nearly the same mass ratio for all 
tteee values of fragment number n. Evidently, as the mass of the contaianent 
ring decreases in comparison with the fragment mass, a point is reached where 
the deformations surrounding the "impact point" dominate the Impact- interaction 
and lead to a local tensile failure without greatly deforming the surrounding 
oontalnment structure. There is no fixed energy ratio for the second trend, 
but a different energy level prevails for each fragment configuration. This 
ia because the magnitude of the total translational energy stored in the 
fragments is dependent upon the configuration (or value of n) . Prom Pig. 33 
it xs noted that an individual tri-hub fragment possesses the greatest trans- 
lational kinetic energy, but collectively the six-fragment configuration 
possesses the largest total .translational kinetic energy (whereas an even 
larger total translational kinetic energy level is attained by a nine-fragment 
Impact configuration). A hypothetical containment limit as n- would be at 
an energy ratio of 1.0 if 0^ were known accurately, because this is the point 
at whxch xt xs assumed that the maximum energy can be absorbed in just reach- 

xng tensile failure simultaneously at all circumferential stations in the 
contunmen^ structure . 

This containment limit is based on the absorption of translational 
kinetic (impact) energy by the containment rings through membrane stretching. 
However, if ben^ deformations occur at a non-negUgible level, the contain- 
ment Ixmit would Shift to a value less than 1.0 because an equal volume of 
material would absorb twice as much energy by pure membrane stretching than 

xt wxll by pure bending before reaching the -local jupture" condition. The 

extent of bending deformation for the-2, 3. and 6-fragment rotor burst was 

mentioned earlier in this section and is borne out iri Pig .’- 42 as a shifting ‘ 
containment limit. 

. of aispl., 1 ., th, Hwre 

on TS, turMn. tri-hnb (n-3) ronor burnt =ont.lnn»„t „n unnful nnd are she™ 
en frgn. 43 and 44„only for the »F natarial for clarity. Fioure 43 shown 

(KE)*"^^° r \ to rxng inner-surface radius r; (h /r) versus 

(KE)ot/°c of the width ratio w^/w^. Figure 44 shLs the 

ratxo of containment ring thickness h^ to ring width versus the energy 
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( 2 ) 


(3) 


ratio (KE)^^/U^ for fixed values of the width ratio w^ Confining 
attention to tri-hub bursts (n-3) , an examination of Figs. 42 , 43, and 44 
reveals the quantitative and the expected quaUtaUve reUtlonshlp between 
the thickness of the ring and the -containment potential- of that structure 

TO achieve containment by the HP 4130 cast steel material (for the present 
limited data only) ; 

(U *2 toilet., ttot . rlt, ,t i.et 0.7 ttoe to. tot.1 

fragment mass should be used. 

Figure 43 shows that a ring thickness h^ at least 0.08 times the 
ring inner-surface radius should be chosen, and/or 

A ring width at least as great as but not much greater than the 
fragment width should be selected, according to Figs. 41 and 44, 

however, the data are too sparse to permit making a more precise 
choice. 

These criteria for estimating the containment bounds do not guarantee 
containment (as is indicated in Figs. 42-44) , however, all of the rings which 
contained fragments that were released from rotors at the design rp „ „eet 
all of the above-mentioned criteria, only test 177 fails to meet these bounds 
but the rotor failed at an rpm significanUy lower than the design speed. 
There appears to be a minimum width, thickness, and total mass for the 

containment ring of given material in order to insure a high probability 
of containment. “ 

These observations are also discernible from similar HAPTC plots in 

“ clearly apE^aat_as_in_F-igs^-44. i„ addition, the— 
nondimensional plots in Fig. 42 readily separate differences based on 
material propertiesr show the relation of the data points to an "ideal curve" 
for that material, impacted by a fragment released at the design rpm. and 
indicate the pertinent gemetfic variables. 

Th. toove brief Ultotratioa of the os. of hoodtoen- 

Stooel eoalysis for dispUylo, toe t»)or f.otois l„fi»,„ei„, o,, cootsltoeo. 
or a rotor fra,«ent if ail appropriate data had been obtained. This process 
can be applied to any set of experlitontal data in order to qoantlfy mjre 
clearly toe tofortant interactions and paranetrio effects, particularly if 
to. disenslonless combinations are chosen oito some prior knooled,, of to. 
principal physical processes present. 
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4. 2. 2 J65 Rotor Burst C ontaiment 

The Curtlss-Wright J65 turbine rotor eapK^ed in this group of NAPTC 
ei^erimental fragnent-containnent studies has a 30. 64- in diameter, a 1.25- in 
axial length, 72 blades, a weight of 127.75 pounds, and was to 

rupture into 2, 3, or 6 equal-size bladed-disk fragments at a 8500 rpm 

The-pre-impact fragment data characterizing these tests are given in Table 7. 
To date, the wa^er of fragment-containment tests using this "large* rotor is 
** similized in Table 8 where the specific tests are idratified by 
test nunber, the number of equal-size fragments^ and the widtii w of the 

4130 cast steel containment ring used in each test. Hote that all of these 
tests (except two) utiUzed rings having w^-w^-1.25 inclt also, e^ of these 
rings had an inner-surface radius r-15.82 inches. Listed also in Table 8 are 
the aechanical property data (37) for the ACIPCO-s<«>plled 4130 east steel 
batch from which these containment rings were made. 

Since the geometry of the J65 rotor is different from that of the 
previously-discussed T58 rotor (that is, a simple scale factor multiplying 
the geometry of the T58 rotor does not result in giving the dimensions of 
the J65) , one can not make direct eoe^arisons between these two sets of test 
results in a single straightforward dimensionless variable sense, in addition, 
the rotor-blade and disk materials for these two rotors are different. Also, 
the available data are too sparse to make any critical comparisons between 
these two sets of data- (T58 and J65) . Hence, the present J65 data are dis- 
played in plots similar to some of those discussed earUer for the TS8, as 
follows: 


Figure No. " Quantities ’ 

45 STFE vs. n ■ 

46 “/(nm^) vs.- (KE)^^/0^ for fixed n 

vs. (KE)^^/D^ foi: fixed n . ; 

For these J65 slower, more-massive fragments and more-massive containment 
rings (than for the TS8 tests) , the value used was estimated by assuming 


a (smallar) effective strain rate of 500 sac**^ for the "mild steel**^ 

ACIPOO 4130 cast steel indicated in Table 8.: Analogously*, 

= 60 ,cee (i + j-^p) = ISf,oeo ^.i 

ze o»hs<l,oce + - "'^ I « 14 320 
'AC~^S I * J ' 

Then, 0 is found by nultiplying this value for U by the initial volume of 
c c 

material in the contaiment ring. 

Because of the relative paucity of the present contaiment ring test 
: data for the J65, one can make only a few tentative observations, parallel- 
; ing those made concerning the TS8 test data: 

(1) There is a minimum ratio of container thickness to ring inn«^r 
radius (h^/r) below which no contaiments are expected. This value 
appears to be 0.27 for the c-W J6S rotor fragments (Fig. 47). 

(2) A ratio of container width to fragment width of one (w • 1.0) 
appears to be optimal for contaiment purposes. 

(3) There is a maximum value of trams lational fragment energy to 

maxim um energy absorption capability of the container ()S) /D 

c 

above %rhich no containments can be expected* This value vuies 
with the number of fragments involved in the impact but in all 
cases must be less than 1.0 (Fig. 46) since this is the theoretical 
limit where pure membrane behavior (maximum energy absorption) occurs. 
It is apparent from the brief study of the C-W J65 rotor impact data and the 
previously-reviewed 6E TS8 rotor impact data than an "optimal experimental 
process" can be carried out to find' a ‘ contaiment threshold. First, a contain- 
ment ring equal in width to the fragment axial length should be chosen. Then ‘ 
the ring thictoess should be varied increasing from a value at which the mass 
ratio is one (since this is the theoretical, lower limit of contaiment ring 
weight — pure membrane behavior) until a' containment is achieved. This process' 
includes all three of the above-mentioned impact-interaction observations and 
should lead rapidly in experiments to the Selection of a "light-weight" fragment 
container. 

Note that under these assumptions (admittedly rough), this ACIPCO material 
under the J65 test conditions appears to be less tough than for the T58 
cases with either NF or ACIPCO material. It would act even les^ tough if 
the "effective strain rate" were sm^ler. See Tables 6 and 8. 
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con^iamMt stroctureo is concerned. Aioost withoiit ejcception i' these oladed-^iisk 
fragments act as "readily deformable" , during their early stages of impact- 
interaction, and deceleration before the "effective fragment" reaches a more 
compact, bird, and stiff state; the containment ring becomes loaded over a 
"fairly large" circumferential region by each fragment and subsequent structural 
response ensues. 

Some supplementary information to substantiate the expectation that 
shear plug failures are unlikely for the postulated types of fragments and 
range of normal-component iiig)act velocities may be of interest from impact 
studies reported, for example, by Rinehart (31], Hagg and Sankey [43], Recht 
and Zpson 1441, and Letl^y and Skidmore 145] . Both tensile rxipture and 
shear-plug failure modes are encountered in the "sub-ballistic velocity 
range", and are discussed in those studies. It has been known 131J for a 
very long time that a "plugging failure" will occur very early in the ia^act- 
interaction-response process, whereas tensile rupture occurs very late in 
this process when the gross structural response of the slowly- responding 
structure "reaches a peak", roughly speaking. The following discussion, 
however, will be confined to shear plug feu.lure and criteria for its 
occurrence, and follows the Hagg and Seuikey 1431 description. Of particular 
interest is the normal conq>onent impact velocity v^ below which no shear plug 
failure can occur and is given, according to simplified theory, by the follow- 
ing inequality: 

PT -h cr^€^\/ ( 4 . 6 ) 


If v^ is large enough to violate this inequaUty, s. 'tr plug production is 
possible. Expression 4.6 means that plugging failure will occur if-the energy 
abso'tb<a-By the target structure’s "impact-affected" region is les^than the- 
energjQLo^t during an inelastic impact (e=0) , where 

K 




^n experimen^lly-determined factor related to the shape' of. 

* " ^ X 1 -i 

‘ -the impacting face of the missile < J K < ^) 

; is the thickness of the target structure (in) 
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P is the shearing perimeter of the contact region between the 

impacting fragment- and' the target structure (in) • 

T is the dynamic ultimate ' shear strength of the target (psi) 

®u^m represents approximately the energy absorption capacity of the 
structural material to failure, per unit initial volume 

is the volume of the "impact-affected" portion of the structure 
is the mass of the fragment 

®c the "impact-affected" portion of the structure 

Vj is the impact velocity component perpendicular to the impacted 
surface 

f(mj,m^) for an inelastic impact t23,43J is given by m /(m +m ) 

c f c 

Restricting attention for present illustrative purposes to a 2-d containment 
ring or beam for which the width of the impacting fragment is equal to 
the width w^ of the containment structure, one may write Eq. 4.6 (as an 
equality) as 


K 






W 

c 


hi = 


-i-m, 
2 -f 




h i 


(4.7) 




+ /c'^c 




where 0^ is the mass per unit initial volume of the containment structure and 
I is an as-yet-undetermined impact-affected spanwise length of the structure. 
Collecting terms and rewriting, Eq. 4.7 becomes 


Bli + C].= o 

where -x ' 

A = 2 i 

B = 2 ^ 


(4.8) 


(4.8a) 


(4.3b) 


32 















<4.8c) 


C = "c / ^ ">f '>/ I 

The-soIixti<»^of Bq[* 4.8-is-gtven-^: 

(1) hH}i Trivial Solation 
or (2) Ah^4fih+C*0 or 


>•'7 




(4.9) 


B 

Since — > 0, a pfaysically-valid solntlon (h ^ 0) is possible oiay for C-0: 




m. 



O 


(4.9a) 


Thus, the plugging threshold iapact velocity is indq>endent of both T- and t,~ 


and is given by 






(4.9b) 


Applying the Table 6 values for the NF-si^lied 4130 cast steel, the 
plugging-threshold perpendioilar-iapact velocity v^ is estiaated by 
single theory to be: 


Z(iiifOooX-os) 

(2.S3)/(3S4,.0^) 


= 4^22 


in /sec. 


if the static properties are used or 

''f = 

if the previously-discussed estimated (iynamic properties are es^loyed. The 
a^licability of this analysis reqpiires the missile to be essentially non- 
defonaable; thus one can not apply this vhen considering the "readily 
deformable" blades of T58 bladed-disJc fragments. How eve r, if one ct^iders 


(.z83)/C3e&.o=i) 


I 

IT 


= 7277 m/sec. 
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to Steel->Sphere Tmpac^ 

Fxrst, it should be recalled the subject MIT-ASRL experiments involved 
perpendicular impact of a 1-00-in diameter steel sphere at the midspan-mxdwidto 
location of 6061-T6S1 aluminus -cpntainaent- beams of thickness h-0. 100-in, 
width w^-l.SO in, and span t-8.00 in with both ends ideally clanged. In these 
experiments only the initial impact velocity of the steel sphere was varied 
in order to produce various degrees of peak structural response including 
nq>ture; as reported in Ref. 16, experiments and test specimens CB-9, CB-13, 
CS-18, C3-16 and CB-14 with steel sphere initial impact velocities ranging 
from 1900 t« 3075 in/sec permitted determining the initial impact velocity 
(or initial fragment kinetic energy) for threshold containsent. The pertinent 
geometric, mass, and material property data are sunaarired in Table 9. Under 





tiMM fra^Mt 9«cMtcy and baas gmemmtxy iapaet omditions, tha stroctiiral 
raaponaa of tha baaa waa of 3 -d charactar at and closa to tha point of 
lapact but appaarad to ba of 2-d charactar owar tha raaaliiluj (approKiaataly 
^0 P**^ cant voIuMwiaa)--portion-of-tha-atructura. Rvptura at tha p ld tr a n 
iapaet station appeared to ba govamad by tha 3-d structural re^onsa 
behavior produced there; this rupture was of tensile character rather 
a shear-plug failure.* 

With respect to a dlaensionlass-par a aa ter , d u i r a ctar ixation of this sphere- 
beaa iapact-intaraction and response problea, 6 im aight eaploy the set of 
diaensionlass paraaetars discussed in Subsection 4 . 1.2 for containent- 
structura response to riqid-fragaent attack — since tha steel sphere 
reasonably quaUfies as a ccaparativaly rigid fragaent in the present instance. 


Tha 18 paraaetars 8 


^18 Subsection 4 . 1.2 need to be relabeled 


material parameters 8,-e 

1 me 


as appropriate for the present sphere-beaa iapaet situation. First, since 
the saae bean material was used throughout this test series, the diaensionless 

®2“^c' ®3“^c' ®4“^' W^f' 

remain valid but are not variables . Meat, the geooetric parameters nay be ° 

recast by using (a) the sphere diameter d to r^lace the previous (r ♦£ ) 

OS f 

and/or and (b) the beaa span t to replace r ~ to relabeled 

parameters 83 ... 8^^ as 63-n. B^-iD/d. 6^^-w^/d, 813-h/w^. 

834“^^* here again all but the last paraawter remain meaningful but none 
are variables. Finally, the "kinematic" set of variables can be recast 
by using the normal-to-the-surface impact velocity v of the non-rotating 


sphere to replace and thereby rewriting 

B^5-(vt)/h, 6^3-v/(h3^), 8^7-(a^/e^)/(p^ v^). »*ere 6^3 is an 


® 18 « 


"effects parameter" — the ratio of the naximia deflection to the bent thick- 
ness. 

Note that the only independent variable is the ini tial i^act velocity v. 
Bence, one may express 8j^3, for example, as a function of the dimensionless 
variables 8^^^, 8^3, and B^^t 


1 /: 


• 20 e 

Note that the previously-cited criterion for shear plug failure v^ > (— H— ij 

and the pertinent data cited in Table 9 indicates that for shear plug failure, 
the fragment velocity v^ should be > 8400 in/sec — well beyond the range of 
test conditions needed to produce threshold rupture of the *<»»» (ifij , 


85 




fin) 






(4.10) 


Itet. i( o,» i. «, in tt. tl« of oecnnonc of th. -..-n 

diWl.o-.nt Sj, c- b. nli.an.t,l. ^ 

-d on. fr^t nnt-Ui hn- b-n bou, 8., «« 8„ dnfin.8 .b- 

» 1. w-lflw, h—. on. o- tw~»i S„ - being tb. nnlj .Igniflo-t 
variable, thereby obtaining 

^ - ^(fn) <«*10a) 

for a final cause-and-effect display. 

Alt«™.ti,.l,, ,.tb- tbnn u.ing bw«t-lnduc.d 

effeot.. on. -gbt .lect to -piny th. .-1-. longitodinnl m>p.r--rf.o. 

“•‘"’““‘I' -y fit. - . fnnotlo. of tt. only 
Significant variable 8^^: 


= 


»*«x 




,)/(p^ v“), it is 


4 

I-t-d-of diwuyin, - . foppoi^ ^ o ^ ^ _ ^ 

c^ny-t-t to twuo. S„ by th. „tio of th. inlti J“hi“ti.- L-gy 
ofthefr—nt ^ »,<0, y*./2 - a« -U-Uz- -iS. ' 

W-wtion Ind- for th. c-tnl— t .wnotnr.- „ „ di,oo.,«, i„ 

4.2.1, %diere 0 *0 tises th» — i 

c c “*® initial volume of containnent 

material and 


(4.11). 


i;=€ 

c •% 


[• 


<^ + 


£5 


“energy capacity” per 
^uut inxtial volume 


(4.12) 


Hence, one may %n:ite 


is the volume of the fragment (steel sphere) material. 
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-m#x 


= % 




(4.13) 


In the sphere-beam impact ejqteriments under discussion (161 , transient 
strain measurements were oude at various, spanwise stations for the upper 
and/or the lower surface of the beam. Steel-sphere i^ct occurred at 
midspan (x-0) on the lower surface. Although the spatially-maxijmmi strain 
occurred at the upper surface at or near station x^, transient strain 
measurements were not made successfuUy there; however, successful measure- 
ments of upper-surface transient strains were made, for exaog>le, at 
|x| ■ 0.6 in and |x| - 1.50 in on the beams cited in Table 9. Accordingly, 
the measured peak upper-surface longitudinal strains (denoted as e for 
convenience) at stations |x{ » 0.6 in and |x| - 1. SO in are shown in Pigs. 48a 


the test/specimen# fragment 

velocity (KE)^, 

and <KE)^/tJg are tabulated below: 

Test 

and 

Fragment 

Velocity 

(KE)^ 

(KE)^ 

~ 


Specimen 

V (in/sec) 

in- lb 

c 

Remarks 

CB-9 

1896 

686 

.0657 

Small permanent deflection 

C3-13 

2490 

1183 

.1133 

Moderate permanent deflection 

C3-18 

2794 

1489 

.1426 

Large permanent deflection 

C3-16 

2868 

1569 

.1503 

Slightly beyond threshold con- 
tainment; specimen ruptxired 

CB-14 

3075 

1804 

.1728 

Well beyond threshold con- 
tainment; specimen ruptured 
very early 

where 0 * 6 x 
c c 

(beam volume) 

= 8700 

(1.20 in^) * 10,440 in- lb. Note in 


Figs. 48a and 48b that the measured peak strain at each of these stations 
increases with increasing (KE)^^ — for specimens C3-9, CB-13, and CB-18. 
However, under more severe impact, specimen rupture occurs and thereby 
"releases the loading" which results in smaller peak strains at "remote" 
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•tatioas |x| - 0.6 ia and |x| - 1.50 in» thia ia particalarly pronrancad for 
apoc inan CB-14 which nptared v«ry aarly aftar liaitial ii^aet. 

-By-analo9y-with-tha aannar of praaaotiag tha RAPTC data in 4.2 

to idantify thraahold containaant, Pig. 49 dapicta tha praaaat lOT-ASItt, aphara- 
bean, ia^aet taat reaulta in tana of tha ratio of bam (containant atruetura) 
naaa to fragaant aaaa varaua (1S)^/D^} fragaant eontaixawnt ia r^raaantad by 
tha ayal^l "O" and non-conf a i naant by "x*. 

Bacauaa fragaant valocity waa tha only indapandant variabla in tha taf . 

16 aphara^aa iapact axpariaanta, thia aaaantially coa«>lataa tha diaplay of 
tha par^nant diaanaionlaaa data, ttaxt, nuaarical paranatric 2-d atudiaa which 
ainulta^ in a reatricted aanae tha aphera-baaa isg>act problea will be 
diacuaaed. 

4.3.2 Paraaetric WMBerical Pr«»«<< ctions 

It is convenient to use . the sphere-beam impact problem to illustrate 
the use of the CIVH-JET 4B code for parametric numerical predictions 
of (a) the small number of variables which can be eng>loyed usefully in 
studying this well-defined pr^lem and (b) the relatively economical confuta- 
tions required to obtain pertinent results. Recall, however, that the 
ClVM-JEr 48 code aif>Ues strictly to problems involving 2-d behavior and 
structural response. Therefore, in applying this code to simulate the 
sphere-beam impact problem, and particul 2 u:ly when one varies certain 
9 ®o®®tric and other problem parameters, changes in one parameter may <"«piy 
changes in a number of other parameters. This matter oust be kept in 
in defining the characterizing non-dimensional variables and qtiantities 

for the following four illustrative cases considered here: 

Case A: Only the diameter of the isfucting fragment is changed in a 

of calculations (i.e., d-1.00, 1.25, and 1.75 in) but the fragment 
mass And the init ia l impact velocity are from CB— 18 

sphere-beam ii^ct condition. ' 

Case B: The fragment mass and initial intact velocity are fixed at the CB-18 

values and d*1.00 in; the mass and the scan of the al«inca 
are fixed at the CB-18 values. Vuied in a sequence of «»lCTilatl on 9 
is the ratio beam width to beam thickness to assess this 

effect on the peak ii^ct- induced strain. 
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Casa Ct Tha fragment mass and initial ia^act velocity are fixed at the CB-18 
values, d«1.00 in, and the span and width of the alu&inum beam are 
fixed at i-a.op in and w^«1.50 in. Varied in a sequence of calcula- 
tions is the thickness h*of the beam to assess this effect on the 
peak ia^ct-induced strain in the 2-d beam structure. 

Case D: The fragment mass and initial in^ct velocity are fixed at the CB-18 

value, d«1.00 in, and the span and the thickness of the aluminum 

beam are fixed at 1*8.00 in and h*0.100 in. Varied in a sequence 

of calculations is the uridth w of the beam to assess this effect 

c 

on the peak i^act-induced strain in the 2-d beam structure. 

In all of titese calculations, the beam with both end s ideally clasq>ed was 
modeled, in accordance with the modeling guidelines established in Section 2, 
by 43 equal-length finite elements with 4 DOF/node, 3 spanwise aind 4 depthwise 
Gaussian integration stations per element, and a fixed beam span of 1*8.00 in* 
Also, the uniaxial stress-strain behavior of the beam material was represented 
via the mechanical sublayer model by the following Fit A stress-strain pairs: 
(C,e) * 41,000 psi, 0.0041; 45,000 psi, 0.0120; and 53,000 psi, 0.1000. The 
material %ras regarded as behaving in an EL-SH fashion — insensitive to strain 
rate; hence, the material strain-rate parameters O and p are deleted. The 
beam initial mass per unit volume was taken as 0.25384 x 10 (lb-sec }/in . 
The impacts are treated as frictionless and locally perfectly elastic 
(coefficient of restitution e*l) . ^ 

Before discussdlng the numerical pr^ctions obtained for Cases A, B, C, 
and O it may be useful to review in a summarized fashion the pertinent 
dimensionless parameters and variables. First, returning for a moment to 
the eag>erimental sphere-beam inpact results discussed in Subsection 4.3.1, 
the previously-discussed parameters for that 3-d structural response problem 
are listed for convenient reference and comparison in Table 10 where the set 
of 18 "usually significant p£irameters" has been reduced to 17 since the 
attacking fragment has only one significant dimension: its diameter for 

that experimental situation wherein only the peak stradn (or ^reshold 
rtpture) was of interest, only one significant variable was present %diile 
all other parameters remained at fixed values. 
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Show U.O la TttU 10 m th. t o t iei r o MlT^; aimMlcalm 

W„ ... «- « • 9 >-^ i-PWt pwbl- » ^ - 

tt. M CI.1I-JB 4S «p.t« <»d. te -oh Jt th. — . », B, C, -« B.r.r 

.1-a.tlo., «. B-S— -W b. . rlold elrc— — 

0 -d -4th -U-r — h— . -h—. H— , U « -ri- — 
p... 4i-.ur Wil. ht«i- tBo ttttmmt — oo— nt (w i. oo- ». B. c 
D), ow — ; wrtoioo thlo M h -«9 wo-plloh- by (*> p-oeilhliw th. 

-4th to -toh th. h— -40. -4 t— rl.4 th. fr-wot — • 

p., o-t -I— 0, to chw,. -cor41o,ly or (h) —i- 0, »i«- •«* 

1 ,., to chug, wr.gci.t-y — .tlp-.tiog .t th. w O- 
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W il.u,io— . p— o-t. Wich r— o fl-4 W th... Wich wry. 1. 
.tctOuc. -th th. pr—iog 41—00. th. folio— — th. ouo-by-ou. 
variables and the result of primary interest: 

Variables 


Casa 

A 

Primary 

6, - 1/d 

Secondary 

Primarv Result 
®18 * ^mar 

B 

®11 ■ >^0 

»10 - 
“l2 ■ 

*18 ■ t«o. 

C 

Sij ■ VI 

»io - 
*11 ■ *^c 

Sl8 • ^max 

D 


*u ■ "o^f 

fi « e 
*^18 na* 

results 

for each of these 

four cases are 

discussed separately 
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shom in Fig. 50 for lowor-surface ii^act oceorrlzig at the midapan 
station x«0i a calculation time step size At of 0.5 microsecond %fas used. 
Shoim in Pigs. 50a akd 50b are the predicted upper-surface longitudinal 
strains~Y^^ at~ stations x"0~ and x»076~in, respectively# as a function of time 
after initial impact. At each of these two illustrative locations there is 
very little difference between the predicted peak strains produced by impact- 
ing fragments of postulated diameters of 1.00# 1.25# and 1.75 in — because 
of the fact that for the given level of initial KE attack# the peak response 
of the beam near midspan is such that the smallest radius of oirvature of 
the beam exceeds that of each of these attacking fragments. Sowever# under 
more severe impact attack# it is evident that at and near midspan the beam 
could deform enough so that the beam's least radius of curvature could 
approach or attes^t to become less that that of the non-defozmable attacking 
fragment; in the latter instance the fragment could interfere with or inhibit 
the beam from incurring such large bending (and membrane) strains as would 
occur under the sane severe initial KE attack from a "sufficiently small 
diameter" fragment. 

For the CB-18 level of initial KE attack# one would need to postulate 
an idealized fragment diameter equal perhaps to the beam span I in order to 
influence the paak iispact-indxiced strain or deflection response of the beam. 

4. 3. 2. 2 Effects of Width-to-Thickness Ratio on the Response of a 
Beam of Fixed Span and Mass 

For this Case B study, the reference alumimaa containment beam has 

dimensions h»0.10 in# w «1.50 in# and 1^8.00 in and is subjected to the CB-18 

c 

initial kinetic energy conditions of rigid- fragment in^ct attack. Examined 

here is the effect upon the beam's predicted tremsient and peak responses of 

varying the ratio w /h of the containment beam width w to the thicJcness h 
c c 

while preserving the cross-sectional area so as to m 2 Lintain fixed the total 
mass of the beam. The following summarizes the geometric combinations 
studied# the reference calculation time step size (At)^# and the At actually 
used: 
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Case 


h(in) 

^,/h 

(1) * ^o”^ 




e- 

■ax * 




. 

(rad/sec) 

B1 

.50 

.30 

1.67 16 

.577 

B2 

1.00 

.15 

6.67 8 

.289 

B3 

1.50 

.10 

IS 5.33 

.233 

B4 

2.00 

.075 

26.67 4 

.233 

B5 

3.00 

.050 

60 2.67 

.233 


Tl^ Step Sire 
(At)- 

(0»8f.2 


(I) 

a 

(Usac) 
.27 
• 55 
.68 
.68 
.68 


At Osed 

.30 

.50 

.50 

.50 

.50 


-tat „ti, ^ «tatai-»t potatatal 

(or the smallest peak strain)? 

r« tta Ptatutatad C8-M ta„i Of tap«:t .ttaol, «a „ 

« 1.0 io. tta ptodictad pota tadtatf.0. .ttta. (.1.0 taU.d~.«tato. 
> tad tta prtalctad pota oppot-tarfta. loogittalota .trtao (y ) .t 
•ld.Pta taatio. ta. «-«. 1 . Pi,. 51. „ . TiSiSfta- 

«ta tat-ta. ttata tta talota Ptattao... 5 .. 0 .^ otatrltaUta to 2 

^.tafta. .tttao. ft 1, .Ota ttat tat oolP dta. tta -.pta-taPtao. pota 
ncrease as decreases (tonard narrower and thicker beans) but 
^ tta ttadio, taotaltatlo. to ttat Ptal tppoP-taPfpo. .ttta. tatin 
Itapta fPtatlOO Of tta tatal, fopptta. tt. p... 

It fp„ Pi,. 51. ttat tta Ctataitatat ta«. .tpootop, -tact 

l.ta. to tta tataltat Ptat .tPtao 1. tt.p tapt tta lta,c,t p.tio . /h _ 

which means a wide and thin beam. However «•»,« « ^ ^ 

However, the present CIVM-JET 4B calcula- 

^ta ™ tat -taid fop Uta. OP «»u t/-^ tac. 0.0 Pta .PP.PP.P. 0 ...,^ 

tta ctaputta 

cod. tapptaitat. ta.«» . toptalta .tat. of taitaita .tP,,.. tata, „ 

- ta Titatatato ,4d, foP .Itatic tatatacP .ta oota, foP .lo^tao-P^L" 

Of taifoPtap-iotata ttata.. tta defltatiot tata-icP dtatatta L,tat 

icantly from that for a "narrow >»>».- e siHM-i 

tor a narrow bean for aspect ratios l/^ less than about 

2. Hence, the Fig. SI predictions of upper-surface tv 1 .r n 

ptad Uttl. ttad fop - /t > 40 . 

c — 
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Aaoth«r "liaitatioa” on the wide-beaa reeolts aay be deef r»hrt am 
foUoMS. For the pnaeat 2-d siouletion of the sphere-bea ii^ect ptoblea, 
the ^ee llt e d fra^BMot of given diefMter may be vlaaaliaad moat properly aa 
a eircolar cylinder of width w^ matc h i ng the beam width w^, with a density 
changed appropriately to preserve the given total fragment mass. On the 
other hand, one may visualize the idealized fragment to consist of a circular 
disk of given material density, diameter, and fixed width w^ such that the 
total mass matches that of the reference fra^Mnt. However, there are no 
means by which the "aiall" tridth w^ of the latter fragment can be recognized 
by the 2-d CIVM-JET 4B code, although this latter fragment nore 

closely the physical sphere-beam is^tact situation %diich-ezhlbits decisive 
3-d features. Thus, one can not expect a 2-d structural response 
as CIVM-JET « to provide realistic predictions of threshold ctmtaiment for 
fragments which produce significant 3-d structural behavior of the 
structore. 

Shown also^ln Fi^* 5xa is the maxinnTm iq)per**surface strain at 
station x-1.50 in as a function of w^; the peak deflection at station x»0 
is shown in Fig. 51b as a function of w^/h. Both of these quantities are 
seen to be rather insensitive to w^ for this fixed initial-kinetic-energy 
impact attack. 

Hote that in this Case B study the cross-sectional area is kept constant 

varied. Thus, the extensional stiffness is the same for all cases 

but the bending stiffness increases r^idly as w^/h decreases. Thus, it may 

be of some interest to examine the time histories of the predicted upper- 

surface strain at the midspan iapact station x-0 and at a station a short 

distance from the impact station: x»1.50 in; these results are shown in 

Figs. 52a and 52b, respectively. It is seen that at station x-0, the 

strain is predicted to occur eulier for the narrower and deeper beams. A 

similar but less distinct trend is observed for station x-1.50 in; however, 

the maximum strain predicted at "rerote" station x-1.50 in is almost 

independent of w^/h *Aereas the peak strain at x-0 (in the iaanediate impact- 

*^^®cted region where the beam acquires a velocity increDent frcm each of a 

succession of fragment u^acts) depends significantly upon w /h. This 

c 
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responM iasensitivity to w^/h is ei^ected at stations |x| > 2 for the 

prMKtiiiiifjMi fixed klMtic energy inpart attack, in accordance with 

St« Venant*s Principle; station x*l«50 in qnalxfies in this sense since fron 
the equation * At (E/p) ^nff ^ astiaated to be approxiaately *060, 

• 099, .099, *099, and *099 in for Cases BI through BS, re^ectively. In view 

of these L values, only the nodes boandinq the center inpacted eleaent 
ox x 

receive iapact-ia^arted velocity increaents^ there are 7 nodal stations 
between the nearest node receiving iapact- induced velocity increa en ts and 
station x*1.50 in. 

It is perhaps useful at this point to remind the- reader of the necessity 

of eoploying a sufficiently saall tiae step sixe At for numerical predictions 

of nonlinear impact- induced transient st ructural response — such as the use 

of the dVM-JET 4B coiQ’uter program. In carrying out the predictions for 

Case B4 (w /h»26.67, w «2.0C in, and h»0.075 in), a At value of 0.75 micro- 
c c 

seconds was tried; this is soneidiat larger than the nominal guideline for 

selecting an appropriate At for this cospoter program: (At) <0.8 (2/ai^^) 

= 0.68 usee. Hence, At*. 75 ysee corresponds to 0.88 (2/u ^ ). For this case: 

max 

(1) an unrealistically large value of '*** predicted. at x»0 and (2) at 

and near the time of this (Y, , ) , there was a severe spatial variation of 

11 p a x 

upper-surface Yjj along the center impacted element, %diereas a slow spatial 
variation is-es^ected physically and is observed in well-behaved cco^uter 
runs. Thus, the use of a too large At was suspected. Accordingly, the At 
was reduced to a convenient value of 0.50 microsecond, and the calculation 
was repeated; this resulted in well-behaved predictions, and the peak at 
x=>0 was plaxisible as seen on Pig. 51. 

In view of the erratic restilts obtained for Case B4 when At»0.75 ysec 

• .88 (2/o) ), the extensive printed-out results for Case Bl were studied 

max 

carefully since At»0.30 ysec * .89 (2/a> ) was used. For the latter case 

pax 

no evidence of erratic behavior was found. However, based xspon the rather 

extensive set of calculations carried out in the present study, as wll as 

in prior studies, it is concluded that restricting At such that At £ 0.8 

(2/0) ) will lead to numericallv reliable transient large-deflection elastic- 

max 

plastic structural response predictiions of simple beam/ring structures when 
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2 T"'” 1 . =iv».^ «. ^ .« 

^ JT 30 T» d.«, ■ 

one Aoses At < 0.8 ( 2 /u ). 

tttBX 


of Bean Thickness on the Response of a Beam of Fixed 
Width aivd Span 

To be disctissed here are the results for Case c wherein the 
fragment has a fixed mass and fragment dimeter (d- 1.00 in) and a given 
initial impact velocity which matches th^ CB-18 sphere-beam impact experiment 
The span and the width of the aluminum containment beam are held fixed with 
^=8.00 in and w^=1.5C in. Varied in a sequence of calculations is the thick- 
ness h of the beam to assess this effect 'on the peak impact- induced strain 
in this 2-d beam structure; a time step size At of ^.5 microsecond was used 
for all calculations. When h is varied, the result is a change in (a) the 
beam mass per unit span, (b) the extensional stiffness, and (c) the bending 
stiffness of the beam, but these stiffnesses do not change in the same propor- 
tion. Table 10 indicates the dimensionless parameters which remain fixed and 
those which vary under these Case C conditions. 


If one confines attention to only peak-response results (irrespective 
of when the peak occurs) , the principal dimensionless variables reduce to 
three; S^Q-h/d, and 6 ^ 3 =h/w^. since d, I, and w^ are held constant, 

each of t-hese variables varies proportionally to h. Hence, for convenience, 
the peak response results: (a) at the midsurface and at the upper 

surface and (b) peak deflection w/f, both at the midspan impact station x =0 
are shown in Fig. 53 as a function of beam thickness h. For this selected 
limited range of beam thickness, the peak values of total strain, membrane 
strain, and midspan deflection all vary linearly with h. 

The result for h=0.10 in represents a 2-d simulation of the CB-18 sphere- 
beam impact result; the experimental CB-18 sphere-beam impact case involved 
an impact severity and accompanying significant 3 -d beam response (near the 
impact station) verging on material rupture. . For the 2-d simulation by the 
CIVM-JET 4B code (and t.he modeling employed) , the predicted peak strains must 
be less than those produced in the actual CB-18 3-d response case. However, 
for discussion purposes it will be assumed that the peak i^per-surface strain 
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predicted by this 2^ aodel r^resents the "contaioBent threshold” value that 
can ^ used to judge idiether or not a given 2*d beam under 2^ fraqaent attack 
can or can not contain the attacking fragment. Accordingly, that "critical” 

Y . value is selected as 16 per ceht. In turn this means, of course, that 


the thinner beam (b>. 075 in) would have xxxptxired but the beam (h«.15 in) 

%a>ttld not since its at m»0 did not exceed 12 percent. 

One can regard this set of Case C calculations as simulating numerically 
the procedure used in the NAPTC oontaiment tests wherein contaiiment rings 
of fixed inner-surf auire radius and fixed width were varied in thickness and 
then subjected to a given fragment attack (fixed number, mass, and initial 
kinetic energy) to determine eaq)erimental ly the thic)cness for threshold contain- 
ment, Hence, if desired, one may di^lay the present results in a fashion 
analogous to Fig, 42 for the HAPTC results; namely, displayed in Fig, 54 are 
the "present predictions” in terms of the ratio of beam mass to fragment mass 
versus the ratio of fragment initial kinetic energy (1489 in-lb) to 

the containment structure energy absorption index O , As xx>ted earlier 0 

times the volume of the beam material and 0 is 8700 (in-lb) /in*^. For conven- 

c 

ient reference, the quantities diaracter i zing the Case C calculations are: 


Case 

h 

h/d 


h/1 

V 

(in^) 

m 

c 

“f 

(KE)j 

0 

c 

Remarks 

Cl 

,075 

.075 

• 050 

• 0094 

0,90 

,6 

• 190 

Did not 
contain 

C2 

,10 

t 

H 

O 

• 067 

• 0125 

1,20 

,8 

• 143 

Containment 

Threshold 

C3 

,15 

• 15 

o 

H 

t 

• 019 

1,80 

1,2 

,095 

Contained 


Included for convenient reference in Fig, 54 are the data points corresponding 
to the Ref, 16 sphere-beam impact experiments; the symbols O and X mean, 
respectively, that the fragment was contained or not contained. Similarly, 
the symbols O and + mean that the predictions (Cases Cl, C2, and C3) indicate, 
respectively, containment or non-containment for the CB-18 level of Jdnetic 
energy attack. 
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It is evident that one could carry out CIVM-JET 4B calculations to 
determine the beam thickness h required for threshold containment of other 
values of single-fragment kinetic energy midspan attack such as, for example, 
the CB-9 or the CB-14 values of 686 or 1804 in- lb on aluminisa containment 
of 1.50- in width and 8. 00- in span. For these cases it %iould be 
postulated that the selected critical oi 16 per cent for threshold 
containiient applies. 

It should noted that the Case C ranges of beam dimensions are such 
that the structure should exhibit beamlike behavior (rather than platelike 
behavibr) as the coog>uter code employed assumes. 

4. 3. 2. 4 Effect of Beam Width on the Response of a Beam of Fixed 
Thickness and Span 

To be disoissed here are the results for Case D %rherein the impacting 

fragment has a fixed mass and diameter (d»1.00 in) and a given initial 

velocity which matches the CB-18 sphere-beam in^jact experiment. The span 

and the thickness of the alumina# containment beam are held fixed with 

t-8.00 in and h»0.100 in. Varied in a sequence of calculations is the 

%fidth, w , of the beam in order to assess this effect on the peak impact- 
c 

induced strain in this 2-d bean structure. Note than when w is changed, 
the monbrane stiffness and the bending stiffness are changed by the same 
percentage: also, a change in w^ results in a change in the beam 

mass per unit span. For these calculations a time step size At of 0. 5 micro- 
second was «nployed. Table 10 indicates the dimensionless parameters which 
roaain fixed and those which vary under the Case 0 conditions. 

If one confines attention to only peak-response results (irrespective 
of %*en the peak occiirs) , the principal dimensionless variables reduce to 

two: 6 .*w /w^ and S,-*»h/w ; since only 2-d behavior is dealt with in the 
11 c f 13 c 

calculations, S ,»h/w becomes the only meaningful dimensionless variable. 

13 c 

For convenience, therefore, the peak response results: (a) the midsurface 

and the upper-stirface values of (Y,,) and (b) the peak deflection w/i, 

XX DUUC 

both at the midspan iogxact station x*0 are shown in Fig. 55 as a function of 

containment beam width w . For this limited range of w values (from 1.12S in 

c c 

to 2.25 in), the peak value of midspan deflection varies linearly with w^. 
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but the membrane and outer surface peak strains show a concave upwards 

relationship as the wid^ decreases. This effect occurs because of the 

increase in curvature at peak deflection as w^ decreases, -and-wi-th-a— constant 

thickness, h , the bending strain will increase, 
c 

As w varies, the mass per unit span of the aluminum beam varies. Thus, 
c 

for a given fragment ii^act velocity, the velocity ia^arted to the impact- 
affected spanwise region of the beam will vary with w^ accordingly, with the 
result that the "rate of absorption of the fragment’s kinetic energy" as a 
function of will also vary significantly with w^ as sho*m in Fig. 56. 

Both tUs subsection and Subsection 4. 3. 2. 3 have dealt with the effects 
of varying the beam mass by varying" either the beam width w^ or the beam 
thickness, h , respectively,’ while holding the fragment properties constant. 
Figure 57 is a ce«?>osite of Figs. 53 and 55 showing the effects on peak mid- 
span deflections and peak outer-surface and membrane strains as the mass per 
unit span is varied (by varying either the width or the thickness) . At 
equivalent values of beam mass per unit span, the t>;o results show the effects 
of varying the beam cross-sectional area. Referring to Fig. 51 in which five 

values of w /h are plotted for a constant beam mass, the effects of varying 
c 

the cross sectional aspect ratio are illustrated and apply directly to 

the results of Fig. 57. In Fig. 57 there are two mass-variation plots, one 
for a variable thickness and the other for a variable width as tabulated in 
the following: 


Thickness Variation Width Variation 

w ■ constant » 1.5 in h = constant » 0.1 in 
c 


h w /h 

c 

w*h»A 
c c 

w 

c 

w /h 
c 

w-h-A 
c c 

0.075 

20 

.113 

1.125 

11.25 

.113 

0.100 

15 

.150 

1.500 

15.0 

.150 

0.150 

10 

.225 

2.250 

25.5 

.225 

tabulation 

it is 

noted that 

there 

are three 

values of beam cross 


sectional area, A ; at each value of A^ there are two values of w /h corres- 

c ^ . ^2 

ponding to a variation in thickness or a variation in width. For “ • 113 in 
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tlw — valu« of produced by varying only the width produces higher 

P^ fr strainl, which is expected fro* ex i mln i n g Fig. 57. For an increase in 

A to .225 in^ (or an increased waaa), the s««llest value-of-w-/h-i)s-attainad 

by a thickness variation. The peak strain is slightly larger than that 

attained by varying the width but the increase in strain is not as largo as 

be expected fron an examination of Fig. 51. An e^lantion of this is 

readily attainable by noting that the relative change in slope of 

vs. w /h as seen in Fig. 51 is dependent upon the magnitude of A^(A^ - w-h). 

C — 

It is likely that as A^ decreases. tl» change in slc^e will increase, while 

for an 1 n ^T Aaning A the effects of varying the aspect ratio will not be as 

c • V • 

great as that shown in Fig. 51. 

All of the parametric variations examined above were performed nneri- 
cally on containment beams «ihich were expected to demonstrate largely 2-d 
deformations. The principal result of the studies conducted in this subsec- 
tion and in Subsection 4. 3. 2. 3 is the demonstration of the effects of varia- 
tions both (a) per unit span and (b) aspect ratio, and their interde- 

pextdency. 

^ ^ t r««.Ants on Scaling Effects on Beam 2-D S tructural Re^onse 
Imt attention be restricted to a "containment beam" of appropriate 
structural proportions and/or pit^portions in relation to the impacting 
fragment(s) such that 2-d structural response will ensue. Further, let the 
coB^wnent of the velocity of the impacting fragment normal to the impacted 
surface be large enough to produce large-deflection elastic-plastic transient 
structural response (vp to,' perhaps, tensile npture) but small enough to 
avoid the -shear plug" type of failure. Within this 2-d structure and 2-d 
iapacting-fragnent framework, the CIVH-JET 4B couputer code can provide 
reasonable estimates of the transient structural re^onses (within its 
specified of validity) of single-layer 2-d con t a inmen t beams. 

lapact attack against a given target structure nay be characterized by 
(1) the time-zero location or the spanwise station and time of initial inpact 
of of the n attacking fragments as well as (2) the following data on 

each impacting fragment: (a) mass, (b) mass moment of inertia, (c) ccmponents 

of translational velocity, (d) rotational velocity, and (e) size of each 
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I 



■ I 

j 


id e allx ed aoadafoiaable fragnent) thaae fragmants are required in the present 
code not to interfere with each other. The resulting structural response for 
a given (postulated) attack can be computed readily, but the effect of -chang- 
ing one of the aany variables in this n-fragnent attack is not readily deter- 
mined in this highly n o n li ne a r problem eicc^t from direct calculations. 
However, if one confines attention to a aioole-fraoment 
attack at a givMi spanw ise station of the structure, one may assess more 
readily the effect upon the n a Kimum predict strain In the containment 
structure of moderate changes of (a) the cross-sectional of the 

2-d structure and/or (b) the kinetic energy (via mass and/or pre-i^>act 
velocity) and/or size of the Idealized attacking fragment. Hence, it is to 
this restricted set of conditions that the following observations are directed 
For the postulated single-fragment Iag»act attack (that is, a fragment of 
given diameter, mass, and perpendlcular-laqwct velocity) against, for example, 
a contaiment beam of given material, of given span, width, and thickness, 
and having both ends ideally clamped, the location and value of the 
extensional strain predicted will depend upon tte spanwise location of the 
initial fragment-beam contact. If d«<l, midspan impact will result in having 
*^11^ max midspan — and -nq>ture" at that station can be expected if 

the pre-impact (KE)^ is sufficienUy large. However, if initial fragment- 
beam contact occurs between midspan and a claaq>ed end,_the strain y 

will, in general, no longer occur at the spanwise station of in itial izq>act; 
as the initial-isq>act spanwise station sqiproaches one end from midspan, the 
location of (Y^j^)^ is expected to shift -rapidly- toward and to remain very 
close to the clamped end.* Accordingly, it is necessary to consider single- 
fragment iaqiact at only a fixed spanwise station in order to assess readily 
the effects of beam and fragment parameter changes upon the peak strain 
produced in the structure. For this purpose, therefore, let attention be 
confined to midspan perpendicular impact against a given beam by a single 
non-rotating idealized fragment. 


w 

It should be noted that the 43 equal-length finite element modeUng of the 
beam by elem^ts with 4 DOF/node is inadequate to represent accurately the 
strain behavior near the clamped end. A much finer finite-eleotent mesh is 
needed near a clai^>ed end to provide reasonable strain predictions in that 
region; this, however, will result in a much larger uhajj and an attendant 
required smaller At for the timewise solution by CIW-JET 4B. 
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With rwapwet to the fragment itself, its modeled behavior may be 
described in terms of 3 qnantitiess diam e t er , mass, and velocity. First, 
assumii^ that the contaisment ben has been-modeled-by-a-suffieient-nunber~of 
(uniform- length) finite elsnents so that "converged striictural response" can 
be realized both from (a) basic structural modeling requirements and (b) an 
adequate description of the l oc al impact-interaction "monehtum transfer" 
behavior, it has befm shorn earlier that fragment diameter has a "negligible 
effect" upon small hovever, for a large enough 

to produce large deflections and substantial changes of curvature; of the beam, 
a large-diameter fragment can constrain the structure aiad thereby ; change 
significantly both the spanvise distribution of and the value-and- 

location of On the other hand, if the idealized fragment diameter 

is very small, neither the spamrise distribution of peaJc nor 
will change since no fragment-beam constraint occurs and the basic finite- 
element modeling of the structure will leave unchanged the value and the 
spanwise distribution of locally-iaparted velocity changes esq)erienced by 
the beam from a fragment-beam impact — for sufficiently small diameter values 
of the idealized fragment. 

Hence, let the structure be modeled adequately and the idealized fragment 

diameter be properly small in the above sense. In that event, the remaining 

"fragment variables" (mass and velocity) can~be xepresented by only one 

significant (amd scalar) quantity: (KE)^^, the pre- impact Jcinetic energy of 

the fragment. Thus, one can evaluate (Y. ,) as a function of (KE)^ A3 
^ 11 max fo c 

for the postulated containment beam of fixed h, w , and i. 

c 

It would be interesting and instructive to see if the energy ratio . 

(KE). A3 for threshold containment (based upon producing a (y, .) of 16 
to c 11 max 

per cent) by a beam of fixed span i«8.00 in and width w^>1.50 in but changed 
thickness would be the same for vairious values of initial fragment kinetic 
energy would ctumge somewhat, since the only changes would be in 

(a) the beam thickness h and (b) the initial 2-d fragment kinetic energy or 

velocity (since a fragment of fixed mass is employed) . Accordingly, the * ^ 

1 

following calculations were carried out by modeling the beam with 43 equal- i 

i 

length 4 DOF/node elements and using an idealized fragment diameter of 1.00 in j 

-3 2 ‘ 

with a fixed mass of 0.315 x 10 (lb-sec )/in: 
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Case h(in)* 


C4 

CS = C2 
C6 
C7 
C8 


.046 

.100 

.121 

.200 

.200 


'“'fo 

‘“'fo^O 

Used in 43- 

(in-lb) 


Element Model 

686 

.143 

.186 

1489 

• 143 

.186 

1804 

.143 

.186 

2978 

.143 

.186 

2978 

.143 

.400 


Spanwise Length (in) of 


^eff 
Plausible 
Value - 2h 


.092 

.200 

.242 

.400 

.400 


rtote that case C5 (the sa«e as C2) represents the naainal C8-18 bea- 

and i„pact coalitions and also serves as the -2-d threshold ^ 
reference- corresponding to (y , x-0 of ^ containment 

For *-hio 4 ^IVmax * 16 per cent and /D «n 143 

For thi, .M tt. IMlcM ; 

686 (CB-9) to 2978 In ik / u* v • from 

- c„e, „ ^ - - ... 

«<u=h th, tM, th, Jin^ 

in thesA ^ regxon being used 

iTpatr^LTl^ir -e distance fl the 

^pact poxat to the midspan of the next element - which ic , 

0.186 in. on the other hand, as argued in Subsection 22 h 

pUuslbl. for this iswct..l,.ct«i 

IT. “ ™tr « - - — 

c... c%, ,r;.“ r "" " 

-.1 (.) UPPC sutlsc. str.1. « th. .ahsp.. 

energy loss ratio ((ke) - (ke) i/fKEi ^ ragment kinetic 

eh * ‘ represent, respectively, the initial and 

e current kxnetic energy of the fragment. Both of these time histories 
ar. v»v stsat„ Costs CS ... CS. hot those for Csses C4 sod C7 2 
distent, the «ilu. of •»> Its ti™ of occ«re»:e are ^^te sllIL 

_e~ee. Csses CS eod CS hut er e distinctly different for Cases C4 and C7. 

»cnl«a valees for h are used here for the sahe of co„y.„ie„ce. 
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The *ii^roper effective- length Dodelisg” for Cases C4 C7 be 

remedied by using finite elements of different and appropriate length for a 
spaiwise region centered about the midspan impact point. For Case C4 this 
would mean using elements of about- 0.092-in length in the central region and, 
for efficient, longer elements for the remainder of the beam. On the other 
hand, case C7 could be -remedied- without changing from the present 43 egual- 
length element model by distributing the iog>arted velocity increments to the 
4 rather than the 2 nodes nearest to the midspan i^act point; rather 
convenient and sisq>le diange has been used in carrying out Case C8 tdiere the 
‘eff is: 0.40 in compared with the physicaUy plausible estimate- of 
0.40 in (=2h) . The results from Cases C7 and C8 are coog>ared in Pig. 59 
(which corresponds to Fig. 58a). 

In Fig. 59 the peak strain of Case C8 is greater than for Case C7 (idilch 
are identical nuoerical examples except for a variation of the -iaqpact-affected 
length-) . it would be expected that as the ingwict-affected length is increased, 
the effects of the iag^ct interaction would be spread over a greater region 
and reduce the peak strain noted at the point of Impact. However, as the 
iagwct-affected length increases, the mass of the beam in the impact- 
affected region increases; for the CIVM (coUision imparted velocity method) 
scheme of subimpact momentum transfer (Appendix A) , the amount of kinetic 
energy transferred to the beam increases during each subir?)act. This may be 
seen more readily from the following tabulation and sketch which depicts for 
the present 43 equal-*element example: 

(a) impact occurring at the center of the element centered at beam 
midspdQp 

(b) the adjacent equal-length elements, 

(c) various values of assumed impact-affected length e^ressed 
in terms of elenent length, 

(d) the number of associated impact-affected nodes on one side of 
midspan according to itself and according to CIVM-JET 4B logic, 

(e) the total kinetic energy transferred to the beam in a single (first) 
impact ratioed to the pre-impact fragment kinetic energy (XE) for 

( 1 ) the entire impact-affected beam region and 

(2) only the center 2-element region associated with the two nodes 
on either side of the midspan impact point, and 

* 

See pages 17-19 and Fig. 10. 
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(f) the ratio of Itoa (2) to ites (1) expressed as per cent: 




No. of Affected. 
Nodes on One 
Side of Midspan . 

Ratio of KE transferred- 
to-Beam to (KE). for 
First Impact ® 

Per Cent to 
Region of 
Two Nearest 

liengtny 

^ef f 

By Code 

To Entire 
Beam 

To Two Nearest 
Nodes 

Nodes 

0 to 1 

1 

1 

.250 

.250 

100 

1 to 1.5 

2 

1 

.250 

.250 

100 

1.6 

2 

2 

.289 

.265 

91.6 

1.8 

2 

2 

.340 

.276 

81.2 

2.0 

2 

2 

.370 

.278 

75.0 

2.2 

3 

2 

.390 

.276 

70.8 

2.4 

3 

2 

.401 

.272 

67.8 

2.5 

3 

3 

.406 

.270 

66.6 

2.6 

3 

3 

.430 

.273 

63.6 

2.8 

3 

3 

.463 

.273 

59.0 

3.0 

3 

3 

.490 

.271 

55.4 

3.2 

4 

3 

.508 

.268 

52.8 


It is readily seen that the central portion of the beam surrounding the 
point of iB 5 >act received more kinetic energy for Case C8 during the first 
subimpact when the effective length is chosen as two times the beam thickness 
than for Case 7. This tends to increase the peak strain predicted at the 
impact location. 

A similar kinetic energy transfer calculation was carried out for Case C4 
with the present element length and for a reduced element length. The initial 

* 

Note that this value for Cases Cl through C7 is .532 and for Case C8 is 2.09. 






energy absorption ratio for the central two element region with the present 
element length (0*186 in) is 0*066 trtiile for the reduced element length 
(0*092 in) the ratio would be 0*034* This would result in a large 
in the rate of absorbed kinetic energy, which will lower the peak strain* 

The ratio of container width to thickness, was not preserved 

during this study* It is conceivable for scaling purposes to vary the mjia« 
per unit span by varying the beam width and holding the thickness constant* 

For a given mass or cross sectional area, A^,j the containment beam will have 
two different aspect ratios, w^/h, for a single mass value by the variation 
of either the thickness or the width' and holding the other two dimensions 
constant* Referring to Pigs* 51 and 57, it ^Id be eiqtected that Case C4 
would attain higher peak strains if the width had been reduced rather than the 
thickness* However, from Pig* 57, a small decease in peak strains is expected 
for Case C7 with an increased width* 

To this point in Subsection 4* 3.2*5, a change in (KE). has been assumed 

£0 

to occur by a change in the fragment velocity v; the fragment mass m^ has 

been assumed to be fixed* Alternatively, one could change (KE)^^ by holding 

V fixed and changing m^* Sv^pose, for exan 5 )le, that is increased by 

75 per cent over the CB-18 ve^ue by maintaining v»2794 in/sec but increasing 

by 75 per cent* On the basis of calculations already done for a contain* 

ment beam of fixed span (l»8.0C in) how can' one choose beam dinwng< r. n < ; w and 

c 

h such that at midspan will be a pre-selected value such as the 16 

per cent found for the basic CB-18 beam-fragment intact condition (Case C2 or 

C5)7 For this endeavor, let it be assumed that a matching of (KE)„ /U shall 

fo c 

be required — thus defining the cross-sectional area w^h as being 75 per 

cent greater than the CB-18 value of 0.15. This means, in turn, that the 

ratio of beam mass per unit length to fragment mass is exactly the same as 

in Cases C2 and C5 for the CB-18 example itself. Further, if one requires 

that the ratio w^ to be the same as in the CB-18 beam (i.e., 15), one obtains 

haO.132 in and w^*l*98 in* Hence, all of the dimensionless parameter values 

of the postulated example match exactly those of the CB-18 problan except for 

h/l and w^/£; these two parameter values are sufficiently close to those for 

CB-18 that little effect is expected. Accordingly, one expects (v ) for 

11 max 

the proposed beam-fragment iaq>act problem to match rather closely that of the 
CB-18 example (see Pig. 51a) . The predictions of t^per-surface 


as a 



function of at the midspan impact station x*0 for these two cases au:e 

shorn in Pig. 60, and confirms this e^^ctation. 

It spears, therefore, that the scaling of one conta inm ent configuration 
to another* in order to contain a geOTetrically similar fragment can not be 
accooplished by merely scaling the volume of the beam to preserve the ratio 
of fragment kinetic energy to beam energy absorption capabili^ . 

Insteeul the proper scaling or preservation of the ratios of beam mass to frag- 
ment mass and beam width to beam thickness must also be taken into account. 

4.4 Conments on Container Design 

The effects that various material, geometric, and k inem atic parameters 
have on a fragment-container impact-interaction have been discussed briefly 
in the preceding subsections of Section 4. Two sets of experimental data 
(NAPTC multi- fragment impact on steel rings and MIT-ASRL single-steel-sphere 
icpact on alimiinum beams) were anatlyzed in order to identify some of the 
primary parametric effects. In addition, a p l ausible 2— d numerical analysis 
was applied for both of these experimental sit\iations for a more extensive 
paurame ter- variation anaiysis. 

There is a basic difference between the numerical analysis caurried out 
on the NAPTC fragment-ring ispact data (Section 3 and Subsection 4.2) and the 
MIT-ASRL fragment-beam irpact data (Subsection_4.3) . The NAPTC fragment's 
initial impact amgle is very different; all of the fragment-ring inpacts 
occurred at oblique angles %diile all of the fragment-beam impacts occurred 
at normal OO"*) incidence angles. In order to make certadn conparisons 
between these two parametric studies, it is necessary to ex a m i n e the effects 
that ispact incidence angle has on the impact interaction. 

Figures 61a and 61b show the results of a fragment-beam irpact inter- 
action for a 60® initial impact incidence angle. The fragment and beam 
are geometrically identical to the CB-18 fragment-beam inpact experiment 
described in Section 2 and in Subsection 4.3. The fragment is a rigid disk 
with a 1.0- in diameter. The double claop^ 6061-T651 aluminum beam has a 
length of 8.0 in, a width of 1.5 in, and a thickness of 0.1 in. The CIVM-JET 4B 
CQoputer code is used for the analysis %#ith 43 equal- length finite elements 
to represent the beam and an EL-SH material approximation for the aluminum 
* 

Based on a structxire of fixed span as stated on page 105. 
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via the mechanical sublayer model and stress-strain Fit A (see page 8) • The 
fragment was allowed to ia^ct the beam ait the midspan location %rith a vertical 
velocity of 2794 in/sec and a tangential velocity of 1613.12 in/sec (total 
velocity of 3226.23 in/sec, a total Icinetic energy of 1985 in-lbs, and an 
Initial inc idence angle of 60* measured frccn the tangent to the beam lower 
sxirface) . The coefficient of friction, Ur was set equal to zero (as was the 
case %#ith all calculations reported in Subsection 4.3). Therefore, this 
fragment's translation^a kinetic energy caaponent normal to the beam (1486 
in- lb) is identical to the total translational, kinetic energy used in the 
earlier normal ixopact calcxilations. 

Figure 61a is a time history of fragment trajectory angle 9. One set of 
data points corresponds to the present exampie of a 60* initial incidence 
angle. The second set of data points corresponds to the CB-18 normal inpact 
incidence case. Note that both examples have identical pre impact normal 
velocities of 2794 in/sec and that fragment rebound occurs almost at the 
same time ^ter initiad. i^act (approx. TAII»650 microseconds) . For the 60* 
initial incidence angle the fragment is 0.97 in from the beam midspan location 
at the time of rebound. 

Figure 61b is an ipper s\irface strain profile at TAIXa655 microseconds 
for both the 60* and 90* initiad incidence angle cases. This corresponds 
to the time of peak transient response for the 60* caise, is about 25 micro- 
seconds before peak response for the 90* case, and is the approximate 
of fragment rebound for both cases. Note that the 60* case has a much wider 
region of significant beam deformation, but the peak strain is greater for 
the 90* case. 

The absolute peak strcuLn predicted for the 90* case is 0.16 i^ch is 
the assumed threshold failure strain. In the 60* case the beam permamently 
absorbed a total of 1660 in- lb of energy while the beam in the 90* case 
absorbed permanently only 1411 in- lb of energy. It is apparent from this 
study that the normal pre- inpact velocity conponent is the most critical 
(for an interaction in which a negligibly small vadue of friction coefficient, 

U, is used), and that an anadysis aissuming a normal inpact for the same 
initial total kinetic energy will be conservative . Therefore, parametric 
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that vere identified froa the analytical 
tudy of fragment-beam noanal impacts should also be appUcable 
ring oblique impacts. ^ fragment- 

It should be noted that there' are some basic. • . 

studi., of Sutaection, 4.2 4.3 (or ^ Wr««tiic 

i.cor.«foo4. „2 Of r 

or wldf. .uodod fo . foooffo. Of oUcZ.r::f::: (“tT 
MditxooaUy, oon. of U» «u41m coDsld«i«l tb. poMibiUw of f 
d.flectioo (. 0 . soction 5) r4tl»r th«, tom conttiomt 

tion) both assume that a soeelfic , geometry or fragment deflec- 

specxfxc tone on the container will be imacd-m, 

::r: r.rjrrr.r: “ 

ro ft ™ot be 444^ tbet -oerr^orrftrLs^i:^ " 

protected, consideration concemina what mi thf k 
probebUlt, b.p.tt-t.,fo. or tbe ptot.otfo. of oof, speclfio 1 ‘ 

-ru my .fdoly besed o. fre,o..t type, cooteloer location 

■bd my other ptoblo. departent paraoetets Tb, aoal ■ ' 

«0 4.3 14 «4.t to be -,eo„,i. .44 applTc^le L“ IrTb"” 
rib, type cobtaloers. »ot cbeldered. bcmer, m " 

iramot d.flectlo.,or .elective fr.,„e„t peaemtL 

Bavlb, ldeotlfl«i tbe mt problneot para«ter 4 affectlb, tbe reek 
sponse of a general ring-like or beam-like container t ■ ' 

tions, and restricting oneself to small strain theory) ' ! 

^Tl *»-^niP.Vb-.rlcal, be appUed^to'^otT.r'^^^’ 

gn irst, two basic criteria t^Kjn which the desion of an 
container might be based should be identified: 

(1) A s^cified region needs to be protected. A fixed container 

surface area is defined and for this geometry a minimum weight 
container needs to be designed. 
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( 2 ) 


A "costing function" is applied to the air vehicle system and 
a fixed mass is allocated for containment purposes. A container 
and geometiy must then be decided t^n to maximi ze the 
containment potential, of this fixed-mass structure. 

A design program will most likely utilize both of these criteria during 
different stages of the design process. First, a fragment generation problem 
is identified and a specif ied' region to be protected is then defined. An 
analysis will then proceed to determine a first estimate of cost (weight, 
labor) to n various sizes and configurations of fragments for that 

particular aircraft engine. ? Then, a "costing function" based on the cost 
of protection versus the level of uoncontainment that can be tolerated (or 
regulated) is defined and an ;upper limit is placed on the mass of the container 
and on the costs of manufacture, installation, and maintenance. 

jn each phase of this design work, a set of empirical and/or analytical 
impact data is needed. A design procedure which depends totally on en?>irical 
data of adequate scope will arrive at a design which might be viewed as having 
been "proof tested", but likely will be costly and time consuming. If there 
is an attempt to decrease the experimental costs by ^lrbitrarily or "intuitively" 
limiting the scope of the experimental program or by using "scaled" experi- 
mental models, there is a probability of alternate containment designs being 
undiscovered or of the erroneous assessment of key parametric effects (because 
of Lnsufficient data or inadequately scaled experiments). Conversely, a 
design program relying solely upon ^ analytical/numerical prediction data 
should be less ejqjensive but may not address the actual containment problem 
adequately because of the basic restrictive assumptions of the analytical 
method, and will not provide a "proof tested" design. 

A more rational and cost— effectxve design process might consist of a 
otnaiT preliminary set of "full-scale" experiments to assess the effects of a 
few basic parameters and to obtain sufficient d at a to calibrate a numerical 
analysis to "compensate" for tl*e "restrictive assumptions" in the analysis. 

The analytical technique (once correlated with preliminary data f.om the 
pilot tests) could then be used for more extensive parametric studies. Also, 
it could be used to aid in the design of an "adequately- scaled" model of the 
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to 0. o«a i. .s.rl-rt.1 

— .Oeqoacy of O. eontaiit~t-«™ctor. th« a.»i9-« 

verified by a saall selective' set of full-scale proof tests. . . , 

tte paraisetric effects that are identified in Subsections 4.2 and 4.3 
could be as a guide for detetaining the original pilot experiaents 

and the series of analytical studies to be used in the design process. 


SECTION 5 


ON THE USE OP A DEFLECTOR FOR FRAGMENT CONTROL 
5. 1 Problen Outline 

TtM preceding sections have dealt with the cootainBent of a fragment, ^ 
%ihich the impacted structure is required to absorb essentially all of the ; 
fragment's translational Icinetic energy. For many practical aM>Ucations, 
only a particular region is required to be free from fragment intrusion, and 
a structure could be utiUred to deflect the fragment from a path that would 
otherwise enter this region. A deflector can be described by the same set; 
of dimensionless parameters that has been used already to characterize a 
container, plus three additional interactive parameters; a, the angle change 
of the path of the fragment; 9 . the angle of initial impact; and the point 
along the structure at id»ich initial ijapact occurs. Table 10 Usts the 18 
dimensionless parameters used to describe a beam containment of a rigid 
spheroidal fragment. To this list will be added o (measured positive 

clockwise) , ® (measured positive counter clockwise) , and 621 -^ (fractional 

along beam). An illustrative cantilever-type deflector is depicted 
in Fig. - 62 with the appropriate approach and deflection angles (9 and a) 
indicated. 

For the deflector system, there are three equally-important dimensionless 
parameters which describe the effectiveness of the deflector; e^, the 
maximum strain attained by the deflector; z/h, ratio of maximum deflec- 
tion (z or w) to the beam thickness h which indicates the extent to which 
the deflector impinges upo:; the space that is being protected; and a the 
change in flight path of the fragment. There are. therefore, three means by 
which a deflector can -fail"; ( 1 ) the deflector remains structurally intact 
but the fragment is still allo%#ed to enter the protected region (insufficient 
path-angle-change a) ; ( 2 ) the deflector structurally fails and it either 

proceeds into or allows the fragment to enter the -protected region (e^ is 
too large); and/or ( 3 ) the deflector remains structurally intact and the 
magnitude of a is sufficiently large to deflect the fragment but the transient 
deflection of the deflector carries it into the region to be protected (z/h is 
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too large), in each of these cases the protected region suffers an intrusion 
by either the fragment or the deflector. 

The primary purpose of this brief section is to study the response of 
an illustrative deflector type structure to the variation of two of the 
diaensionless parameters discussed above: 9, and 5; the results will be 

assessed on the basis of the transient values of a and */h, and of e . 

There is an insufficient amount of experimental data on fragment-deflIjL,r 

and response to constmet -effects plots-, but based on the evaluation 
set forth in Sections 2 and 4 of this report, a brief n«Berlcal parametric 

st^rfy was performed and will be des«ibed. The remaining dimensionless 

parameters. 8^ - 8^g, wiU not be addressed in this section. These parameters 
have been discussed in Section 4 and it is believed that their roles will be 
similar for a deflector as they are for a container structure. 


— MuPeJ^^cal Model of a Cantilevered-Beam Deflector 
In order to keep the parametric evaluation of a deflector as simple as 
possible, a nmnerical model nearly identical to that used to model the -beam- 
containers of section 4 is utilized. The deflector is modeled simply as a 
cantilever beam of length 8.0 in, width 1.5 in. and thickness 0.1 in composed 
of 6061-T651 aluminum. The CIVM-JET 4B computer code is used to analyze the 
beam as 43 equal-length finite elements (44 lumped mass nodal points) . The 
material properties of the 6061-T6S1 aluminum are assumed to be elastic-strain 
hardeiung (EL-SH) approximated by the mechanical sublayer model using Fit A 
(see page 8) . A single attacking fragment is modeled as a 1.0-in diameter 
rigid cir^ar disk with a translational velocity of 2794 in/sec and a mass of 

.3815x10 (lb-sec )/in. This model is identical to the model CB-18 conditions 
used in Section 4* 


The boundary conditions of the beam deflector represent a crxtxcax tactor 
influencing the effectiveness of the deflector, m Subsection 4.4 a numerical 
study of an identical beam that was doubly clamped was performed for a frag- 
ment incidence angle of 6-60- . A fragment deflection occurred during that 
analysis; however, for the purposes of ms study, a deflector will be 

defined to have one end free. Therefore, the analyzed deflector has one end 
clamped and the other end free. 
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5,3 Pa raaetr ic Evaluation of a Cantilayer -Beam Deflector 

In analysis# parameters (Table 10) are held constant and 

only parameters 6^^ and varied. Thus, the evaluation of 

a fragment-deflector impact is based on an examination of the foUowing: 


(a) 

®18- 

z/h • fj^(6,5) 

(b) 

®18- 


(c) 

®19“ 

a ■ 


Deflector Response 
Fragment Response 


There are an i^inite number of values that could be chosen for either 
5 or e. however, for this analysis, only a few selected values will be 
examined. If the b4am %«re freely su 4 >ported (resting on wires, for instance) , 
the best possible impact location for maximum energy absorption would be at 
the centcur. Because the beam is clamped at one end and since the beam 

have been chosen to be equal to the -threshold containment size" 
for a doubly-clamped beam vhen injected at midspan Iv the same fragment size 
and KE as utilized here) , the probabiUty of a successful deflection/contain- 
ment should increase as the value of C decreases. Note that the CIVM-JET 4B 
analysis can not properly evaluate the 3-d stress state that will be found 
at the clamped end so that impacts vdiich occur close to the clamp wUl provide 
meaningless predictions near the clasq>ed end. Therefore, 5 will be selected 
as 0.5 < C < 1-0. For 5*0.5 the impact- interaction will initially be identical 
to the containment studies of Section 4, and this value (5*0.5) is chosen to 
be the first point of evaluation; the subsequent structural response will be 
different since now one end of the beam is free while the other is clamped. 

For of the 5 values examined, there are three values of 9 that are 

analyzed; 9-30*, 60*, and 90*; (where 0»9O* is an initially normal impact). 

The values e-90* was chosen to determine the amount of "tangential velocity" 
that the deflector can impart to the fragment as the beam deforms; also, "cross- 
correlations" can be drawn to the normal- in^iact containment sttidies of 
Section 4 for beams wit); both ends claaped. Corresponding studies of 9-120* 
and 150" were not performed because the deflected fragment would travel 
towards the claa^jed end in these two situations, rather than towards the free 


end emd# hence , would not "escape'*. 

^Assuming that transverse shear- induced failure can be ignored. 
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1 


Figures 63a through 63c shm the results of. a fragment la^act upon a 
cantilever beam for and 6«30*, 60*, and 90*. Figure 63a is a plot of a 

resulting from these three is^act configurations. Note that for all three 
situations, including the relatively; "shallow* iaq^act angle 6*30*, the fragment 
has become "embedded" irithin the confines of the deflector. Figure 63b depicts 
this situa t i on more clearly as a plo^ of beam profile and fragment location at 
905 microseconds after initial impact. Note from Fig. 63a that this corre- 
sponds to the time at which the fragments are all being deflected towards 
the clamped end since a is decreasing (the 8*30* case does not show a negative 
value of a during the tiM span presented) . The beam has deformed greatly 
and for each 6 case there is a portion of the beam whi^ is below any portion 
of the disk-like fragment. The impact-interaction with a doubly-clamped 
containment beam (Section 4) is generally finished by 900 microseconds after 
initial contact. It is interesting to examine the energy remaining in the 
fragment at this time after impact with a deflector : 



INITIAL 

IMPACT LOCATJON: 

5*0.5 

1 

Inpact 

Fragment Kinetic Energy 

Fragment Kinetic Energy 

Incidence 



■ 


Angle 

Before In^ct 

At .905 Microseconds TAII | 

e 

Total (KE)^ 

Normal (KE)^ 

Total (KE)g 

(KE)^ Parallel 
to^z axis 

(de?) 


(in-lb) 

(in-lb) 

(in-lb) 

30* 

1489.1 

372.3 

912. 8 [61]* 

47.4(13]* 

60* 

1489.1 

1116.8 

459.9(311 

296.6(27] 

90* 

1489.1 

1489.1 

544. 6 [37] 

534.8(36] 



The fragments for ea^ value of 6 possess the same pre-ia^act total 
kinetic energy? after 905 microseconds, the fragments still possess a signifi- 
cant level of kinetic energy, and will continue to impact the deflector. The 
peak transient strains, however, have already occurred by this time and 
Fig. 63c sho%rs a plot of i^per- surface strain versus deflector location at 

* 

The bracketed numbers are the percent ratios of the energies at 905 micro- 
seconds to the corresponding pre- impact energy values. 
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905 microseconds after la^wict, plus the value of the peaJc strain ^rtiich 
occurred during the intact interaction. Mote that the peak strains 

for 0-6O* and 90* occur at or very near the point of initial inpact. For the 
0»9O* case the fragment remains at the same x location (Fig. 63b) and there 
is only on^ significant peak for the upper surface strains. The 0-60- case 
shows the fragment traversing the span of the deflector* creating two 

of significant peak upper-surface strains. The first (and largest 
peak) occurs at the point of initial inpact whUe the second peak occurs at the 
spanwise location where the fragment has become "eirt)edded" (Figs. 63a and 
63b). However, the 0»3O« case shows a very small i^ak at, the point of initial 
inpact and a significantly larger peak »pper-surfaee strain at the spanwise 
location of fragment "embedding". The 0«3O" fragment still possesses 60 per- 
cent of its initial kinetic energy at 905 microseconds after initial inpact,- 
as can be seen from Fig. 63b, the fragment will continue to impact the deflec- 
tor at this new spanwise location creating a large peak ipper-surface strain. 
The fragments for 0-6O* and 90® have expended most of their original kinetic 

energy during the early impacts while for 0*30® the fragment still possesses 
most of its original kinetic energy and will expend it during a later succes- 
sion of impacts. 

-Further examination of Fig. 63c shows a predicted large concentrated 
strain at the clamped end for all three 0 cases. However, as was pointed 
out earlier in this subsection and also in Section 2, the 2-d CIVM-JET 4B 
analysis can not accurately evaluate the large 3-d strains which occur near 
a clasped boundary? also, even within the 2-d framework, a finer finite 
elosent mesh would be needed near the clasped end to permit a reasonable 
approximation to the spatially rapidly-varying strains near the clamped end. 
Therefore, care must be exercised when e xamin i ng the strain distribution 
along the deflector, and strain predictions at and near the clasped end 
should be viewed skeptically. 

Overall, the peak strains for these three impact conditions (0*30®, 60®, 
and 90® “ all at 5*0.5) are significantly lower than the peak strain pre- 
dicted for 0=60® and 0=90® for eux identical doubly-clasped containment beam. 
However, the beam deflections are very large and could result in either (1) 
the excursion of the deflector into the protected region or (2) the eventual 
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failure of the deflector at the clan^d end because of the large Ifi dnc e d 
strain and/or bending mooent. It appears that the fragments will continue 
to remain in contact with the deflector and will probably not enter the 
protected region without a tensile failure of the deflector. 

TO complete this brief examination of fragment-deflector impact.interaction. 
it is desired to consider a value of iniUal-impact-station locaUon K that 
wiU allow the fragments from all three 8 cases to escape from contact with 
the deflector. In order to -guarantee- this occurrence, a value ^-.875 was 
chosen for 8-30-, 60-, and 90- with all other parameters identical to the 
C-.50 situation. Figures 64a through 64c show the results of thi^ in 5 »act 
scenario. 

Figure 64a shows that Uttle fragment deflection (o) has occurred for 
any of the impact 0's. The 0-30- and 60* cases show identical final values 
for o; in all three 0's, the fragment has ceased impacting the deflector by 
300 microseconds after initial contact. Figure 64b shows the large (continu- 
ing) excursions of the deflector's free end; for 0-3O« and 60«, the fragment 
has completely left the vicinity of the deflector by 905 microseconds after 

initial impact. Examining the fragment kinetic energies, the following is 
noted: 



For each 0 case, a significant portion of the fragment's kinetic energy still 
remains at 90S microseconds after initial impact and it is evident from 


See footnote on page 114* 
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64b all three cases that the Cra 9 sients either have escaped or will escape. 
Figure 64c depicts the «^per-surface profile at this tijne with only one peak 
registered for each value of 0 and each peak i«)per-surface strain occurring 
at or near the point of initial impact. Note that very low strain levels nave 
been recorded at the clamped boundary because the rotations at the clai^ at 
this time (Fig. 64b) are also very small. It appears for all three values of 
9 that the deflector hM sustained a "permanent hinge" near the midspan location 
(Fig. 64b) and indeed this is a point of peak v 5 >per-surface con?)ressive strain 

(Fig.r 64c) . , , 

Overall^ an impact in the vicinity of C“0.875 appears to result in the 
"failure" of the deflector. The deflector has failed tiy allowing the frag- 
ment to escape with most of their original kinetic energies along paths very 

i 

close to their original path. The deflector also sustains extremely large 
deflections of its free tip, the excursions of which may easily enter the 
"protected region" . It appears that tensile failures of the deflector are 
not imminent or likely for large values of The deflector has formed a 
secondary plastic hinge at a location remote from the clamped end, resulting 
in relatively low bending strains . The peak strain at the point of inpact 
for the 0-90® case is relatively large (0.14 conpared with a "failure strain" 
of 0.16 deduced from an identical beam that was doubly clamped and subjected 
to the CB-18 impact conditions); however, since the fragment has escaped, a 
tensile fulure at this location would be inconsequential to the flight path 
of the fragment. It is interesting to note that the peak upper-surface 
strains for 0-60® and 90® at 5-0.875 are larger than those recorded for 
5-0.50. This is because of the larger levels of local bending that occurs 
for 5-0.875, which can be seen from a cooparison of Figs. 63c and 64c. For 
the 0-30® case the relatively low level of total energy t? ansferred from the 
fra^ent to the deflector for 5=0* 875 prevents upper-surface strains from 
becming larger than those predicted for 5=‘0*50. For 5=0*50 the 0-30® case 
had its peaJc \pper-surface strain located at a spanwise location removed 
from the point of initial impact. For 5=0.875 the fragment has escaped, 
which precludes the transfer of kinetic energy to the deflector by "later" 
subimpacts . 
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5»4 Cocnents on Deflector Design 

A brief overview of deflector structures that are assumed to be 
geometrically identical (except for the boundary conditions) to the previously- 
analyzed container beams was presented in Subsections 5.1 - 5.3. A simple 
parametric evaluation was performed on the effects of (1) the spanwlse 
location of i n i t ial impact and (2) the initial fragment incidence angle. 

The deflectors analyzed were assumed to have all geometric and 
parameters independent of spanwise location. 

The user of: a deflector type structure assumes that only a specific 
region is to be protected, thus the criticalness of a given fragment attack 
against a given <tef lector structure of this type is dependent \ 5 >on the span- 
wise location of toe initial Impact. It is conceivable that a systematic 
variation of toe Spanwise mass distribution of the deflector (in the region 
near the clan^, between the clamp and the initial Impact region, in the 
impact region, ud between the in^ct region and toe deflector's free end) 
will lead to a container of minimum total mass and an optimal mass distribu- 
tion for fragment-deflection purposes. The question of mass-distribution 
effects was not studied in this report, but could be carried out readily 
with the CIVM-JET 4B code. 

The boundeiry conditions for the deflector ue themselves extremely 
critical in determining the effectiveness of a fragment deflector. An 
actual deflector may be attached at one end to a "rigid" structure but may 
^U.so be attached at a number of spanwise stations to a relatively flexible 
surrounding structure (e.g., the thin engine casing, stiffeners, etc.) which 
will tend to limit the deflector's deformation and also absorb some of the 
fragment's kinetic energy. Considerations of this type are highly problem 
dependent and would require an extensive set of numerical or en^irical 
studies in order to assess the parametric effects adequately. 

The purpose of Section 5 was not to identify all of the parameters 
affecting a fragment-deflector impact-interaction nor to arrive at a deflector 
design by optimizing the mass distribution or the applied boundary attach- 
ments. Instead, only a brief illustrative examination of a structure that 
is identically geometric^aly configured for either a containment or a 
deflection role is intended. Recognizing the limitations of the previously- 
discussed deflector analysis, it is still possible to make some overall 
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cGments cmcernixkg the use of a deflector structure, for fragment control. 

In Subsection 4.4 a possible process vas reviewed for the design of a 
fragment containment structure using a combination of eapixLcal and analytical/ 
nisnerical data. A similar process is recommended for the design of a deflec- 
tor. A design is guided by cost (weight, manufacture, ins tallation, and 
maintenance) considerations. A careful decision must be 'made concerning the 
extent of the fragment control (the maximum size and energy of the considered 
fragment, the size of the region to be protected, and the implications of 
selective fragment penetration or deflection) based on these cost constraints. 
The final design is bounded by these considerations, and the allowable mass 
must be optimized to control the widest variety of possible fragment-structure 
iapact interactions. A deflector may be the "first choice" as a fragment 
control structure. However, if a container^ is designed, there is the possibil- 
ity that a portion of the container support structure might fail (either 
early in the intact- interaction or due to a direct impact by a. second frag- 
ment) • This situation would necessitate the treatment of the combined 
container and its failed support as a single "deflector type" structure. 
Therefore, a deflection analysis may be needed in a wide variety of fragment 
control designs. 

Given an adequate set of pilot deflection experiments for correlation 
purposes, the 2-d impact- interaction analysis described in this report should 
give a reasonable evaluation for a deflector design. The last stage of a 
design process for a contaiiner, deflector, or dual container/deflector would 
require full-scale proof- tests because of the inherent limitations of the 
present numerically-generated code predictions or of enpirical "scale-model" 
generated experimental data. 


That is, a partial ring or a panel which is intended to contain rather than 
to divert or deflect the fragment. 
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SECTION 6 


exploratori modifications and extensions of the ahalvsis 

6.1 General Conalderafcion«» 

in the present calculations discussed thus far, the fiaite-eleoeht 
analysis procedure described in Refs. 4 and 17 for the large-deflection 
elastic-plastic transient responses of BemoulU-Euler beaas and rings has 
been utxlired. The retention of the classical Bernbulli-Buler assu» 5 ,tion 
^t normals to the reference axis of the beam undergo no change in length 
during deformation implies a restriction that strains are not large, since 
large straxning will produce thinning (in tension) or thickening (in ccmpres- 
sron) of the body-s depthwise dimension. Also, the change-of-curvature part 
of the straxn-displacement relations used in Refs. 4, 17, and 20 (and named 
stram-displacement relations A. B. c. D. and E in Ref. 20; also, see Eg. B 5 
of Appendix B) are restricted to small membrane strains, since the membrane 
straxns have to be included in the Lagrangian bending strain expression when 
they are finite. The term -finite strain" (as compared with "small- or 
"infinitesimal- strain, refers to strains that are not negligible when conferee 
With unity, considering a curved beam, for example, as depicted in Fig. A.l 
with curvilinear coordinates C, >1. C in directions 1, 2, and 3, respectively 
xf the strain measure is the Green ("Lagrangian", strain tensor in nixed 
component form* e;q,ressions of the type (I+Iy") , for example, oedur in the 
constitutive equations (see Eq. B.21 of Appendix B) ; if Green-strain- 

tensor mixed component ^^e reference axis of the beam (or "membrane strain", 
expressions of the type (I+ 2 Y 2 , ?°-dependent part of the strain- 

displacement relations (see Eq. B.2 of Appendix B) . If so that l+Iy^-i 

one usually refers to the strain as being "infinitesimal" or "small"; otheLse 
the strain is^referred to as being -finite" or "large", since it is always 
true that Y 2 tr 2 ' ^he small-strain approximation breaks down first in the 
constitutive equations and later in the strain-displacement equations. 

T^or mix^ cemponents (identified by both subscripts and superscriots) ar. 
used here for general applicabiUty since in curvilinear co^JSS^S 
ormulation takes a much simpler. form when formulated in terms of the te 
mrxed co mponents (for example, yj instead of y^., and instead ofsij“ “ 

in rectangular Cartesian coordinates, is satisfied identically. 
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In addition to the BernouUi^EuIcu: assumptions, the mechanical-sublayer 

model of the material was used in Refs* 4 and 17; this model, in tarn, requires 

that the stress-strain curve being modeled must be monotonically increasing 

— the stress associated with this stress-strain curve must not decrease with 

increasing strain — and unloading must proceed elastically at the sam^ 

modulus as the original elastic modulus* Thus, since the stress «■ — a 

K p T 

(see Ref* 21 page 23, where it was assumed that the original mass density 

p was equal to the final mass density p so that t. was equal to the **true 
o . ^ : 

stress** a^) versus axial Lagrangian strain tensor component approximately 

exhibits this type of monotonic behavior (whereas vs* Yjj^ does, not) , a 

mechanical-sublayer fit to this uniaxial tension stress-stradn curve was made 

and the res\ilting values as reported in Subsection 2*2 (termed Fit A) amd in 

Subsection 2*4 (termed Fit B) were employed* The resulting stresses were 

used in the JET 3 and in the CIVM-JET 4B program as playing the role of the 

second Piola-Kirchhoff stress which appears in the basic governing equations 

[4, 17 and Appendix A of the present report] upon which the above programs 

are based* In view of the above considerations as well as the data scatter 

in the experimental ^measurements of uniaxial stress-strain behavior of a 

number of 6061-T651 aluminum test specimens, the adopted compromise procedure 

was believed likely to provide reasonable predictions of structural response 

involving small levels of strain* At what strain levels these approximations 

lead to unreliable predictions remains quite uncertain* Therefore, some 

additional exploratory studies of limited scope designed to remedy, in laucge 

part, these noted deficiencies were carried out and are discussed in this 

section* 

The modifications made consist essentially of five aspects as described 
concisely in the following, and discussed further in Appendix B: 

1. The proper second Piola-Kirchhoff stress tensor conponents 

appearing in the basic governing equations are employed in the 
computer program by making proper transformations in the constitu- 
tive equations of certain stress and strain measures ezrployed in a 

^In Ref. 21 rectangular CartesiM coordinates are used; for these coordinates, 

^ ^ However, in this section and Appendix B the “axial direction 

for a beam or the circumferential direction** for a ring is regarded as the 
n or •*2*' direction* 
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new mechanical sublayer fitting of uniaxial 

6061-T651 aluminum [21] as cited under item 2. ““ 

stress-strain datk f21J of the 6061-T651 aluminum 
material e^ressed in terms of the Kirchhoff stress t f • 

T for uniaxial) versus £* has been fitted ' ^ 

fashion by means of the me^^-^ " Pxecewise-linear 

ans Of the mechanical-sublayer model, where 


^ s Arr =— i- - ^ /j r \ 

e^-kf = U(uEj 


(6ela) 


(6.1b) 


aiid where 


A^(A) 


Z^(Z) 

Po(P) 


£* 


- lo«i WUed to th, t,,, 

reformed) specimen 

- th. ori,ital 

- POT orioiiai oodofootM 

volume of the material 

Z~Z 

- ^ - lon^ltuan .1 oloopoUo. .troio 

= f = Cauchy (-^e") stress = i ^ , 

i to f - icoritht^c J 

- Po 


P = Kirchhoff stress = a (i+g ). a = JL 

“ E A 

^o^o^o ~ ~ mass conservation 


[P/(l+E )]/A = a^(l+E )? P foj^ce)/ (undeformed area) == 
u o t q' • 
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It turns out [46, p. 29] that the uniaxial stress-strain curve for 
ductile netals such as aluminum is essentially perfectly antisymmet- 
for tension con^ression) and monotonically increasing 
when expressed in terms of and e*. Hence, the mechanical-sublayer 
model is well suited to describe this type of behavior. In particular, 
the u niaxi al stress-stradn behavior for the 6061-T651 material of 
Ref; 16 (as given in Pig. 29a of Ref. 21 for longitudinal tensile 
specimens 6f circular cross-section and converted to T vs. e*) is 
given by ^e following T^,e* pairs for mechanical-sublayer piecewise- 
linear fi€^ng: » (0,0); (44,200 psi, 0.00442); (49,200 psi, 

0.075); aiid (76,400 psi, 0.615). 

3. It is assumed that strain-rate effects can be approximated by an 
equation of the form of Eq. 2.1, eis follows'*’: 


2 ) 






i. 

f 


( 6 . 


^ere 




<Vk 


“ static yield stress of the kth eleistio, perf c ;tly- 

plastic mechanical sublaver (supersaript s ar J. subscript k 

» rate-dependent yield stress of the ktt elastic, 
P®^^®®^y~Plsstic mechanical sublayer ’(subscript k) 

• dg* Z 

^ ^ *dt" ^ X ^ longitudinal component of the rate-of -deforma- 

tion tensor (also called D^. in Eq,, B.IC of App. B) 

For illustrative purposes, the material strain-rate constants d and 
p for aluminum cited in Ref. 25 are used: d=6500 sec"^ and p»4. 

4. Thicicness changes are accommodated approximately. 

5. The change in length of the reference axis because of finite membrane 
strains is included in the expression for the changes-of-curvature 

in toe ^^angian strain-displacanent equation. 


Equation 6.2 is similar to Eq. 47 of Ref. 47 whito dealt with multiaxial 
conditions. 
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To illustrate the effects of these aodiflcations qpon the predicted 
transient 2-d re^Mutses of sijiple s t r uctur es, calculations were carried out 

two severely- loaded beams 1161 : (a) erplosively-impulsed CB-4 *na 

(b) steel— sphere-impacted beam ^ecimen CB-18. 

6.2 Explosively-Tirmised a«»«n CB-4 

Specimen CB-4 was an i n it i a l ly -flat 6061-T651 beam of 0.102-in 

thickness, 1.497-in width, and 8.005-in span with both ideally clasiped, 

and subjected to u n i f o rm ia^ulse loading over its entire width for a 

1.8- in spanwise region centered at midspan; ea^erimental data for test 

^>ecimen are given in Bef. 16. For this case, the iapulse loading produced 

essentially a u n i f o rm initia l lateral velocity of 10,590 in/sec over the 

i^mlsively-loaded portion of the beam. For anal^is, the half span of 

specimen CB-4 was modeled by using 20 equal-length 4 OQF/node finite elements, 

and synm»try conditions were i nqa osed at mid^pan. Four spanwise ca,ie,f;^ n 

stations and four d^thwise G au ss ian stations were used for the volume 

mmierical integration of tte finite-element equations. A consistent mace 

(CM) matrix was eoployed for each finite element. A rima increment sire At 

of 0.25 microseconds was used (equal to 1.6/w , wbere u is the 

frequency of the discrete numerical model) . The aluminum material was treated 

as behaving in an EL-SH-SB fashion, with material rate const;mts d=€,S00 sec~^ 

and p=4. The mass per unit i n i ti a l volixae p vas taken as 0.25384x10*^ 

2.4 ® 

(Ih-sec )/in . Response predictions for ^>ecimen CB-4 were carried out for 

the above modeling and conditions for: (a) the new mechanical sublayer 

fitting procedure of vs. c* (Subsection 6.1, item 2) and strain-displace- 

ment relation Type F cf Jqipendix B and (bj the former sclieme (aiploying stress 

(Tj^) vs. strain (Y^j^) Fit A of Subsection 2.2 and strain-displacement relation 

Type C of Appendix B and Ref. 20). 

I n d ic a t ed in the following tabulation are the cenpaurisons of these two 
P^®^crtions with each other (emd/or versos ea^>eriiiiental data) as shown in 
the indicated figures for the time histories of the longitudinal Green strain 
tensor can?)onent? on the upper (non-loaded) and/or the lower (in?)ulsively- 
loaded) surface at variems spanwise stations of beam specimen CB-4: 


On these figures is used to denote the longitudinal con^Jonent of this strain 
since this is the notation used in corresponding earlier figures in this report. 
However, this would be in the more systematic notation of Appendix B. 
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Tine Histories of*Yn Surface 
OPPer ( 0 ) and/or Lower (L) 


Figure 

station lx 1 (In) 

Predicted 

Measured 

6^ 

0. 0 '(Midspan) 

0 and L 


65b 

1 . 4 ' 

0 and L 

0 

65c 

2.2 

P 

0 

65d 

3.0 

0 

0 


With the exception of the oidspan station (x-0) , all other stations 
coincide with finlte-elasent nodes, and the associated strain value plotted 
is the average value given by the two elements at the nodal-juiUion location. 

It is seen that, of the spanwise stations shown, the major differences 
between the two procedures occur at the siidspan station in, where the 
finite-strain formulation shows that between 150 usee and 500 ysac the lower 
(loaded) surface experiences larger strains than the upper surface while 
the former "small-strain" formulation indicates the opposite behavior. Also, 
at this midspan station, the strains predicted by the finite-strain procedure 
are considerably larger than the strains predicted by the small-strain pro- 
cedure. At the other stations, where smaller strains occur, the differences 
***^®*® predictions are correspondingly smaller. 

Shown in Pig. 65e is the spanwise strain distribution at t-300 usee 
from x-0 in (midspan) to x-4.00 in (clamped end), of the upper (non-loaded) 
surface. This time instant is taken as typical, since the strains have 
already achieved their peak and about 97% ^the initial kinetic energy has 
been transformed into strain energy by that time. The strains predicted bv 
the finite-strain formulation are larger than those predicted by the smalll 
strain-formulation with the exception of a region-at £he' 
loaded zone (x-0. 9 in) and a region at the middle of the half-span (x-2.0 in 

tox-2.4in). The nodal strain discontinuities typical of the 4 DOF/node 

finite element (employed in the JET 3 and CIVM-JET 4B programs) are evident 
from the graph. This assumed-displacement finite-elaaent model involves 
cubic polynomials in the assumed-displacement field fer v (the axial displace- 
ment) and w (the lateral displacement) . The degrees-of-freedom (OOP) involved 


See the footnote on the preceding page. 


125 


at each end of the finite-element are thm 

, are the displacements v and w and the 
<iisplacement gradients" y » — > !L mnA rfi « ^ ^ 

^ These degrees-of.fr-ed«» 

provide continuity of disDlscenijin^ ^ ^ ^ -edcm 

(?" - X * V2 yC ' (V and w) and continuity of «»pteane strain 

/2 X V2 X + V2 ♦ ) but the be^ag strain lz< - ci(- #) (Uv) ♦ * iX,, 

IS not continuous at the nodes since and h are not h ^ 

s^nem ot^.< • 3 t) *“* ?n <*egrees-of -freedom. 

«■ »■».!» ooi, ^ 

«»« .. a., nod., onl^r cowiouit, dlspUc-... d..l„U« W 

a * i ^ - 5 “- «»*- d.ri..tiv.. - 

X ♦ Is of th. ,soo»i ord« (or ,o.dt.tlc) for «, initlsUy.,tr.l 5 ht 

sgrM of ths polyooiiial iimalrsd in th« reprs.Mt.tion of the first 
deriva^ves of the displacement gradients ^ > l5l ^ j 

ii . ij; . ± . Si 3,2 - 5K <r' “> 

9n^ *'> “ “* «" l«itiau,-.trai,ht 

^s. ^.in, the . noP/noi, oMio-Mhio ei-ent,. fr- Pi,, .s, u i, ohMrsM 

n2i?”V ““ -«»«i^.tion 

IS (mainly) either quadratic or linear. 

It i. .i» Observed that the lar,«,t discontinuities occur at locations 
where bending stradns are largest* at th» an /4 ^ -i. • 

(«-0 9 inl „d ,r .. . M the Of the is^ulsively-ioaded ton. 

to ti! r r *“ ““ “ ““ “* '«•» t" 

J4.0 i«). At the cianped Mne, a very iar,e strain discontinuity is 
ent. The reason for this is that this regioST in«,ii,s-hl<,h'£;;;i; " , 
»nlinea:it,. the strain discontinuity at the cias^d ton. i. significaitly 
larger ,»th the finitMstrain foMulation, rt^ch involve, a „„ ^ 

representation of the behavior tl»n the -ssall-strain- fonsulation. it is 
evident that a fin» „sh of finit. eiasents is needM i, this clanped-tl 
tMion to represent accurately this nonlinair behavior. Hoover, tin. and 
^ restriction, have prevented a Mr. thorough study of this at this 


'fScS^; Sn£Mrty-rr^r;:,%'M.2T^^^^^ 
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Finally, the predicted transient midspan transverse displacement w is 
shewn in Fig. 66 . It is seen that the finite-strain fomulation and small- 
strain formulation predictions are in fairly good agreement :*ith each other. 

Thei casting time required ^or the two formulations for «q>losively- 
iaqnilsed bea CB-4 is displayed conveniently in the following tabulation 
(for eOOO time steps with a time step size of 0.25 microseconds; all runs 
were conducted on an IBM 370/166 cooqniter): 

**<>• of No. of Gaussian Total No. 

Formulation FB Sta. per Elem. of Onknown 

Spanwise Depth 00 ^ 

Small Strain 20 4 4 79 

Fiip.te Strain 20 4 4 79 


Strain-Displ. Hass 
Formulation Relation Matrl. 1 i 

Type 

Small Strain C CM 

Finite Strain F CM 


No. of 

CPU 

CP0(min) 

Cycles 

Time 

(min) 

(DOF) (Cycles) 

4000 

8.63 

27.3xl0~® 

4000 

11.07 

35.0x10"® 


The Meets on CPO tl« of the length, expteeeionn need end 

reunited for the finite-strain cUenUtions are evident frot an inspection 
of the last colisnn. 


6.3 Ste el-Sphere-lmpacted Beam CB-18 

Steel-sphere-lmpacted 6061-T6S1 aluminum beam "sp^lmen CB-18, the 
geometry and impact conditions fnr wh i ch have beu r#scribed in Subsection 2 1 
has been analyzed only for EL-SH material behavior (a) the (new) finite " ’ 
strain procedure (employing the stress (t^) vs. strain (e*) mechanical-sublayer 
fit of Subsection 6.1, item 2; and strain displacement relation Type F of 
Appendix B» and (b) the former scheme (enqiloying stress (x,^) vs. strain (y ) 

Fit B of Subsection 2.4 and strain-displacement relation Type B of Appendi^ B) . 
The beam was modeled with 43 equal-length 4 DOF/node elements; locally 
perfectly-elastic impact (e»l) was assumed. Three spair^ise Gaussian stations 
and four depthwise Gaussian stations were used for the volume numerical 
integration of the finite-element equations. A diagonalized (lumped) mass 
(DM) matrix for each element was employed. A time increment size At of 



0.5 microseconds was used (equal to. 1.16/u^, where is the 

frequency of the discrete numerical modal) . Sunmarlzad*L the foUowing 
tabulation are the comparisons of the two predictions with each other (and/or 
versus experimental data) shown la the indicated figures for the time histories 
of the longitudinal component of the Green (-Lagranglan") strain tensor y on 
the upper (non-lmpacted) ai^or the lower (impacted) s^ace at various ^ 
wise stations of beam spec^eh CB>18: 


Time Histories of y. . on Surface' 
n (Upper and/or L (Corner) 


Figure 

station |x|'- (in) 

Predicted 

Neasiired 

67a,b 

,0.0 

U and L 

— 

67cpd 

0.3 ^ 

U and L 

_ 

67e 

0.6 

D 

D 

67f 

1.2 

0 

0 

67g,h 

1.5 

0 and L 

0 and L 

67i,j 

3.0 

0 and L 


67k, i 

3.7 

0 and L 


67nip n 

4.0 

U and L 



Location x=0 in (at the midspan of the beam) coincides with the i" id « q > a n 
Gaussian integration station of a finite eleanent. Location x « + 4.0 in is 
at the clamped end of the beam and coincides with a finite element node, at 
which claaq»ed-end conditions have been imposed (namely that the displacanents 
V and w and the lateral-displacement gradient are zero) . All other stations 
occur at locations intezmediate between the end and the midspan of a finite 
element, and do not coincide with spanwise Gaussian integration points. Also, 
measured permanent strains are indicated on these figures where avadlable. 

These figures show that the strairs predicted (a) by the current 
"finite-strain procedure" and (b) by the former "saall-strain procedure" 
agree reasonably well with each other and/or with experiment at all of 
these stations except x*0, 3.7, and 4.0 in. Large strains do occur at both 
x=0 2 U)d xM.O in; also, the occurrence of large strains at x=»4.0-in exerts 
a distinct and pronounced effect at "neaurby station" x=«3.7 in (located in 
the element adjacent to the finite elenent at which the claIl^>ed end condition 
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has been ixi^sed) • Although the calculations have been carried out for only 
900 microseconds, it appears that the current "finite strain procedure" would 
provide better permanent strain coi^>arisons with measiirements at all spanwise 
stations (if carried out long enough in time) than by the former "small*strain 
procedure" . 

Figura 68 shows that the time histories of the midspan lateral deflec- 
tion w from these two predictions for beam tt-18 are very close to each other. 

Finally, the time histories of the support reactions M , S , and F at 

y z X 

station ac-4.0 in are shown in Figs. 69a, 6pb, and 69c, respectively, for 

these two-predictions. The agreement between these two predictions is very 

good for the longitudi n al support reaction .force F^ (associated with the 

membrane strains), but one observes some differences in the transverse support 

reaction (shear) force S and large differences for the support reaction 

z 

bending moment M^. These differences are believed to be caused by the fact 
that the expressions of CIVM-JET 4B for the bending part of the strain are 
valid only for small rotations and small strains, while the finite strain 
version of the program does not have this restriction. Of course, the 
st^port reaction bending moment is most influenced by the bending part 
of the strain-displacement relations. 

The ccmputing time required to analyze steel-sphere-inpacted beam CB-18 
by the two procedures, under otherwise-identical conditions, is conveniently 
displayed in the following tabulation (for a time step of 0.50 microseconds; 
all tans were conducted on an IBM 270/168 cooputer) : 


• ' 

No. of 

No. of (Russian 

Total No. 

Formulation 

FE 

Sta. per Elem. 

of Unknown 



Spanwise Depth 

DOF 

Small Strain 

43 

3 4 

170 

Finite Strain 

43 

3 4 

170 



Strain-Displ. 

Mass 

No. of 

CPU 

CPU (min) 

Formulation 

Relation 

Type 

Matrix 

Cycles 

Time 

(min) 

(DOF) (Cycles) 

Small Strain 

B 

DM 

2250 

5.11 

13.4xl0~® 

Finite Strain 

F - 

DM 

1850 

6.81 

21.7xl0~® 
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SECTION 7 


SCMHART, OONCLOSXONS, ANt) COMMENTS 

7,1 siaaary 

-The-type of response exhibited by a given protective structure to engine 
rotor fragment impact attack depends, among other things, t^n the lateral 
(axiail) dimension of that structxire coiopared with the corresponding diioen- 
sion^s) of the attadcing fragment(s) . If these diinensions are coo^»arable, 
the structural response of the protective structure is essentially the same 
at all of its axial stations at a given circumferential location; accordingly, * 
this structural response is termed "tw^diaiRnsional" (2-D) • However# if the 
axial length of the protective structure is substantially greater than the . : 
corresponding dimension of the attacking fragment# the structure exhibits a 
three-dimensio]^ (3-D) type of structural response. The present report is 
devoted to considering only the former (2-D) category of structural response 
problems. In particular# the present report describes studies performed to 
aissess the applicability# feasibility#^ and. utility of employing a 2-D struc- 
txiral response conputer code (a) for predicting the trauisient large-deflection 
elastic-plastic structural responses of 2-D ductile metal protective struc- 
tures to fragment upact and (b) for carrying out variation-of-parameters 
studies to assess the consequent inpact-induced structural responses of such 
fragment-impacted structxires. These studies involved the use mainly of the 
CIVM-JET 4B conputer program which permits analyzing the 2-D structural 
responses of single^ layer cxirved Bemoulli-Euler beams and/or rings which 
are subjected to inpact attack by one or more idealized non-deformable 
fragments * Some calculations were also carried out with the CIVM-JET SB • 
program which has featiires and capabilities similar to the former code except 
that it can represent Bemoulli-Euler 2-0 beam/ring stinictures that consist 
of (a) a single layer or (b) several layers of different materials which are 
hard bonded at each interface. CIVM-JET 4B employs the central-difference 
timewise operator while CIVM-JET SB uses the Houbolt timewise finite-difference 
operator. The spatial description of the structures in both codes is 
accomplished by the use of finite elements. In general, one seeks to enploy 
the smallest number of finite elements# which will result in converged or 
accurate transient response predictions. 
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value for the idealized fra 9 ment radius one can reasonably escpect that 
the predicted peak response should coa^are favorably with that measured; this 
was the principal aim of this particular study. Examined also was the effect 
upon the predicted responses of the containment ring and the fragments of 
friction between each fragment and the containment ring. 

Dimensional analysis considerations to identify pertinent dimensionless 
parameters were en^loyed in examining experimental data obtained (a) by the 
NAPTC on engine rotor;^burst inpact tests on steel containment rings and (b) 
by the MIT-ASKL on aluminum beams subjected to steel-sphere ispact attack. 

The use of the dVM-JCT 4B conputer code for carrying out parametric 
structural response studies has been illurtrated for both fragment-containment 
structure and fragmentrdef lector structure. 

Finally# since the CIVH-JET 4B conputer code as actually isplemented 
accommodates large-deflection elastic-plastic transient structural responses 
but only small strain, some modifications to the analysis/conputation proce- 
dure have -been developed to alleviate this restriction. Some preliminary 
results from this study were obtained and are included. 

7.2 Conclusions 

Within the context of enploying the present type of conputer code(s) to 
represent: (a) 2— D Bemoulli-Euler large— deflection, elastic— plastic transient 

structural responses such that the maximum strains are not laurge (i.e., less 
than roughly 8-10 per cent) and (b) the fragment (assumed to produce the 
impact-induced structural~respolTse) idealized as a non-deformable circular 
2-D fragment, the following observations are made from the present study .and 
from related similar MIT-ASRL studies [4,5,15,17,18,20]:* 

1. For the analysis of single-layer** structures subjected to fragment 
inpact, the CIVM-JET 4B code is more efficient that CXVM-JET 5B. 

Since CXVM-JET 4B uses the explicit conditionally- stable timewise 
central-difference operator together with the very compact and 
efficient unconventional form of the equations of motion whereas 

The listed conclusions are based upon the present study unless indicated 
otherwise. 

» ’ 

However, for the analysis of hard-bonded multi-layer Bemoulli-Euler 
structures, only one of these two codes (CIVM-JET SB) is applicable. 
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CIVM-JET SB uses the unoondltionally-stable (for linear systems) 
timewise Houbolt finite-difference operator together with the less- 
oompact conventional form of the equations of motion, it was specu- 
lated that the latter co ntoi n ation might be more efficient mainly 
because (hopefully) of a l^er allowable time step size At than for 
the former. This has turned out not. to be true because the At is 
found to be limited by the "momentum— transfer and correction proce— ■ 
dure" employed to account for momeiitum transfer from the fragment 
to the target structure at j each of a succession of ia^cts. 

All further comments pertain to the - ^plication of only the ClVM-JET 4B 

coo^uter code. 
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TO insure numerical stabiuiy, the time increment size At for solu- 
tion of typical nonlinear transient response problems be 

selected such that At < 0.8(2/u ) where u is the hiqhest 

natural frequency of the mathematical model of the system for small- 
displacement linear-elastic behavior (17] . No numerical instabilities 
have been encountered when this guide has been followed; however, 
either local or global numerical instabilities have been observed 

for At > 0 . 88 ( 2 / 0 ) ). 

max 


3. At each of a succession of inlets of the non-deformable fragment 
against the 2-0 target structure, inqjulse-momentum relations are 
enqjloyed to predict the velocity decrement suffered by the attaeJeing 
fragment and the associated velocity increment acquired by an intact- 
affected region of the target structureT Alternate schemes have 
be^' used to estimate the size of this intact-affected region. 

Since physically the size (or spanwise length) of that intact- 
affected region must be independent of both the time increment size 
At used in the calculation and the size (length) of the finite 
elements used to model the structure, elementary stress-wave-propaga- 
tion arguments have been used to estimate the distance I. from 

eff. 

the in 5 )act station to each end of the in5>act-affected region to be 
approximately twice the thidmess h of the structure. Accordingly, 
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in regions of expected iispact, the finite element lengths should not 
exceed 2h; smaller finite-element length values should be used, if 
needed, for other nonr impact considerations such as solution 
convergence. 

4. TO be consistent with certain assumptions eii 5 )loyed in the approxi- 
mate procedure for predicting the velocity increment received by the 
impact-affected region of the structure from a given ispact, a 
diagonalized mass matrix of the finite-element TOdel of the structure 
is employed* Various means of achieving the diagonalized mass matrix 
were explored. The main effect of these alternate choices was to 

affect 0 ) which, in turn, influences the selection of the size of 
max 

At required for numerically stable calculation:^, (see Subsection 
2.5,2). 

5. The inclusion of plausible values for the coefficient for friction y 
between the target structure and the idealized fragment reveals very 
little effect upon the transient response of the inqpacted structure 
but can affect significantly the rotational motion of the idealized 
fragment [15, 18, aind the present report]. 

6. Conclusion 5 ia?>lies in part that the response of a fragment- 

impacted structure is only slightly affected by friction and/or by 
the rotational motion (or rotational kinetic energy of the attacking 
fragment) but that the structxiral response is dominantly affected by 
the translational velocity (or tr 2 Uislational kinetic energy) of the "" 
attacking fragment. This observation is consistent with the find- 
ings in this study. ~ « 

7. From the present series of calctilations of the transient responses 

of fragment- imparted 2-D structures, it is noted that the peak 
structural response is affected: (a) only slightly if one ignores 

strain-rate sensitivity effects for relatively rate-insensitive 
materials such as 6061-T651 (or T6) aluminum* but (b) very signifi- 
cantly for rate-sensitive materials such as mild steel. Hence, it 

^However, for aluminxjm beams ind/or rings subjected to severe impulsive load- 
ing, strain rate effects are found to have a very significant effect upon the 
predicted responses of these structures. In these cases a laurger volumetric 
proportion of the structure experienced higher strain rates and for a longer 
period of time than in the present fragment- impacted structures. 
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is isqportant to determine vrtxether or not the mechanical behavior of 
the material under consideration exhibits significant rate sensitiv- 
ity and, if so, to include that description quantitatively. Further, 
uniaxial material stress-strjaln measurements (or data) are required 
from material coupons whose lengthwise axis is rUgned with the main 
direction of tensile or compression stress e:qpected in the 2-D frag- 
ment container or deflector structure ~ since mechanical anisotropy 
present could lead to the inadvertent use of iaproper mechanical data 
in the analysis with attendant misleading results. 

8. The presently-ayailable experimental data on the responses of both 
beam and ring structures subjected to fragment Impact are inadequate 
to represent both (a) 2-0 structural response and (b) 2-0 fragment 
impact attack — up to the threshold rupture condition of tte impacted 
structure; conditions short of this rupture state are widely-agreed to 
be Inadequate for containment assessment purposes. Since suitable 
techniques are now available for conducting such tests and making the 
essential measurements, experiments to ramedy this data deficiency are 
recommended. 

9. The fact that the present conputer codes (4,5] as act^^Tiy impio- 
inented can accommodate Is :ge-def lection, elastic-plastic transient 
2-0 structural responses but only small strain means that although 
large numerical values of transient stradn can be and eire predicted 
from these programs, those predicted large strain v^ues (i.e. , 
greater than about 8-10 per cent or even less) can not be correct. 
Steps are described in Section 6 to remedy this deficiency and • 
preliminary results show that significant differences in transient 
strain predictions in certain portions of the structure are revealed 
by this more comprehensive treatment. However, the predicted 
transient deflections (displacements) are almost unaffected by those 
ixtprovements . 

10. The present calculations and comparisons with NAPTC experimental 
measurements of a steel containment ring subjected to T5 8 turbine 
rctor tri-hub burst attack in NAPTC Test 201 indicates that the use 
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of an id galized non-deformriale fragment of plausibly-selected size 
(as an approximation for each of the deformable bladed-disJc fragments), 
for use in the CIVM-JET 4B program can provide predicted peak and 
permanent ring deflection configurations in reasonably close agree- 
ment with experimental observations. The current CIVM-JET 4B pre- 
dictions of large (or finite) transient and permanent strain for 
. this case can not be correct but the breadboard CIVM-JET 4c program . 
now available should remedy this deficiency. However, the available 
measurements of transient and permanent strain for this case are 
' inadequate to permit a definitive assessment of the accuracy of 
either of these predictions. 

Finally, the foUowing two additional observations are offered: 

11- In the present model of fragment/structure impact and interaction, 

impulse-momentum relations are-used to estimate the momentum transfer. 
Hence, the velocity increment imparted to the "impact-affected por- 
tion" of the structure for each of a succession of ingjacts depends 
upon the ratio of the mass of the fragment to that of the impact- 
affected region of the structure. Hence, in seeking to define 
scaled conditions to preserve the peak structural-response strain 
which is induced by fragment in^iact, one must preserve certain energy 
absorption and "incremental momentum-transfer" parameters as noted 
in Subsection 4. 3.2.5. 

12. examination of NAPTC experimental data in which spinning rotors 

were caused to rupture into n equal-size bladed-disk fragments 
(n=2,3,4, or 6) and subsequently to impact' against containment rings, 
shows clearly why the bi-hub (2 fragment) attack at a given burst 
rpm produces more severe ring response than observed fcr 3-, 4 -, and 
6-fragment attacks ~ even though a 3-fragment burst results in 
fragments with the greatest amount of translational kinetic energy 
and hence were often viewed as posing the most severe threat to a 
containment ring. For n=3 and greater, the deforming blades 
cushion" the impact, and peak ring response is reached while only 
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the "smooth contour" of the fragment is in contact %rith the ring; 
aH-of-these-transiently-and-pexiaanentiy deformed rings 

exhibit a smooth contour* However, for the 2-fragment case, a sharp 
comer of each of these 2 fragments attacks the ring and 
extroaely severe and XocaLized bending (and perhaps some shearing) 
of the ring — resulting in local fatilure at significantly smaller 
attack kinetic energy levels than for the n>2 cases* 

i ■ 

7*3 Ccaanents 

Included here are indication^ of some analysis development needs , 
recommendations for experiments to remedy is^rtant deficiencies in e^eri- 
mental data on severe fraginent- impact-induced structural response, and some 
related observations* 

Analysis Needs 

Within the framework of 2-D structtiral response, of containment or 
deflector structure subjected to 2-D non-deformable fragment ia^jact attack, 
the type of nonlinear transient response analysis and computer code discussed 
herein which includes large-deflection, elastic-plastic behavior but only 
smal l strains (as actually implemented in the computer code) should be 
extended" to "incliide finite strain effects; the preliminary work to g end 
as discussed in Section 6 and Appendix B should be carried forward to 
completion, including cosparisons with appropriate well-defined experimental 
structural response data for a range of maxinuTin responses extended to and 
including material npture* 

Although under the conditions representative of engine rotor fragment 
impact, the type of ductile metal containment structure npture observed has 
been, in almost all cases, noted as being "tensile rapture" without clear 
evidence of significant transverse shearing effects, it may be useful to 
extend the theoretical prediction method to include transverse shear deforma- 
tion effects so as not to preclude the proper treatxDent of this behavior 
where and when it occurs* This inclusion, however, will tend to make the 
predictions more conplex and expensive* 
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The necessity of talcing into account the fact that fragmts such as 
rotor blades, bladed-disk segments, or bladed-rim. segments often deform 
severely during in?>act attack against containment or deflector structure 
is a matter %ihich has been weighed and debated considerably. One school 
of though holds that taking into account transient fragment deformation would 
involve a great amount of extra confutation and book-keeping because, it is 
argued, one can predict the peak cont^u.nment-struicture response with adequate 
accuracy (in view of the various uncertain^es present in many other aspects 
of the impact/response problem) by representing each attacking deformable 
fragment as a simple non-deformable fragment of plausibly-selecrted size, mass, 
as well as pre-impact translational and rotational velocity — although there 
is clear recognition that the t ime history details between the actual response 
and that predicted by the idealized model are certain to be different. Another 
school, while conceding the plausibility and the analysis feasibility and 
efficiency of the above-described viewpoint, argues that fragment deformation 
and possible subsequent multiple fragmentation will alleviate and distribute 
the ispact-induced loads on the structure such that the actual response of 
the containment structure will be less severe than one would predict by the 
above-indicated idealized rigid- fragment emal'ysis. However, no proposals for 
deformable-fragment models to represent specific categories or conditions of 
importance for containment-structure design are offered, but sometimes there 
are i^lied wishes and suggestions that engine rotor failure presents no real 

safety problem and should be ignored. Although aircraft losses and fatalities 
caused by or stemming from engine rotor fragments have been small to date, the 
annual rate of uncontained engine rotor fragments in commercial avaiation has 
persisted at about one per million engine hours since about 1962 as reported by 
M£uigano et al. [37], McCarthy [42], and Gunstone [55]. 

It should be noted that both deformable-fragment and rigid-fragment 
analyses have been carried out successfully for impact-attack by a single blade 
of a T58 turbine rotor against an aluminum containment structure [15,17]; the 
amount of computing required and the type of impact/interaction analysis used 
for these two cases are very different. Thus, either approach can be employed 
if required. Thus, if a persuasive case can be made for including deforming- 
fragment behavior in the analysis of aircraft engine rotor fragment impact 
attack against containment or deflector structure, this can be done — but with 
an attendant addition in computational complexity and cost. 
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It should-be noted that methods and computer programs for ore^'etln 
l»~l. «hl,h under,, 3-0 .rmct„u r..p,n.. ■ „-, . ,„ 

as actually implemented these coa,)uter codes acconmod «. ^ ^ * However, 

strain n«nr.<. < 4 . i ^ accommodate properly only small 

in ! "C-mnded thnt thl. «,rl. b. «««d«, t, finit. nttni., 

d.f,™7^ *” ••P*«3 ,d th. rr«.l„t Ur,. 

«1«9U-Uy,r ducriu „tn3 pr«.dU„ ,rr,c- 

«a r.«>i,«j,. u n. ^ ,9 

. ^ r, predicr me rru.i»r rerpon... ,rl7Zr,::r' 

17“L “tncb. incudin, conditions 

Which prodoc 5«,„.„ti.i Uilor. of th. rnrloos str«:ur.i Uy.rs - 

7lT^;7 T““' ■” “<• W“h uy„ cn h... .sbstnn. 

Ooth.»pit. ud hocc. Of th. conpunitp .f tbis 7.^.7 7.’ 

uTrTTi ^ P--»d witho,; d.Uf. 

P With .nd closoly coordinnt.d ,ith a„t effort rt»ad u „ 

-Ui i^ct-mdccd r.,po„. pro,r„ for ™u-d.fi.ed .np,rU»t,i oJT 
tions and careful detailed measurements. 

Ejqperimental Needs 

Of 77“”'^ -W> deflniti.. dat. for th. .fflof„r 

in d nscurecy of proposM „thods for predictln, the Upact- 

i~ac« responses.., incUdin, thr.shoid roping, ductU.-„t .3 proTul 
structures, th. follc«i„, erperUent, are r«»-.nded. 

(.) TO ,«..rat. well-defined 3-D structural 

perpendicular 2-D Upact b, a soUd steel cylindricei bod. 

Transient strains should be oeasured on each surface at various 

spenwis, stations. T»t « suitable r^ion centered at .idspen and 

for reoions close to each clpsped end of the bene. , pattern of 
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mechanically Ilghtly-scirbed closely-spaced grids should be provided 
to pen&it the determination of the values and distribution of the 
permanent strain since these are regions where significant strains 
are expected* 

(b) To generate well-defined 3-D structural response, data for efficient 
theoretical-experimental correlation evaluation of the adequacy of 
the proposed finite-strain extension of the analysis of Ref. 10, 
6061-T651 aluminum beams like those of Ref* 16 should be subjected 
to midspan perpendicular 3-D impact by a solid steel sphere. 
Measurements similar to those described under (a) should be made; 
of special inq>ortance are the .proposed permanent strain determina- 
tions from the mechanically-scribed grids* 

This currently xmavailable experimental information is essential for a mewing- 
ful evaluation of the adequacy and accxiracy of the cited predictions methods 
to predict structural responses up to and including material rupture . 

More cCT^rehensive experimental data than are now available on contain- 
ment ring response to engine rotor fragment in^ct~ would be valuable- both in 
its own right and for checking the adequacy of prediction methods such as 
CIVM-JET 4B or alternate approaches. In particular, it is recommended that 
4130 ceLSt steel containment rings much like that in NAPTC Test 201 be subjected 
to T58 tiirbine rotor tri-hub burst intact attack such that the peak intact- 
induced response will be (a) much closer to ring rupture than in NAPTC Test 201 
and (b) slightly exceed the threshold ring rupture condition. This could be 
achieved at the nominal rotor-burst kinetic energy level of NAPTC Test 201 by ' 
using somewhat thinner containment rings. The recommended measurements 
include: : . ' . . 

(a) Transient strain at a number of circumferential locations on the 
ring's outer-surface midlength location, and at the midsurface on 
both end faces at several circumferential locations. 

(b) Finely-spaced lightly- scribed closely-spaced lines on the entire 
outer circumference of the ring should be provided, and both their 
pre-test and post-test spacings measured to permit determining the 
permanent outer-surface strain. 
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(c) 


al»b-.p..d photo,r.ph 5 rtould b. »d. to obootv. u„ .... 

ttioo. ct tb. „d of tb. ci.^, pboto,r»bf. 

to tbot «0blt«d io NUTC Toot «7 should bo sought. 

Ibis «uld ptortd. tspras«>t.ti» «>d dofllod dsts o» to»tiOo.«ot-tlo, 
rsscoosu to OO.T1SU dsfontsblo blsdsd-dlsb ftsgbsnt ij^uct sttuck 

rot osd, Of tbss. tl.p.. Of ««,l.».ts, U.1..1U1 st.«o to»Us sttsss- 
tss^ito n,tut. Should b. »„1« iu. «, co.tbl»..„t-.tr»ctut. coupoo. 
.bos. ^ otl«.t.tlo. Of «ob 1. ,1, pstuuoi so tb. spso^s. ot olsouJ^ 

.twtl.1 dlt.cUo», (2) psq^oulsr t^otl«t.Uon (1). m (1) .t 
i«t«^.t. o.l«,ftio» sb,!., tb... ^ 

not tb. «torl.l 1, isotrspic. Such tMt.. «^d b. c«ri.d out .t (.) tb. 

’n”Trj“ nnot-bo.d sp«d, tb... .ps.d. 

° t.cotdM. UgbUir-sorlbn) si.olume.1 grids should b« provldsd on 

~cb sp..^ to p«.it d.t.„inlb, tb. distribution „d vUu. of tb. p.t«„.„s 
ou«r-.urf.o. strsl, on „cb spools,... „.„or««.ts should «so b. ,»l. to 
P.mit dstor^Mo, tb. pr.-t.st „d post-tsst or„....ctionU sr.« of .«b 
specioen along its axial test-region length. 

General 

Although b.pro™.«.ts ii tb. cpobiliu., of tt. CIV»-JEI 4 B cot=ut.r 
«._n..d«i.«d vr. pl«u.«i. It is b.li.v« tb.t tb. CIVM-OI, « ot.,.,.^ 
nod. b. Of sigiUflcMt „.i,t«o. to industry for t»ny psr«»trio «d 
pr.Usun«y d.,15-., puspo... 1. rotor burst protsction studios sod in tb. 
dtsign of frog«nt-oont.i™»„t «d/or fr.gM.t-d.fl.ctor structuros 

it is .ugg..t.d tb.t M offoctiM prooMur. for tt. d.sigMng of fr»- ' 
~nt contM.,M.t or fr.gt«,t dofloctor structur. to . giys. typ, of „,i.. 
rotor fr.„».t .tucb »uld b. to co«iuct . f« ..ll-iostruMntM full-sd. 
tost.. Mtb tbM. prolisdnor, d.M m guidMoo, . cospaUr cod., such m - 
«. for sMitplo, could b. usm to curry out p«s,*tric „d trsdMff 
studios. Of courss. .11 otb„ rolovMt «tplricl «d s«,u-«.piric.l info™- 
Uo. sb-uld a lso bo .ppliod i. dSMloping . proposod protMtivMstructur. 

■inlsnss of suS^tSL?' to^ddiU^nf'SjrsSSd “* 

tb. str^-t.t.-d.pMdMt bohavior of tb. Mtortif^J^ “ ?*“t»in. 

•ultlMisl stress conditions, including tb. mSm sSfS "■* 

ond stru. r.eos up to .t lout 3.000 I.c'* — m f 

.doget. ..r,ss, strUn, stroi.-r.t. doscription fH “ 
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design. Then several carefully-*selected full-scale proof tests should be 
conducted to demonstrate the adequacy of the final design and to confirm 
that important effects had not been overlooked. Such a procedure should 
be much-nore-tiae-and-cost effective than a largely e3q>eriittental approach. 
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TABLE 1 

COMPARISONS OP RESPONSE PREDICTIONS FOR VARIOUS MASS MATRIX 
- MOD E UN GS-OF-IMPOLSIVELY-LOADED CLAMPED-ENDED BEAM CB-4 AT 
100 MICROSECONDS AFTER APPLICATION OP THE INITIAL VELOCITY 


Mass Matrix 
Modeling 

Consistent 

Mass. 

Diagona.'. 
Method 2 

Diagonal 
Method 3 

Midspan 

Displacement 

0.7099 ; 



W (in) 

0.7141 

0.7114 

Max. Gaussian 

i 



Si.a« Upper*Surface 


_ 


Strain (per cent) 

10.04 

10.12 

10.07 

Max. Upper-Surface 




Nodal Sta. Avg. Strain 
(per cent) 

10.39 

10.49 

10.39 

Kinetic Energy (in- lb) 

809.1 

813.2 

812.5 

Elastic Energy (in- lb) 

30.5 

29.3 

31.0 

Plastic Wbrk (in-lb) 

1118.9 

1161.4 

1162.0 

Initial Kinetic 
Energy (in-lb) 

1958.5 

2003.9 

2005.5 











AAOIA A 


DATA W^CTERIZiNG NAPTC TEST 201 FOR T58 TOPFU^b ROTOR 
TRI-HUB BURST AGAINST A STEEL CONTAINMENT RING 


Containment Ring Data 


Inside Diameter (in) 
Radial Thickness (in) 
Axial length (in) 
Material 

El^lstic Modulus (psi) 
4130 Cast Steel 


15.00 

0.625 

1.50 

4130 cast steel 
29 X io®r 


Fragment Data* 


t 58 Tri-Hub Bladed Disk Fragments 
***^®*^^^ Disk: A-286 Blades: SEL-15 


Outer Radius (in)- 

Fragment Centroid from Rotor Axis (in) 
Fragment Pre-Test Tip Clearance from Ring (in) 
Fragment CG to Blade Tip Distance (in) 

Fragment Weight Each (lbs) 

Fragment Mass Moment of Inertia about its 
CG (in lb sec2) 


Rotor Burst Speed (rpm) 

Fragment Tip Velocity (ips) 

Fragment CG Velocity (ips) 

Fragment Initial Angular Velocity (rad/sec) 
Fragment Translational I® (in-lb) 

Each Fragment 
Total for Three Fragments 
Fragment Rotational Iffi (in lb) 

Each Fragment . 7 ^ 

- • Total. -for Threes Fra^ehtS’ ' ' ' 


Applies to each fragment unless . specified otherwise. 


7.00 

2.797 

0.50 

4.203 

3.627 

0.0666 

19,859 

14,557.2 

5816,7 

2079.6- 

158,922 

476,766 

144,018 

432,054 









TABLE 4 

SDMMRRy OF dVM-JET 45 CALCULATIONS MADE TO ANALYZE 
4130 STEEL CONTAINMENT RING RESPONSE TO T58 TURBINE 
ROTOR TRI-HUB BURST ATTACK IN NAPTC TEST 201 


Calculation^ 
Identif. No. 

r^(in) 

W 

Simultaneous 

Impact 

THR-1 

2.555 

0 

Yes 

THR-2 

3.360 

0 

Yes 

THR-3 

4.203 

' 0 

Yes 

THR-4 

2.555 

0.3 

Yes 

THR-5 

2.555 

0 

NO^ 


a: In all cases: (a) 48 equal-lengtii c\ibic cubic elements with 

4 DOF/Node were used and (b) the steel ring material be- 
havior was modeled as EL-SH-SR. 

b: Fragment 1 was released first a nd-6p -mi-croseconds later, 

fragments 2 and 3 were released; unsymmetric 
response subsequently occurred. 




TABLE 5 

DATA CHARACTERIZING THE PRE- IMPACT FRAGMENTS IN THE NAPTC T58 
TURBINE ROTOR BURST TESTS AGAINST 4130' STEEL CONTAINMENT-RINGS" 


Quantity 


Mass of Fragment 
(lb-sec^) /in 


Mass Moment of Inertia 

2 

about Fragment CG, I^(in-lb*sec ) 


Distance from Axis of Rotation 

of Rotor to Fragment CG» r_(in) 

CG 

Distance from Fragment CG 
to Blade Tip, il^(in) 


Naninal Pre-Impact Translational 

Kinetic Energy per Fragment^ 

(KE) , (in- lb) 
ot 

Nominal Pre- Impact Rotational Kinetic 

Energy per Fragment^ 

(KE)^ , (in- lb) 

or 


Number of Equal-Size Fragments^ n 


6 


Ncaainal Translational Velocity 
Fragment CG, (in/sec) 

Blade Tip (in/sec) 


2 

3 

1.43xl0"^ 

9.22x10"^ 

.120 

• 0666 

2.266 

2.797 

4.734 

4.203 

161,050 

158, 200 

262,216 

115, 125 


- it 

4,746 

5,858 

14,661 

14,661 


Nominal Failure Speed ; - 

= 20,000 rpm = 2094 rad/sec 









TABLE & 


SOMMJUOf OP CONTAINMENT RING TESTS, CONDITIONS, AND MATERIAL PROPERTIES 
FOR THE NAPTC T58 TURBINE ROTOR BURST CCaTEAINMENT TESTS* 


Rin 

7 Width Conditions Studied | 



Test Identification 

Axial Width. 

Nunber of 


w^(in) 

Fragments 

NAPTC Test Number 

0.5 

2 

135,137 

1.0 

2 

139,131,133,142 

0.5 

3 

168,169,172,177,178,189,207 

1.0 

3 

126, 127, 128, 130, 132, 138, 139, 141,188 

1.5 

3 

192, 193, 195,196, 197, 199,201, 202 

2.0 

3 

170,173,174,176 

1.0 

6 

134,136,140 

Xden 

tification of Containment Ring Material Supplier 



Test Identification 

Supplier 

Ntanber of 



Fragments 

NAPTC Test Number 

National 

2 

1-29, 131T133V13S7I37, 142 

Forge (NP) 

3 

, 127, 128, 129, 130, 132, 138, 139, 141,168 



169, 170, 172,173, 174, 176, 177, 178 

NF 

6 

134,136,140 

• . ACIPCO ^ 

.3' ;; .c 

r * •* . -> 

r 188, 189, 192, 193, 195, 196> 197, 199 



201,202,207 

Nominal Mech 

anical Properties of 4130 Cast, Steel Cohtainmenc Rina Material"** 


a 

o 

\ V , E P V 




Si^)p^er 

(ksi) 

(ksi) (in/in) (psi) Ib/in^ 

NF 

80 

111 0.08 30x10® .283 0.3 

ACIPCO 

83 

121. 0.15 30x10® . .282 0.3 

Ref. 37 — subscript "o" refers to static properties. 






TABLE 7 


DATA CHARACTERIZING THE PRE-IMPACT PRAOfflJTS IN THE NAPTC_J.6S 
TDRBINE rotor BURST TESTS AGAINST 4130 STEEL CONTAINMENT RINGS 
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TABLE 9 

SDMMASy OP TESTS, CONDITlCaiS, AND MATERIAL PROPERTIES 
FOR THE MIT-ASRL SPHERE-BEAM IMPAGT-TESTS* 


Test Identification and In^ct Conditions 


Beam* 4uid 

Fragment Data^ 

Remarks 

Test NOe 


Initial Kinetic 
Energy (KE) (in-lb) 


CB-9i 

1896 

666 

Small Perm. 
Deflection 

CB-13 

2490 

1183 

Moderate Perm. 
Deflection 

CB-18 

2794 

1489 

Large Perm. 
Deflection 

CB-16 

2868 

1569 

Just Beyond 
Threshold Rupture 

CB-14 

3075 

1804 

Well Beyond 
Threshold Rupture 


a: Nominal Beam Dimensions are h=0.10 in, w =1.50 in and i=8.00 in> 

both ends 2 ure ideally clanged. ® 


b: Fragment is_a 1-inch diameter steel sphere of nominal mass 

0.3815 X 10 (lb-sec )/in; midspam Dn pact. 

Nominal Mechanical Properties of 6061-T651 Aluminum Beam Material 


a 

• o 

a 

u 

' "m- ./ 

P ' 


P ^ V 

(ksi) 

(ksi) 

(in/in) 

(psi) 

* ih/in^ 

42 

45 

. 0.20 

lOxlO®" 

^0.098 '! 0.3 

• . * #c- • 


*Ref. 16 subscript; "o" refers to static proper^es. 






TABLE 10 (Concluded) 



159 





STRUCTURAL MODELS 




DIAMETER SOLID STEEL SPHERE 
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FIG. 15 POST-TEST PHOTOGRAPH OF STEEL-SPHERE-IMPACTED 6061-T651 ALUMINUM 
BEAM CB-16; IMPACT VELOCITY = 2870 IN/SEC 
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43 BLEM, EI.-SII BIT A 
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43 ELBH 







DISK MATERIAL : A«286 

BLADE MATERIAL: S^-IS 


(a) Schematic of Rotor Modified for a 
Three-Fragment Burst 


FIG. 22 GEOMETRY AND VIEWS OP'tSB TURBINE ROTOR TRI-HDB-BORST FRAGMENTS 
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(c) Photograph of Post-Test T58 Turbine Rotor 
Tri-Hub Bxirst Fragment 


FIG. 22 CONCLUDED 
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FIG. 



(a) Pre-Test 



(b) After Rotor Burst but Before Impact 

23 DYNAFAX PHOTOGRAPHS OF T58 TRI-HUB ROTOR-BURST ATTACX AGAINST A 
STEEL CONTAINMENT RING IN NAPTC TEST 201 
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(c) Time After Initial Isqpact, TAII-200 ysec 
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(i') 


TAII=1J58 usoc 


FIG. 23 CONGI.UPKD 
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TAPE CH. #1 
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#37 (180*) 
SCOPE CH. #4 
TAPE CH. #12 
€ PERM - -33093^ 
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TAPE CH. #7 
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SCOPE CH. #3 
TAPE CH. #11 


#1 ( 0 *) 
SCOPE CH. #1 
TAPE CH. #9 


#5 (20*) 
TAPE CH. #3 


#21 ( 100 *) 
TAPE CH. 


#17 (80*) 
TAPE CH. #5 


SCOPE CH. #2 
TAPE CH. #10 
€ PERM - -32266-C<€ 


#9 (40*) 
tape CH. #4 
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FIG. 25 POST-TEST 201 CONFIGURATION AND LOCATIONS OF STRAIN GAGES 
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LOAD OP 7000 POUNDS WITH A 150 MICROSECOND DURATION 
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tip 
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jT FRAGMEUf WITH IDEALIZED FRAGMENT-SIZE CHOICES 










FIG. 28 GEOMETRIC, TEST, AND MODELING DATA FOR, THE 4130 Sl^BL CONTAINMENT RING SUUJECTEO TO 
TRI-HUB T5B ROTOR BURST IN NAPTC TEST 201 »w»«rr Mnu ouiMEin-fcD TO 



fig. 29 THE EF FECT S OF IDEALIZED FRAffifENT SIZE ON THE PREDICTED 
extreme response of the steel containment ring IN NAPTC 
TEST 201 
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FIG. 34 l>R£OIC*rBO AND MEASURED KING OUTER-SURFACE PERMANENT STRAINS ON THE NAPTC TEST 201 RING 



(b) Cyllodrical Shell 



(c) Curved Shell 
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rfirfVr, 




SWB X 10 (IN- LB) /LB 



FIG- 39 
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(a) Small Rotor 2-Fragment Containment Post 


Test Results 



FIG. 40 


(b) Small Rotor 3-Fragment Containment Post Test Results 


ILLUSTRATION OF POST-TEST CONTAINMENT RING AND 
FOR T58 2-FRAGMENT AND 3-FRAGMENT ROTOR BURSTS 


fragment CONFIGURATIONS 
(REF. 37) 
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SPECIFIC TRANSLATIOllAL FRAGMENT ENERGY, STFE X 10 , (IN-LB) /LB 



FIG. 41 


NAPTC IMPACT TESTS OF THREE EQUAL-SIZE T58 TURBINE ROTOR 
BLADECHDISK FRAGMENTS AGAINST 4130 CAST STEEL .CONTAlNrlENT 
RINGS AS A FUNCTION OF RING AXIAL LENGTH (WIDTH) 
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RING MASS 

TOTAL FRAGMENT MASS 



FIG. 42 NAPTC T58 TURBINE ROTOR BURST IMPACT DATA OF RING-TO-FRAOffiSIT 
MASS RATIO VS. (KS)ot/°C NUMBERS OF EQUAL-SIZE 

BLADED-DISK FRAGMENTS IMPACTING 4130 CAST STEEL CONTAINMENT RINGS 
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WITH THE CB-18 IMPA(ri’ CONDITKXIS 
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, (b) History of Fragment Energy Loss Ratio 
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t !(.. fil) (.'pMI'AHISON OF CIVM-.IET 4U PKEUICTEO UFFEK-SURPACB HIOSFAN TRANSIENT STRAIN FOR THE CB-18 
UEAM WI’III THAT F<JM A SCAIJ:;i) UEAM SUUJECTEO W MIUSFAN IHPACI* AT THE CB-18 VELOCITY BY A 
2-p FRAGMENT OF 75 FEM CENT GREATER MASS 
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NOTE: a IS POSITIVE CLOCKWISE 

6 IS POSITIVE C30NTERCL0CKWISE 



FIG. 62 SCHEMATIC OF A SIJIPLIFISD CANTILEVSR-BEA:'! DEFLECTOR WITH 
t:4PACT AN IDEALIZED FRAGMENT 
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(b) Dofloctor ProfiX® «nd Praymint Location 
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(a) Transient Strain at Station x«0 


FIG. 65 MEASUREMERTS AHD/OR gREDICT IONS OF TRAHSIEIfr LO HGI T UD IHAL 

UGRAaGIAH (GREEH) STRAIN OH THE SURFACE FOR VARIOUS SPANWXSE 
STATIONS OP EXPLOSIVELY-IMPULSED 6061-T6S1 ALCMINDM SEAN C3-4 
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FIG. 63 CP M T INO EO (C3-4) 
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(c) 0!pp«t-Surfac« Translene Strain at x*2.20 in 
PIG. 65 OQBTIHDa) (CB-^) 
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FIG. 65 CONTEIDED (CB-4) 
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TIIE FINITE-STHAIN ttJHMUlATION 
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if Id. 6 / CONTINUBO (CB-16) 
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ITINUEO (CB-18) 
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TIM^l AFTER INITIAL IMPACT (USEC) 
(n) Station x«4.00 in, Lower Surface 
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AFPENDIX A 


REVIEW GF GO V E RM IR G EQUATICHIS 

A^l Introdaction 

For the reader's convenience, a brief reHew is given here of the basic 
equations xxpoa idiich the ccoqnxter codes discussed in this report for predict- 
ing the large-deflection, elastic-plastic, transient structural 2-d responses 
(of beans and/or rings) to externally-applied loads are based. In particular, 
the equations' of notion based upon the Principle of Virtual Wor)c are described 
for a structure represented as a general solid continuum* These equations 
are then utilized to provide a finite-elem en t analysis of 2-d bean/ring struc- 
tures; are the pertinent str 2 dn-displacement relations, the selections of 

appropriate generalized displacements and assumed displacement fields for the 
elenent, and the derivation of elanent properties* A convenient stress-strain 
representation for the mechanical behavior of the bean/ring structural material 
is described. Finally, an illustrative transient-response solution procedure, 
utilizing the timewise central-difference (explicit) operator for a 2-d 
structure sxabjected to prescribed externally-applied transient loading is 
given; Refs. 3 and/or 5 may be consulted for a description of the use of 
the (ijiq)licit) Boubolt timewise finite-difference operator for solving the 
nonlinear equations of motion* 

TO avoid unnecessary repetition, the procedures used to anadyze the 
responses of 2-d beams ani/or rings to fragment ijjq>act is not included* For 
that information, the reader is invited to consult (1) Ref. 4 concerning the 
CTVM-JET 4B (single-layer structure) program or (2) Ref. 5 for the CIVM^-JET 5B 
(multilayer structure) program. 

A. 2 Formulation of the Equations of Motion 

Consider a deformed continuum (or structure) to be in equ i lib r ium under 
the action of body forces and prescribed externally-applied surface forces 
(or stresses) on the surface area A^ ; let prescribed geometric boundary 
conditions be given on the remainder^f the surface area * Let this body 
undergo arbitrairy infinitesimal virtual displacement 6u^ consistent with the 
prescribed geometric boundary conditions; for this situation, the Principle 
of Virtual Work requires that [17,20, and/or page 248 of Ref. 48): 
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/ -r^ 

X '■* J fa..rfA = 0 a.i, 

.Trl::i:*:r;x'rr ~ « 

- x,‘ a •“•’ 

area of the uadeformed state) y i (“easured per unit 

the u are the displacem^t cop^^nents p 

i. . ,iv« ^ .«r «„rLL r 

surface area a — ^o *“<* the (Reformed 


as 


a. a^a. ^ 


^ ? \^L%* . ~l~ CL //.* ^ 

i'f f,£. ^,£ 


(A. 2) 


»*ere ( ) . denotes covariant differentiation with respect to th 
coordinate Mai 2 n • *^espect to the material 

e t (3-l,2,3) using the metric tensor of the-undei^ 

By e»5,loying the concept of D*Alemberfs Princiole 
PoB^ may be regarded as consisting of the D-Alembert i 

and other body forces p f^ ^ ) h 

00 ^ eta.). Hence, one may write: 

A 6 / = -/ a*’ + f {*' 

® *'« Jo (A.3) 

«»r. O <aota dc«,a «a ^ 

Applying Bq. A.3. Eq. A.l becomes 


V 
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As usual la the assumed-displacemeit version of the finite elenent method, 
the actual continuous structure is modeled as a con^tibly- joined assonblage 
of finite elements. E3q>ressing Bq. A.4 in terms of the contributions from 
each of the H finite elements or regions, one obtains: 


N 

If = o 

Next, for each finite elem^t one chooses a selected finite nun*er of control 
stations or nodes at %*ich ^control generallred displacements" {q} are defined; 
these "control stations" mAy be on the bound;tries and/or the interior of the 
elenent. Osing these control-station q's, one assumes the spatial distribu- 
tion of the displacements in the element to be expressible in of 

appropriate interpolation functions such that one may wite 




where i (i-1, 2, and/or 3) represents a selected direction. Since the strains 
are e^quressed in terms of displacement gradients and combinations thereof 
(see Bq. A. 2)# one nay write 


Renee# 


Hert, expressing Eq. A.S in terms of the elenent* s generalized displacaaents 
and using Eqs. A.6 and A.7, one has 
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.h«. th. fbUoid., for «ua fi.it. .i...,, ,M,.a.ipt -.■) = 

W= J 


(A. 9a) 




f»>}=I /Djs‘^otv 




(A. 9b) 


w=lfU[i^Js‘W 


•A' 'b 


(A. 9c) 


(A.9d) 


throughout the volume 


Note that {p} and [h) involve stress information 
of each finite element* 

Next, it is desired to egress Eq. A.9 in terms of the independent 
global generalized displacements {g*} for the cooq.lete assembled discretized 
structure rather than in terms of the {g} for the individual finite elements 
or every fxnrte element one can taJce this into account and express the (g) 
of each element in terms of the {g*} by 




(A. 10) 


Apply^ Eg. A.IO to Eg. A.9 to describe the system in terms of independent 

global generalized displacements {g*}, one obtains: 

N 


I Uf\(L-ur} * fw-jff*} -{f*}) = 


(A. 11) 


where 


- [Jj [nt][j] 


(A. 11a) 
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(A.Ub) 


{!>•} = [J]>} 


K} = pfff} 

Performing the summation and since the 6q* are arbitrary and independent, 
one obtains; 




- where 








(A. 12) 


(A.I2a) 


(A. 12b) 


(A. 12c) 


(A«12d) 
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Hote that (Ml is the global mass natrix, [H] {q*} represents "forces" 
associated with large deflections (nonlinear terms in the stara< n - <<4 «pi 
relations) and plastic effects*, {p} represents intenua forces,, 

plastic forces, and some plasticity effects associated with the 
terms of the strain-displacement relations*, and {p} represents the 
generalized nodal load vector accounting for externally-applied distributed 
or concentrated loads. 

The equations of motion represented by Eq. A. 12 are usually the 

"unconventional form" but are very convenient for analysis and cooqiuter 
implenentation. On the other hand, one can obtain the "conventional form" 
of the equations of motion by eoq>loying the stress-strain and strain- 
displacanent relations to express the stress in terms of the element 
generalized displacements {q}, using the concept of initial strain: 



(A. 13) 


where yf, is the plastic stradn component of the total strain^ Y and 
ij)cZ « 

E is the matrix of elastic constants. Applying Eqs. A.13 and A. 10 to 

Eqs. A. 9, A. 9b, and A.9c, one obtauns: 


[M ]{ r} -[K]{ fl= [ f] *{ F/1 «... 

where [M] is the global mass matrix, (K) is the usual global stiffness matrix 

(for linear-elastic small deflection behavior) , {p} is the generalized 

load vector representing externally-applied distributed or concentrated loads, 

} represents a "generalized loads? vector arising from large deflections 

and is a function of quadratic and cubic displacaaent terms, {p^} and {p^} 

P P 

+ 

See page 50 of Ref. 17. 
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are the generalized load, vectors due to the presence of plastic strains and 
are associated, respectively, with the linear and Mn^ear terns of the 
stradji-displacenent relations. 

Equation A. 14 represents the ^called "conventional" formulation of the 
equations of notion for nonlinear responses, where the large deflection and 
plastic effects are taken into account through the use of generalized (or 
pseudo) loads which are function of plastic strains and displacements. The 
element stiffness matrix (kj may be readily shown to be 

f^] = J 

AlSOr 

[4*] = [ jf[^][j] 


and 



(A. 14c) 




A more detailed discussion may be found in Refs. 17 and 20. 

Note that the very cocq>act and efficient "unconventional form", Eq. A. 12, 


of the equations of motion is i^ed together with the timewise central-difference 
operator to predict 2-D transient structural responses in the JET 3A and JET 3C 
codes of Ref. 3 and in the CIVM-JET 4B code of Ref. 4. Or. the other hand, in 
an endeavor to enploy a larger allowable time-step-size increnent dt for the 
transient response solution, the Houbolt timewise finite-difference operator 
is used in the JET 3B and JET 3D codes of Ref. 3 and in the JET 5A and 
CIVM-JET 5B codes of Ref. 5; these codes utilize the "conventional form", 

Eq. A. 14, of the equations of notion. 
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A.3 Finite Blanent An alysis of Curved Beams 

consider an Inltlally-nndefomed, arbitrarily-curved, variable-thickness, 
single-layer bean or ring subjected to prescribed transient extemally-appUed 
surface loads and to only D'Aleabert body forces (inertia loads). Let it be 
assuoed that the ring consists of ductUe netal and that a large-deflection, 
elastic-plastic transient response will be produced. For analysis the 
structure will be represented by a coa?)atibly-joined assemblage of H finite 
elements, one of which is depicted in -Pig. A.l where its geometry and 
nomenclatnre are shown and where the deformation plane is n, C .• the 
coordinates n along and Z normal to the centroidal reference axis of the 
beam are en 5 >loyed as the reference coordinates for this curved beam element. 

It is useful and convenient to use the foUowing geometry to describe 
this typxcal curved beam element and to approximate the actual given cong)lete 
beam or ring by a finite number of these "typical elements". Note first 
that a global V,z Cartesian reference axis system as *»ell as a local y,z 
cartesian reference axis system are defined; for the latter, the +y axis 
passes through the ends (that is, nodes i and i+1) of the element and makes 
an angle +o (for this ith element) with the +Y axis. The slope, ♦, of the 
reference circumferential axis n, which is the angle between the tangent 
vector i to n and the y-axis of the local-reference Cartesian_frame may be 
approximated by a second degree polynomial in rj, as follow (17J : 

‘P(’l) = b. + b, ’r + 

where the constants b^, b^, and b^ can be determined from the geometry of the 
curved beam element as described next. Assume that the chance in element 
slope ^ between nodes i and i+1 is small such that 

= I (A.16a) 

and 

^i) = ^i4l ” (A.16b) 
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restricts the slope change within an elsjent to <15 degrees. The arc 
length, n., of element i is approximated to be the same as the length of a 
circular «c passing through the nodal points at the slopes and 


hence, is given by 


^ y 


(A.17) 


Sin 


(- 


where is the 

by 


length of the chord joining nodes i and i+1, and is given 


The three 


l:=[(z - zj% (y -r)']' 

constants in Eq. A. 15 are then determined from the relations 


(A. 17a) 


<p(o) = <P, 

(A. 18) 

f =J 

o ^ 

From Eq. A. 18. the constants in Eq. A. 15 are found to be 

t = 4>i 

* . / (A. 19] 

b, = 

KcordiMly. th, radiu, of o«v.tur., R. .t a,. c«.troid.l aai. t. 

R - -(3d/3n)‘* - a.. coordiMtas yin) a.d 

Z(n) of the centroidal axis are given by 

y<;,) = Y. + J -=-*[ + “"I 
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and 


where 


'L 

Zc?i) = Z. +J + 



(A. 20b) 


(A. 20c) 


The thickness variation along the elonent is approximated as >v»-iT»g 
linear in n between nodes; thus. 


*f>t) = hji 



(Ae21) 


This CQB5>letes the needed description of the geoaetry of the curved 
beam element. To be revie%ied next aure the strain-displacenent relations, 
the assxamed-displacement field for the elenent, and the resulting "eluent 
property matrices" identified in Eqs. A.9a-A.9c and in Eq. A. 14a. 


A. 3 . 1 Strain-Displacement Relations 

Let it be assumed that the curved beam d eform^ af? gnr^ j ng to the 
Bemoulli-Euler hypothesis; that is, (a) plane sections remain plane, (b) 
nor mals to the midsurface (or reference sxirface. —• centroidal axis) before 
deformation remain normal to that reference surface after deformation, auid 
(c) these nor m a l s do not stretch. Thus, the displacements v(n,C) and w(n,^) 
at any location within the element can be expressed in terms of the 

displacements v(n) and v(n) at the reference axis (%diere C=0) by 


where 








5 W V 


(A. 22) 


(A. 22a) 
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By applying Bq. A.2 to Eq. A. 22, one can sho^ that the only Lagrangian (Green) 
straon component of iaq«rtance in this problem y (hereinafter denoted .is 


I for sia^licity) is given by (20] ; 

2 




+ ? 


Linear 


Nonlinear 


f-ifi 


(A. 23) 


Linear 


Membrcuxe 


or more conqpactly as 


Bending 


Ki - '5 


%/here 


% = 


V' = 


w 




V 

r" 




/C(^)— - 




(A. 23a) 

(A. 23b) 
(A. 23c) 
(A. 23d) 
(A.23e) 


Note that the part of which does not depend upon C is often termed the 
membrane strain, and is co=n>lete in the sense that all of the linear and the 
nonlinear terms of this portion are present in Eq. a. 23. On the other hand, 
only linear terms have been retained (in*Refs. 3, 4, and 5 and for present 
purposes) in the bending part of this strain; the inclusion of nonlinear 
bending terms is discussed, for exanq)le, in Ref. 20. 

Sine, (1) th. displac«Mt, v(n,SI and l(n,C) can b. « 5 cesa*i in terms 
of . and » and 12 ) tbs ctansional strain is erpressad =ontani«,tlv in 
tamss Of th. ^tiuas X and it is both cnvmUent «si tmtutai to msplox 
those four ^tities as ■,m.«:alitmi displarmnents- {,) at each end. i and 
1+1, Of the ith finite element. 
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A«3,2 Assqned Displacement Field 

To ocooaot tor tt. diopUo«ont. aod ...u ri,ld-tod, 

tt. di.plac=»ot field toter-tto fom , 17 , 101 . 





ft 


cas<^ 

-Sin ^ 

-sift tp 

cas<p 

^ O 

O 


O O o 


This can be written more con^actly 



as 


H=- 


V 

w 


G.C’t) 

L 




The eight undefined parameters fi a v 

. w "1 ®8 ej^ressed in te ts of th« 

.i,^t oeofolieed di,,u„o„fe J,, „ fi..d„ 


elanent — four at each end (or node) . 


f, 


These {q} are defined by 


tsJ 

Hence one may %n:ite 


= 1"‘ f «. 


25 ) 


where 




(A. 26) 
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[A]= 


C*s4i 


O 

O 

o 

o o 

a 

•sin 

cas^. 

O 

o 

O 

o o 

o 

O 

o 

1 

Q 

o 

o o 

o 

o 


o 

1 

o 

o o 

o 


SinS 

A„ 


o 


III 


casi 

^3 

\ O 

- z 
Ti 

O 

o 


O 

1 

^.(*1 




o 

P 

O 

1 

-i(*\ ^7, 



(A. 26a) 


(A. 26b) 


= (X.: X) (tH 

Correspondino to the assumed displacement field Eq. A. 24, one finds 

^ = lo c I ^7 37 .£]{?} = Ig,J{^} (».27a, 


and 


X = [o . O ... 


27b) 


Hence, the quantities C^(n) and ic(n) appearing in the strain 
^rfiich is denoted hereinaifter sin 5 >ly as e are given by Eq. A. 23a and become. 
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(A. 28b) 


OoQbining Eqs. A. 24 through A. 28b, one obtains 
and 

where 

I J[(J][A ‘ ] for i - 1, 2. 3 


(A. 29) 


(A. 30) 


(A. 30a) 


and 


^ ‘ -^v 

IBxJ[u]-[o o 1 


- ^ ,?:f 5 ;^ J 

-2 -3^V-^VJ 


In the process of solution, it is necessary to evaluate the strain 
iaeroKiic from to tii» t^. osi.^ e,,. ,.i 2 ood ».14, on. has 


where 


Also, 


from Eqs. A. 24, A. 23c, 


A* 27a, and A. 29, one ma^ write 
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|:i 


4 

U(^) 


I 

|yj 

i 


Ooo 
'T J 


r ]f?}=[N(,)]f f ] 


(A. 32) 


Accordingly, consistent with the assumed displacement field, one may express 
the -velocities v, w, and ^ (vhere ( ’ ) denotes the time derivative) as 





(A. 33) 


A. 3>3 Finite Element Properties 

The finite element property matrices of interest are (ml* {p}, [hj, and 
{f} given, respectively, by Eqs. A.9a, A.9b, A.9c, and A.9d for use in the 
"unconventional formulation" ~ to vhidx the present review is being 
restricted. In addition, however, as explained later it is useful to 
evaluate the el«nent stiffness matrix Del for linear behavior as given by 
Eq. A. 14a. 


Element Mass Matrix [m] 

Since 

[-]= J 

one needs to form 0^(5^) as indicated by Eq. A. 6. For this case, using 
Eqs. A. 22 and A. 32, one obtains 



1 O -5] 

r 

V 

0 

0 

w 

rj 


Comparing Eq. A. 34 with Eq. A. 6, 



(A. 34) 


[«i - [ : : 
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Htoce, fzan Bq. a. 9a, {m] becomes 






(A. 36) 


« «y •t.tta. n Md tl» i 

for ttln rings, one obtains R 


f"*] “ J f/V/c>t)] [Bcy)] [N<x)] 

^ 8X3 91fa 9 va 


fXJ 3X3 
for a beam of uoiform width b. 


(A. 36a) 


3X8 




I O 
O I 
^^0 




h(^ O 

O hex) 

O o 


(A. 36b) 


fe% ) 

/2 


toat IS the mass per unit volume of the undeformed body and the 
^teqratron rs performed over the undeformed volume, m practi J, the 
^.gxation ia B,. ».36a to fot. i, option, _ w 

^ “S' the eseo«d diepUcel 

c^^ed the , 3 = 0 ^ .ee, .ettth to , h..,Up“ pouted „ttU „oe 

ta SoheeT^r^ ««cee u dUcoaeed 

Element Matrices {p} and pil 

Note that these matrices are given by 


(t)-[yD.Js‘'oo 

i‘j=J fyic;js“a>/ 


(A. 9b) 


(A. 9c) 
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Since the onl^"aatin 9 non-zero strains and stresses" present in this problem 
® indicated sionuttion becomes quite con^ct; for 

convenience, let 5 a and = e. Also, note from Eq. A. 7 that D 
pertains to the linear part of y^^;; hence, from Bqs. A. 7, A.23a, and A^M, 
it is seen that 




Alsop the no h l i n ear-* term contributions provide 


(A- 37b) 


Bence, Eqs. A. 9b and A. 9c become 

ft"}"! •*-“> 

V 

[<-]= J 

Since D^, D^, and 0^ are functions only of n, one can perform first the dC dC 
integration and then the dn integration. Accordingly, let 

and (A.40) 


Applying Eqs. A.40 to Eqs. A. 38 and A. 39, one obtains, respectively: 

T.- ^ 

{ f*} = I LfT) fli] Mc^)^ 

-[*>] °1 <^’1 

C 


Since the Stress a varies, in general, at each time step or instant at every 
location within the volume of the element, {p} and [h] are evaluated by 
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Of a. u-oio. ,of^ 

process* 

Element Applied-Force Vector ff> 

Ihi, Ol^t .ppXiod-foro. «ctor {£) 

t»r. oioce. for m. aiscuMioo. u h« bo„ tlat oo tody forcto f^ 

are present. Hence, Eq. A.9d sljt5»li£ies to 


ff} = i 1/ dA 

^nl 


(A.41) 


For convenience and simplicity, let it be assumed that the surface inteqra- 
txo^ Of the prescribed appU^ surface tractions have already been 

'I - - ^ provide force resultants 
(n) and P^(n), and the moment resultant M (n) each per unit spanwise 
distance _ 


In this case, Eq. A.41 reduces to 

{^} = I [n<’d] ji( I 


(A. 42) 


Elemen t Stiffness Matrijt ficl 

to tto prootot discutoion, tto toootioos of totioo for tto ecptot. 

cture are based on an "unconventional" formulation in which 
the c™nal elastic stiffness matria, tK), for the "complete assembled 
st^c^e does not appear e^^licitly. However, in order to calculate an 
^lowable tome step sire, At, for the conditionally-stable central-difference 
trmewrse operator, the largest natural frequency contained in the (linear) 

todot Of to, ^ ^ ^ J 

^colouon. to. .Itotic toiffotos totrja f„ toe to,»oito stroctor, ,„st 
be cooputto. Th. .leetic etiffo.to totri, ftl for «. 

Eq. A. 14a: ^ 




(A. 14a) 
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For the present problem Is singly the elesttc modolus B and 

[Dj^J = ([dJ + Hence, one obtains, after performing the dC dc 

integration for ^htn rings: 

A. 4 Stress-Strain Description 

For convenieace and sia^licl^, the mechanical sublayer model is 
employed to describe the mec h a n ica l behavior of the structural 
»*lch is assumed to be ductile. In the mechanical sublayer model, tre 
u niaxial tension (or coinpression) stress-strain curve of the material is 
first approximated by (n+1) piecewise-lihear segments which are 
at coordinates ejj)» k * 1, 2, ...nl, as depicted in Fig. A. 2a. Next, 

the material is envisioned as consisting, at any point in the material, of 
n equally-strained "sublayers" of elastic perfectly-plastic material, with 
each sublayer k having the same elastic modulus E, but an appropriately 
different yield stress (see Pig. A. 2b). For example, the yield stress 
of the kth sublayer is 

^Jt ~ (A. 44) 


Then, the stress value, Oj^, associated with the kth sublayer can be defined 
uniquely by the strain history and the value of strain and strain-rate 
present at that point. Taken collectively with an appropriate weighting 
factor for each sublayer, the stress, a, at that point corresponding 
to strain e may be expressed as 


n 



where the weighting factor for the kth sublayer may readily be confirmed 


to be 




(A. 46) 
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%^ere 






(A. 46a) 


The elastic perfectly-plastic and linear strain-hardening constitutive 
-relation -may be treated as special cases. In the case of elastic perfectly- 
plastic behaviorr there is only one sublayer, and in the case of linear 
strain-hardening material there are two sublayers and the yield limit of the 
second siiblayer is taXeh sufficiently high so that the deformation in that 
sublayer remains elastic. Finally, it should be noted that the above rules 
for the mechanic 2 d.-sublayer model require that the stress-strain curve being 
represented must be such that the stress either increases monotonically with 
increasing strain and/or reach a limit but it can not decrease? some of the 
isq>lications of this restriction are disctissed in Section 6* 

From the cooqnitational point of view, the use of the mec h a ni cal -sublayer 
model is very convenient to analyze problens with general loading paths 
including loading, unloading, reloading, and cyclic loading. Its features 
include the "kinonatic hardening rule" lAich takes the Bauschinger effect 
into account (see Fig. A. 2c). Also, this mechanical sublayer model may 
readily acccmniodate the strain-rate effect. Figure A. 3a illustrates 
schenatically the uniaxial stress-strain behavior for a str 2 tin-rate dependent, 
el 2 istic, perfectly-plastic material whose rate depezKience is described by 
[24,25]: I 



while Fig. A. 3b depicts the corresponding behavior for a strain-hardening 
material which is represented by the mechanical siiblayer model, each stiblayer 
of which has the same values for the strain-rate constants D and p. For 
this special type of rate-dependent strain-hardening material, the stress- 
strain curve at a given strain rate e is sinq>ly a constant magnification 
of the static stress-strain curves along rays emanating frcm the origin. 
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Let it be assumed that at time instant t^ the stress and strai ns of all 
s^oblayers and# hence# of that materxal point are known* Betoreen t^ and 
let there be a strain incranent ^£^ 2 ^ ^ associated strain rate (Ae^j^)/(At) . 

The individual sxiblayer stresses are found conveniently from the following: 


A€ 


> o 
= o 

< o 


- 


^ 1 . 




= CT 


% 




t 




(A. 48) 






=: - <r. 




%rtiere a "trial" value of stress at time t^j^ in the kth sublayer defined by: 
t 


(T. = + E 


(A* 49) 


and 


= ^ve ( I + 


('*hH ) 


(A. 50) 


is analogous to Eq. A. 47 and applies to the kth sublayer. Once one has 
determined for all sublayers at time t^j^, Eq. A. 45 is used to compute 

the actual uniaxial stress at that material location at 

A. 5 Treuisient Response Solution Procedure 

For present illustrative purposes, let it be assumed that Eq. A. 12, 

the "unconventional form" of the equations of motion, is to be used to carry 

0 

out a timewise step-by-step solution in small increments At in time by using 
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the timewise central-difference operator. At time t=0, the structure is in 
a known unstressed unstrained state with {q*} = { 0 } and {q*} = {a}, for 
example. Also, the structure is subjected only to prescribed extemally- 
applied tremsient loading. 

Accordingly, Eq. A. 12 is to be solved at a sequence of instants in 
At apart by eiqjloying the following expUcit, conditionally-s table, central 
difference, finite-difference time operator approximation for the accelera- 
tion ^ at time t^; 

" C4t/ 


Where 0(At) means that this finite-difference approximation has a truncation 
error of order (At) . Also, one nay approximate the veloci^ ^ at t-im. t 
by: “ 

• 2 , 

fn. = J7At) * 0(At) 

(A.Slb) 


At time i n sta n t t^ the equations of motion (Eq. A. 12) become; 

[M]{rl^fPh[H]jfi={Fi 


(A. 52) 


In Eq. A. 52 all quantities, in general, except [M] change with time. If the 
solution of Eq. A. 52 has been obtained for eeurlier times, one conq>ute 
from Eq. A. 52 and then obtaOn from Eq. A. 51a. 

Assuming that at t=0 the structure is in a knom condiUon {q*} =0 and 
{q*}^j={a}, for example, one can readUy obtain {q*}^^ at t = ndt for ^1 from 
the following Taylor series eacpansion: 

{ ?*}, = { ?'}. + { fi +7 / f'l Oatf 

since is prescribed and all other quantities are knom. 
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In the timewise stepHby*step solution process involving lerge-def lection 

elastic-plastic transient responses, {p} and [H] change with time and hence 

m m 

must be reevaluated, in general, at each instant in time. These quantities in 

turn are composed by assembling the contributions {p*} and [h*] from each 

m m 

finite element. It is seen that these quantities involve voliaae integrals-o£ 
information- involving, in general, the stress state In practice, these 

evaluations are carried out by appropriate numerical integration; namely, 
Gaussian quadrature. 

At any instant in time t (m»0,l,2, , one needs to solve Eq. A. 52 for 

m 

{q*} , which is of the form: 
m 

[m] f X(i)} — (A.54) 

%diere 

[M] is a known banded positive definite synnetric matrix (the 

mass matrix for the restrained or unrestrauned structure, 
%diichever case is being treated) . 

{x(t)} is a vector of unknowns %diich must be determined by solving 

m 

Eq. A.54. 

(b(t)} is a kno%m vector (representing all terns except [M] {q*} 

m m 

in Eq. A. 52) . 

In principle, one can aJ.ways form the inverse matrix [Mj ^ and pre-multiply 
Eq. A.54 by (Mj*^ to obtain 

which results in the solution: 

since [M] ^(Ml = [I] %rtiere [I] is the unit diagonal matrix. However, it has 
been found that independent of the number of time instants at which one ^ 
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wishes to solve Eq. A. 54, such a procedure is not as efficient as is the 
CholesJci method [49]. 

Briefly, the Choleski method involves factoring the matrix [Ml to form 
a lower triangular matrix iLl and an upper triangular matrix (which is the 
transpose of the former) such that (MJ = (t] (Ll’’ where (Ll’^ is the transpose 
of ILl. Thus# Eq. A. 54 may be reiiritten as 

IL llL f [xCijj == { h(t)l (A.56) 

Next, form an intermediate matrix {y}^ »*ich is defined as 

{ “ f LJ (A.57) 

From Eqs. A. 56 and A.57, it follows that 

[L]{yI= {ut>l 

At each time instant, one solves Eq. A. 58 for {y}^ very readily because (LJ 
is a lower triangular matrix. One then solves Eq. a. 57 for {x}^ very rapidly 
also by algebrcLic back-substitution. 

The following gives a concise step-by-step description of the typical 
pj^ohlmo formulation and solution process. 

from a set of given initial conditions at ^ilno t = t = 0 on 
the generalized displacements ({q*} = {o}, for example) and the generalized 
velocities {q}^. one can solve Eq. A. 52 for {q*}^ at time t^ and then employ 
Eq. A. 53 to compute {q*}j_. A slightly different but simUan procedure is 
then used to advance the solution in successive time increments At. The 
process involved in using the finite-element method and the present timewise 
solution procedure follows [20J : 

Ste£j,: Construct the mass matrix (mj for each finite element and then 

assemble these contributions according to Eq. A. 12a to form the mass 
matrix (M) for the cai 5 )lete assanbled discretized structure. This [Ml 
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represents the "final" mass matrix if the structure has none of its 
generalized displacements constrained (that is, held equal to zero, 
for exan?>le); however, if such constraints exist, one forms a reduced 
or fonstraihed mass matrix (and, in fact, a reduced set of the equations 
of Motion) by deleting the rows and columns of [M] associated with those 
generalized displacements which are prescribed to be zero- Next, this 
conttraipcu mass matrix is factored to consist of a lower triangular 

matrix [L] and an upper triangular matrix [L]^ according to the Choleski 
scheme: 

[Mj=fL][Lf 

Since [Ml does not change in value with time as the transient structural 
response proceeds, one needs to determine [L] and [L]’^ only once — these 
quantities need not be re-evaluated at each time step of the calculation. 

SteE^: The prescribed extemally-appUed transient forces can be employed 

to calculate the generalized applied forces {f*} acting on each discrete 
element at each time instant t^ of interest. These, in turn, can be 
assembled according to Eq. A.12d to form the assembled applied- loads 
vector {f} for the coaq>lete assembled discretized structure. 

Ste£_3 ; Assuming that at zero time (t = 0) , the generalized displacements 
{q*>o » 0, the generalized velocities are nonzero (4*} =. {a}, and that 
nonzero external forces {f}^ are present, Eq. A. 52 becomes 


Min ={F\ 
” fL][LnrHF}, 


(A. 60) 


(A. 60a) 


from which one can calculate {q*}^ by using the earlier-described 
Choleski scheme. Then from Eq. A. 53 one obtains 
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(A. 61} 




where 


hi = { fi - { fl 

{ fi = H = prescribed initial geaeratlized 


velocities 


Also, 


{ fi = {f} fl 


(A. 62) 


(A. 62a) 


(A. 63) 


For this case, however, it has been assumed that {q*} = {o}. Thus, 

the displacanent configuration {q*} at time t = t + At is known. 

-1- 1 o 

step 4 : Knowing the generalized nodal displacement incranents 

{Aq*} = {q*}j^ - {q*}^ and the generalized nodal displacements 

{q*}j^ at time tj^, one knows 2l1so the unstaurxed individual elanent 

quantities {Aq}^^ and {q}^^ via Eq. A. 10. Hence, one can calculate, 

in general, the strain incranent (AY-.), developed from t to t at 

i-3 o 1 

every Gaussian station (or point) required over and depthwise through 
each finite element from Eqs. A. 23a cuid A. 30: 





(A. 64) 


With a knowledge of (a) the stresses at t^ = tj^ - At, and (b) the 

strain increment (Ay..)., one can determine the stress increments 
It ^ ii 

(AS )j^ and the stresses (S -’)j^ at time at each Gaussian station 
by using the pertinent elastic-plastic stress-strain relations. 
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incliiding the yield cxjndiUon and floe rele, in general. However, 
a simpler procedure will suffice for the present probleas*. 

Steg_S: Next, one can calculate {p}^ and [h]^ for each individual finite 

element by using Eqs. A.38a and A.38b. respectively. Assembly of this 
formation according to Bgs. A.12b and A.lfo, respectively, provides 
{P}j^ and Since the prescribed generalized force vector {p} is 

available from known {f}^ information, the equation of motion,* Bq^ a. 52, 
at time insta n t tj^ becomes: 

In the interest of minimizing con^ruter storage and the number of 
manipulations, one first forms for each individual element 

^^n^l ~ ~ ~ Then one forms the right-hand 

side vector of Eq. A. 65 fay 

[J ff b, b, ... b I 

*• J ( -*2 (A.66) 

For clarity of discussion, however, the form of the equation repre- 
sented by Eq. A. 65 is used here. 


SteE_6= Si«=e th. ri,ht-ta»d sid. o, A. 65 i. no. c c«. 

tte Chol.su scha» to sols. th. follosio, .qoatioi. for th. .coU«..ioo 


{q*}. 




(A. 66a] 


Ste£_7: With {q*} now known, one can calculate the generalized displacement 

increnent {Aq*}^ from Eq. A. 51a as 


where 




(A. 66b) 


+ "" 

See Eqs. A. 48 - A. 50. 
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(A. 67a) 



(A. 67b) 


■ li'l - {fl 


Thus, from Eq. A. 67a one has 




The process then proceeds cyclically from Step 4 onwards for 
time steps as desired. 


(4.37) 
as many 


py^r the central^ifference tbne operator applied to a systaa of equations 
such as Eq. A. 12 or Eq. A. 14, it has been shown [50,51] that At most be less 
than or equal to 2/a,^ to avoid esponenUal growth of error [roumloff , 
truncation, gross) where is the largest frequency embedded in the 
mathematical model of the system. The criterion At < 2/u must be 
satisfied for a linear dynamic system such as that involv^^ small^lsplace- 
«ent linear-elastic behavior. However, for nonlinear large-displacement 
elastic or elastic-plastic behavior, numerical experimentation [17, for 
example] has confirmed that a smaller At is usually required to avoid this 
Instability. As a rough guide one may try At < 0.8 (2/«) as an 5mitial 
selection; if this value is not small enough, the calculation will blow up 
(overflow) before many time steps have elapsed - this behavior will be 
readily apparent. In such a case one must choose a smaller time increment At. 

Alternately, one could employ an implicit finite-difference time operator 
such as that of Houbolt [3,17,52] or Park [53], for example. The solution 
procedure is similar except that nonlinear large deflection and plasticity 
effects are handled either by iteration or extrapolation. Further discussion 
of these alternate procedures is beyond the intended scope of this review. 

Finally, this review of the solution procedure for the governing 
equations for large-deflection, elastic-plastic, 2-D structural response has 
has dealt with the case in which the structure has been subjected to 
prescribed transient externally-applied loading. For cases involving 
2-D structural response produced by fragment impact, the reader is invited 
to read the (more lengthy) descriptions-given in (1) Ref. 4 concerning 
si.-»gle-layer structures and the CIVM-JET 4B program and (2) Ref. 5 concerning 
multilayer hard-bonded Bernoulli-Euler 2-D structures. 
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R(n) “ -Oi>/3n) 
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h(n) - h, (1 - J.) + h. . JL 


LOCAL SYSTEM 

5. n» C - COORDINATES 

v.'»»4».X - DISPLACEMENTS 

" ELEMENT GENERALIZED 
DISPLACEMQITS 
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u. 3w V 


i+1 


CARTESIAN REFERENCE 

Tf , S . ~ GLOBAL 

COORDINATES 
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COORDINATES 


3n 


3n R 


fig. A.l »C.^a^URi: for GEOMETRY, COORDINATES. AND DISPLACEMENTS OF 
A CURVED-BEAM FINITE ELEMENT 
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UNIAXIAL STRESS, 
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(c) Schematic of Loading, Unloading, and 
Reloading Paths 


FIG. A. 2 CONCLUDED 
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STRESS 



STRAIN, C 

(a) Elastic, Perfectly-Plastic Material 



STRAIN, € 

(b) Special Strain Hardening Material 


KIo. A. 3 SCatriAM* JF STRAIN-RATE DEPENDENT UNIAXIAL STKi SS-STRAiN vTJKVES 
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appendix b 


® «>« smiH co«.mo«s 

B»1 Ob’jeg’^’ftpoo 

The intention in this appendix is to ««— • 

displaconent and stress-strain i P«rtinent strain- 

transw rT ““ 

^ t^rnr 

T into account*- These 

in detail in Ref. 54. elations are developed 

B^2 Strain-Displacement 

Fi 9 . A.1 with ^ 

respectively) deforms according to the B ' ^ ‘*^®ctions 1, 2 , and 3 , 

J22d^ to account for deformation-induceTri"^”'" hypothesis but 
structure. That is, it is assumed that: thinning of the , 


(a) 

(b) 

(c) 


plane sections remain plane, 
normals to the midsurfaee /«,. 

before deformaUon remain normal ^^ZtZll 7 ” 
deformation, but surface after 

these normals may stretch or contraci- — 
approximation that in the olastl according to the imposed 

bAin, i.^res,^,T,Zs uT “ 

r:L“ ^ - 

trams are assumed to be small) . 

‘r/‘ ““■= «>“ • '‘•it. cuc.,1 

— s <r’ . 0, ^ ^ 

''idth (or i) direction (that is, y^ = 0) . i across-the- 

that the extensional strains * 3" ® show ( 54 ) 

locauon in the structure are JIvenX '' " 


t =C 


(B.: 
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IB. 2) 


n = y‘ 




r. = 


. 0-2 n) 
' 1 


(B.3) 


where 


C is the C-location of the -material" before deformation 
(identified by sv^ierscript -o?) 


« the membrane strain at =» 0 


(B.2a) 


(B.2b) 


“ "diange of curvature" 

The TOMtitles X T. as dsflMa b, E,,. ,.23^ ^ »spscti™i. 

Also, it c» te ,ho». IM) thst tt, tr». cbst,t of ctrvature t^orrod to 
the deformed midsurface is given by 


d9 

ds 


0 nr 


(B.4) 


.hat. 9 i, u» . 0 , 1 . Of rotation, and a i. th. coordinat. alon, the dafor»d 
tateenoa aurfaca. Thinnin, affacts ara alao includad in Eqs. B. 1 -B. 3 . 

It shoald ba notad that Bin. 9.1 tbron,h 9.3 asply for finita aaaalo. 
«.d ara danotad harainaftar. for oonvanianca, aa tba Tfpa r atrain-diaplaca- 
ment relations, however, the structure must be thin such that 1 + ^ r 1 
ahara R ia tba originai undafo™d radio, of oorvntata. 0 . th. otj.r bald’, 
it shouid b. anphaaitad that th. -a«n attain- .apraaaion for x* ^ 
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Eq. A. 23 applies to both strains and small rotations such that 

1 + sin 9 " 0, and cos 0 r 1 where 0 is the angle of rotation. 

B.3 -Assumed-Displacgnent Field • 

For the BemouUi*Exiler displacement behavior cited in Subsection B.2, 
the assianed-displacenent field for a curved beam or ring remains the same 
as that defined in Subsection A. 3.2. 

For convenient reference# the strain-displacement relations termed 
Types A# B# C# D# and E in Ref. 20 are (in the present notation) given by;* 



B.4 Finite-Element Properties 

The finite element property matrices of interest are [m] # {p} # [h] # and 
{f} as given by Eqs. A. 9a# A.9b# A. 9c# and A.9d# respectively, and/or by 
Eqs. A. 36a# A.38, A. 39# and A. 41# respectively. Note that each of these 
matrices is evalxiated by volume or surface integrals over the original 
undeformed volume or surface. These evaluations are straightforward except 
for {p} and (hj given by Eqs. A.9b and A. 9c# respectively# which involves 
the use of the second Piola-Kirchhoff stress — however# in finite- 
strain plasticity one vor)cs with the Kirchhoff stress T versus logarithmic 
strain e*. Thus, one can represent the mechanical stress-strain property 
data from uniaxiaJ. (direction r\ or "2“) static (superscript "s") tensile 
and/or con^ression data in terms of vs. £* since this stress -strain 
information is essentially perfectly antisyimaetric with x monotonically 
increasing (decreasing) as e* increases (decreases) # where 
* 

Note that Eq. B.2 may be viewed as a “modification" of the Type E relation 

given in Eq. B.5. 
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(B.6) 






(B.7) 


and 


u 


^ » engineering stress of a uniaxial static test 

o 

specimen; P is the applied load and A is the 

o 

pre-test cross-sectional area of the specimen 

the measured axial (n-direction) relative 

elongation (also known as the engineering 

strain) of the uniaxial test specimen 

change in gage length 
original gage length 

output %dixch strain gages or extensometers 
can^provide 


Necking effects, if any. should be deleted from the data. Next, one makes 
a piecewise-linear fit of the vs. e* data in terms of n + 1 segments 

defined by the coordinate k - 1,2 n; the coordinates of 

these segments are used in the mechanical sublayer material model to repre- 
sent this behavior. According to the mechanical sublayer Bodel, the material 
is envisioned as consisting, at any point in- the- material, of n equally- 
strained sublayers of elas^c, perfectly-plastic material with each sublayer 
having the same elastic modulus E but an appropriately different yield 

stress. For exaBi>le, the static yield stress of the kth sublayer is 
given by 

(ti - E 

Hence, at any given state of strain 

-u 


(B.8) 

the associated stress is given by 
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(B.9) 


Hi 


where the weighting factor Aj^ for the kth mechanical sublayer may be readily 
ooxifirmed to be: 






(B.IO) 


where 


= E 


for k 




for k « 2,3,...n 


(B.ll) 


^n+l ^ 

The me<*anical sublayer model is very convenient to analyze problems with 
general loading paths — including loadin.. unloading, reloading, and cyclic 
loading; also, it approximates the Bauschinger effect reasonably well. 
Further# strain-rate effects can be accosaaodated by treating the sublayer 
yield stress as being strain-rate dependent according to: 


(c-i J 


(B.12) 


- I 


where (T ) is the rate-dependent yield stress# D = y « rate-of-defornation 
u Jc u ^ 

axial caagonent in a uniaxial test# and d and p are material constants which 
could be evaluated from experiments for each mechanical sublayer. Hence# 
the stress at any given condition of strain t* and strain rate may be 
written as 


X 


= L 

4=1 




(B.13) 


For a curved beam which experiences: (1) significant circumferential 

(q •2" direction) strain and through- the- thiclcness (C or **3" direction) 
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strain but negUgible^lateral (C or "1- direction) strain and (2) non- 


one may ej^ress the stress 


negligible stresses only in the n-direction, 

sublayer in terms of the circumferential component 
^2 the rete of defonmtxon teosojr by 




(B.14) 


if is entirely elastic (i.e., if d| - (D^)*) and hence I(xf)j2 < , j 

A Jt — u 


However, 


if 


and if 


(% > (-‘A . . fe), 


(B.15) 


(B.16) 


Accordingly, the total stress x^ at a given deformation and rate of deforma- 
tion is given by 


= Z (rn 


A-t 


(B.17) 


However, for the evaluation of the element property matrices {p} and (hj , 
one needs to evali^te and use the second Piola-Kirchhoff stress or sf 
and the associated Lagrangian or Green strain y or y?. Thus, one makes 
use of the following relations (54J: 


rl = s‘ 0 ^ Z n) 


(B.18a) 




- S^(i-fZ.r^) +z^ y 


z 

2 


ce = 


0 


(B.lSb) 


(B.lSc) 
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In the timewise step-by-step solution process and associated computer 

2 *2 

program, one has available Y 2 ^ 2 * Thus, since 


and 




(B.19a) 


(B.19b) 


it follows that 


Ko^^K) = Y- K(s:u(+zjr:)] 


= 0^zrl)Z A. (Si), 


Hence, 


^=i 


n 


s: = 


1. \ 


A^i 


(B.19c) 


(B.20) 


Next, one can rewrite Eq. B.14 in tezins of second Piola-Kirchhoff stress 
and Lagrangian (Green) str 2 un information from Eqs. B.18b and B.lSc to obtain 
for the kth sublayer: 




n 


or 


(i + zy') 


/CA _ 


One may integrate this differential expression by means of the trapezoidal 
rule, for exan 5 >le, from time instant i-1 to an incrementally close instant 
i to c^tain: 
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where 


AS - sy'[(i+zr^)-Ai-]f 

^ [0+2^‘)'-0+2y%ri-z(Ar)’] 

^ ( s:)> 


r = 




^A-S^ 

4y= x‘ 


(B.22) 


(B.22a) 

(B.22b) 

(B.22c) 

(B.22fi) 


Of . a, ^ clr=™rf.r^..tl.l 

stress (Sj)^ = s St U«e tj 4t siqr ,1^. tircumferentlsl station and any 

,i»on depthwise station in a onr^d-Oea. eies„n,t iouows; such infotnation 

IS needed since the volnE» inteprals detinin, {p) and (hi are evaloated by 

span.ase/depthvise nnnerioai intepration (in parUcnlar, Oanssian ^dra- 

tnre) . :^icaily, one begins at the station in guestion by tooein, the 

stress (S^ ) and the strain (yj' ) of the kth snblayer at ti™ t. , and 

the strain increment 4y at that station from time t. to time t ' 

1-1 


begin, one assnnms as a trial (superscript T) that the trial str^s at time 
tj= nay be evaluated by follguing an inorementally-elastic path. 

at}] 

[O * 2 rf-(hzr‘)Ar+ 2 CatT] 

A check is then performed as follous to see what the correct value of 
must be: k 


s; = s;v 


(B.23) 
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9 


“ -Cxl ^ Sj0^^r‘) ^ (xl 

t 

“ 5j(i + zr‘)> Cr.)^ 5 

S‘ = ^ X)ji 

(i +zr‘). 

« sj (i + zr^) , 

<ri CrJ^ 

^•-JTTPT 


(B.24) 


(B.25) 


(B.26) 


This procedure is applied to all sublayers at the station in question. 
Hence, the stress at time t^ is given by 


S‘. (sff = t A, Sj 


(B.27) 


4^ = / 


In the above process, strain-rate effects have been taken into account 


by using: 


(x)^ [ 


I + 


Dl IT 


] 


(B.28) 


where D_ is given by Eq. B.18c. Thus, 


N.‘(V. 


I + 




J(i+ 2 ir‘) 


± 1 
t> 


(B.29) 
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